Research on Fluid–Solid Coupling Mechanism around Openhole Wellbore under Transient Seepage Conditions

: Hydraulic fracturing is one of the most important enhanced oil recovery technologies currently used to develop unconventional oil and gas reservoirs. During hydraulic fracture initiation, ﬂuid seeps into the reservoir rocks surrounding the wellbore, inducing rock deformation and changes in the stress ﬁeld. Analyzing the ﬂuid–solid coupling mechanism around the wellbore is crucial to the construction design of fracturing technologies such as pulse fracturing and supercritical carbon dioxide fracturing. In this study, a new transient ﬂuid–solid coupling model, capable of simulating the pore pressure ﬁeld and effective stress ﬁeld around the wellbore, was established based on the Biot consolidation theory combined with the ﬁnite difference method. The numerical results are in excellent agreement with the analytical solutions, indicating the reliability of the model and the stability of the computational approach. Using this model, the inﬂuence of seepage parameters and reservoir properties on the ﬂuid–solid coupling around the open-hole wellbore was investigated. The simulation results demonstrate that, during wellbore pressurization, signiﬁcant changes occur in the pore pressure ﬁeld and effective stress ﬁeld near the wellbore. The ﬂuid–solid coupling effect around the wellbore returns to its initial state when the distance exceeds four times the radius away from the wellbore. As the ﬂuid viscosity and wellbore pressurization rate decrease, the pore pressure ﬁeld and effective circumferential stress (ECS) ﬁeld around the wellbore become stronger. Adjusting the ﬂuid viscosity and wellbore pressurization rate can control the effect of seepage forces on the rock skeleton during wellbore ﬂuid injection. For the same injection conditions, rocks with q higher Young’s modulus and Poisson’s ratio exhibit stronger pore pressure ﬁelds and ECS ﬁelds near the wellbore. This model furnishes a dependable numerical framework for examining the ﬂuid–solid coupling mechanism surrounding the open-hole wellbore in the initiation phase of hydraulic fractures.


Introduction
Hydraulic fracturing, also known as fracking, is a method used to extract natural gas or oil from underground reservoirs by injecting high-pressure water and chemicals into the rock formations to create and widen fractures, thus facilitating the release of energy resources.Hydraulic fracturing technology plays a pivotal role in the extraction of unconventional hydrocarbon resources, revolutionizing the energy landscape [1,2].The significance of this technology stems from its ability to enhance well productivity and unlock vast reservoir potential.In hydraulic fracturing, a thorough understanding of fracture initiation and propagation is essential for designing and optimizing fracturing operations effectively.Analyzing the nuanced processes involved in crack initiation and propagation is crucial in tailoring hydraulic fracturing strategies, ultimately influencing the success and efficiency of well stimulation endeavors.Therefore, a comprehensive investigation into the intricacies of fracture nucleation and propagation is essential for advancing the state of the art in hydraulic fracturing construction design [3][4][5][6][7].
Drawing upon the linear elastic theory and Terzaghi's effective stress theory, numerous scholars have previously conducted in-depth analyses of the stress field around wellbores during the injection of working fluids.Notable contributions include the work of Haimson and Hubbert [8][9][10][11], who systematically analyzed the stress field around open-hole wellbores.They categorized the stress field into three components: one arising from the geostatic stress of the reservoir, another from wellbore pressure, and the third from fluid filtration.Their research suggests that fracturing occurs when the effective total stress around the wellbore reaches the tensile strength of the rock formation, with the wellbore pressure becoming the formation's fracturing pressure.This led to the derivation of widely used formulas for open-hole well formation fracturing pressure, known as the H-W formula and the H-F formula.Subsequent studies have built upon this foundation.For instance, Reisabadi et al. [12] introduced an analytical model that considers the effects of depletion, matrix shrinkage, and wellbore trajectory when estimating the stress distribution around the wellbore.Li et al. [13] developed a heat-fluid-solid coupling numerical model for horizontal wells in NGH reservoirs, taking into account hydrate decomposition.Ren and Cai [14] established a three-dimensional model to calculate temperature and stress distributions in reservoir rock when liquid nitrogen is injected into the wellbore.Raad et al. [15] investigated the dynamics of a CO 2 injection wellbore using a fully coupled non-isothermal transient multiphase wellbore-reservoir flow model.Qiu et al. [16] created a pseudo-3D hydraulic-mechanical coupling finite element model, considering the DM and anisotropy of shale in wellbores.This model allowed for a parametric analysis of the effects of elastic anisotropy, fracture volume fraction, fracture permeability, mass transfer coefficient, and fluid pressure on wellbore stress.Xiao et al. [17] proposed a semi-analytical model for coal wellbore stability during mixture gas injection.Additionally, Rim et al. [18] developed a new seepage model that employs the method of the separation of variables and the Bessel function theory to predict variations in the pore pressure and seepage force around a vertical wellbore during hydraulic fracturing.By integrating fluid dynamics and rock mechanics models, Hou et al. [19] developed a comprehensive data-driven workflow to predict pressure based on fracturing experience from adjacent wells.
Furthermore, by assuming the perforation channel as a small-scale open-hole wellbore, many scholars have analyzed the stress field around the perforation tunnel under the assumption of an impermeable reservoir, deriving analytical solutions for the formation fracturing pressure of deviated perforated wells [20][21][22].Wang et al. [6] simplified the perforation channel as an open-hole wellbore and analyzed the impact of fluid seepage on the effective stress field around the perforation tunnel, deriving the initiation pressure for fractures at arbitrary perforation angles.While existing research has provided initial insights into the coupling effect between fluid seepage and rock deformation around openhole wellbores, most of these studies are based on analytical models established using steady-state Darcy's flow theory.Fluid seepage into reservoir rocks during fluid injection is a transient process that depends on time.Therefore, it is necessary to establish transient fluid-solid coupling models to investigate the influence of factors such as fluid viscosity and wellbore pressurization rate on the fluid-solid coupling around the wellbore.
This study initially established a numerical simulation model for fluid-solid coupling, considering transient fluid seepage effects, based on the Biot consolidation theory and the finite difference method.Subsequently, the reliability and stability of the fluidsolid coupling model were examined using a pore pressure analytical solution model.Finally, this model was employed to comprehensively analyze the impact of seepageinfluencing parameters and reservoir properties on the fluid-solid coupling around the open-hole wellbore.

Fluid-Solid Coupling Model
To analyze the fluid-solid coupling effect around the open-hole wellbore under transient seepage conditions, this paper establishes a two-dimensional planar wellbore physical model.As shown in Figure 1, the wellbore has an inner radius of a and an outer radius of c, where r represents the distance from the wellbore axis, θ represents the angle, the water pressure at the inner boundary of the cylinder is P w , and the initial pore pressure at the outer boundary is P o .It is assumed that the rocks surrounding the open-hole wellbore are isotropic and homogeneous.Under the influence of the pressure difference between the inner and outer boundaries, fluid seeps radially from the wellbore into the pores of the reservoir rock.By establishing a polar coordinate system, this paper analyzes the effective radial stress σ' r and ECS σ' θ formed in the rock around the cylinder due to fluid action.(Note that in this paper, positive values are assigned to compressive stress, while negative values represent tensile stress).

Fluid-Solid Coupling Model
To analyze the fluid-solid coupling effect around the open-hole wellbore under transient seepage conditions, this paper establishes a two-dimensional planar wellbore physical model.As shown in Figure 1, the wellbore has an inner radius of a and an outer radius of c, where r represents the distance from the wellbore axis, θ represents the angle, the water pressure at the inner boundary of the cylinder is Pw, and the initial pore pressure at the outer boundary is Po.It is assumed that the rocks surrounding the open-hole wellbore are isotropic and homogeneous.Under the influence of the pressure difference between the inner and outer boundaries, fluid seeps radially from the wellbore into the pores of the reservoir rock.By establishing a polar coordinate system, this paper analyzes the effective radial stress σ'r and ECS σ'θ formed in the rock around the cylinder due to fluid action.(Note that in this paper, positive values are assigned to compressive stress, while negative values represent tensile stress.)For an open-hole wellbore, the stress equilibrium equation and stress boundary conditions satisfied by the surrounding rock when fluid seepage enters the rock pores are as follows [23]: Here, α is the Biot effective stress coefficient and P represents the pore pressure.As shown in Equation (1), the effective stress field formed in the surrounding reservoir rock due to fluid action around the wellbore can be decomposed into the superposition of stress fields generated by both volume forces and surface forces, including: (1) internal and external boundary surface forces Pw − αPw and Po − αPo; and (2) the seepage force term αdP/dr.
To establish a finite difference format, the first step is to create a stress function Φ to characterize the radial stress σ'r and circumferential stress σ'θ: For an open-hole wellbore, the stress equilibrium equation and stress boundary conditions satisfied by the surrounding rock when fluid seepage enters the rock pores are as follows [23]: Here, α is the Biot effective stress coefficient and P represents the pore pressure.As shown in Equation (1), the effective stress field formed in the surrounding reservoir rock due to fluid action around the wellbore can be decomposed into the superposition of stress fields generated by both volume forces and surface forces, including: (1) internal and external boundary surface forces P w − αP w and P o − αP o ; and (2) the seepage force term αdP/dr.
To establish a finite difference format, the first step is to create a stress function Φ to characterize the radial stress σ' r and circumferential stress σ' θ : Due to the axisymmetric nature of the stress field formed by fluid action around the wellbore, the geometric and physical equations satisfied by the rock surrounding the wellbore under fluid influence can be obtained using plane strain conditions: Combining Equation ( 2) with (3) and substituting them into Equation (1), we can derive the compatibility equation and boundary conditions satisfied by the stress function as follows: To solve the compatibility equation, it is necessary to calculate the distribution of pore pressure around the wellbore under transient seepage conditions.In the Biot poroelastic consolidation theory [24,25], the governing equation for pore pressure propagation under the combined influence of solid deformation and fluid is as follows: In the equation, the parameters are defined as follows: Q is the Biot modulus, when α = 1 and only the seepage force is considered, Q = 1/φc f , where φ is the porosity of the rock, and c f is the fluid compressibility in Pa −1 .G is the shear modulus defined as , where E is the Young's modulus in Pa and ν is the drained Poisson's ratio.ν u represents the undrained Poisson's ratio.t is time measured in seconds.κ is the permeability coefficient defined as κ = k/µ, where k is the permeability in m 2 and µ is the fluid viscosity in Pa•s.ε kk represents volumetric expansion.
To facilitate a convenient solution for the pore pressure propagation equation, Geertsma and Yew, among others [26,27], derived a conversion equation relating volumetric expansion (ε kk ) to pore pressure (P) based on the Biot consolidation theory: By utilizing Equation ( 6) to process Equation ( 5), we can obtain the pore pressure propagation equation that takes into account the influence of solid deformation: Processes 2024, 12, 412 5 of 16 The combination of the Compatibility Equation ( 4) and Pore Pressure Propagation Equation ( 7) forms the fundamental control equations for the transient fluid-solid coupling model around a bare wellbore under seepage conditions.By applying the finite difference method to discretize the above equations, we obtain: where h is the iterative step size, h = c − a/n; ∆t is the time step size; and (r i , t j ) represents the grid nodes.The adopted time scheme is backward differencing.Since the solution area is two-dimensional and axially symmetric, there is no need to determine the unit type, just mesh the solution area.To solve Equation ( 8), we meshed the region within a certain radius around the open-hole wellbore, as illustrated in Figure 2.
By utilizing Equation ( 6) to process Equation ( 5), we can obtain the pore pressure propagation equation that takes into account the influence of solid deformation: The combination of the Compatibility Equation ( 4) and Pore Pressure Propagation Equation ( 7) forms the fundamental control equations for the transient fluid-solid coupling model around a bare wellbore under seepage conditions.By applying the finite difference method to discretize the above equations, we obtain: where h is the iterative step size, h = c−a/n; Δt is the time step size; and (ri, tj) represents the grid nodes.The adopted time scheme is backward differencing.Since the solution area is two-dimensional and axially symmetric, there is no need to determine the unit type, just mesh the solution area.To solve Equation ( 8), we meshed the region within a certain radius around the open-hole wellbore, as illustrated in Figure 2. Taking hydraulic fracturing as an example, according to the fracturing construction curve, the wellbore pressure can be considered to increase at a constant pressurization rate [28].Meanwhile, at the far-field boundary, the pore pressure can be assumed to be Taking hydraulic fracturing as an example, according to the fracturing construction curve, the wellbore pressure can be considered to increase at a constant pressurization rate [28].Meanwhile, at the far-field boundary, the pore pressure can be assumed to be the initial pore pressure.Therefore, the boundary conditions for the pore pressure propagation equation are as follows: where C represents the wellbore's pressurization rate, MPa/s o .
As indicated by the boundary conditions specified in Equation ( 9), the wellbore pressure is directly influenced by both the wellbore pressurization rate and time.Meanwhile, at the far-field outer boundary, the pore pressure remains constant and equal to the initial pore pressure.
According to the governing Equation ( 8) and boundary conditions 9, the tridiagonal matrix required to solve the stress function and pore pressure can be established.The tridiagonal matrix AX = B for pore pressure is: where The tridiagonal matrix that the stress function satisfies is CY = D: Under the boundary conditions specified in Equation ( 9), we simultaneously solve the aforementioned matrices through programming to obtain the pore pressure field and stress function Subsequently, this stress function can be utilized to solve and acquire the stress field resulting from the fluid interactions around the wellbore under transient fluidsolid coupling conditions.During hydraulic fracturing, the effective circumferential stress directly affects the fracture initiation pressure [29], so subsequent chapters will only focus on the circumferential stress around the wellbore.

Model Validation 2.2.1. Pore Pressure Field Verification
When considering the absence of rock matrix particle compressibility (Biot coefficient α equals 1) and neglecting the influence of solid deformation on pore pressure, the pore bnpressure propagation Equation ( 7) can be simplified to: where χ is the propagation coefficient, m 2 /s; this study assumes c f = 4.5 × 10 −10 Pa −1 .
Despite not accounting for the influence of solid deformation on pore pressure, Equation ( 14) remains a transient pore pressure propagation equation dependent on time.Based on Equation ( 14), Carslaw and Jaeger [30] derived an analytical solution for the radial distribution of pore pressure around a hollow cylindrical body under a constant pressurization rate: where: J 0 and Y 0 are the first and second kinds of zero-order Bessel functions.By adjusting the parameter β in the finite difference Equation ( 8) and comparing it with the analytical solution for the transient pore pressure distribution mentioned above, we can verify the accuracy of the difference format established in this work.
Figure 3 compares the analytical solution with the numerical solution of pore pressure distribution around the wellbore under a constant pressurization rate.It is evident from Figure 3 that, starting from the initial moment, the pore pressure distribution near the wellbore undergoes significant variations due to wellbore pressurization.As one moves away from the wellbore, the distribution gradually approaches the initial pore pressure before pressurization.At different time instances, the analytical solution for transient pore pressure and the numerical solution presented in this study demonstrate a close agreement, with the relative error limited to less than 0.1%.This validation confirms the accuracy of the transient fluid-solid coupling model established in this paper.

Effective Stress Field Verification
There is no analytical solution available for the stress field around the wellbore under transient seepage conditions.However, an analytical solution does exist for the stress field around the wellbore in the case of impermeable reservoirs.When the rock is impermeable, At different time instances, the analytical solution for transient pore pressure and the numerical solution presented in this study demonstrate a close agreement, with the relative error limited to less than 0.1%.This validation confirms the accuracy of the transient fluid-solid coupling model established in this paper.

Effective Stress Field Verification
There is no analytical solution available for the stress field around the wellbore under transient seepage conditions.However, an analytical solution does exist for the stress field around the wellbore in the case of impermeable reservoirs.When the rock is impermeable, the ECS σ' θ induced by the fluid pressures (P w and P o ) at the inner and outer boundaries of the wellbore can be calculated using the Lame formula and Terzaghi's principle of effective stress [31][32][33].
When c >> a, Equation ( 17) can be simplified to a 2 r 2 (P o − P w ).By varying the reservoir's permeability, we employed the finite difference scheme detailed in Section 2.1 to compute the effective circumferential stress resulting from the fluid flow around the wellbore during transient seepage.We investigated the variations in the transient stress field surrounding the wellbore as the reservoir's permeability gradually decreased, eventually reaching impermeability.As depicted in Figure 4, the results obtained from the numerical solution are compared with those derived from the analytical solution provided by Equation (17).In the case of an impermeable reservoir, the seepage force term in Equation ( 1) is negligible, and the fluid primarily exerts a surface force on the rocks surrounding the wellbore.As depicted in Figure 4, for impermeable reservoirs, the effective circumferential stress resulting from fluid interaction around the wellbore is tensile stress.However, when fluid seepage enters the reservoir, the effective circumferential tensile stress induced by the fluid around the wellbore, due to the seepage force, increases significantly.A greater effective circumferential tensile stress increases the susceptibility of the well wall to tensile failure.This demonstrates that the seepage force significantly increases the possibility of tensile failure of the wellbore wall.Under identical injection conditions, an increasing reservoir fluid permeability amplifies the seepage effect, leading to a higher effective circumferential tensile stress around the wellbore.
As the reservoir's permeability gradually diminishes, the numerical solution begins to converge toward the analytical solution.When the reservoir's permeability approaches an impermeable state, the analytical and numerical solution curves essentially overlap.This confirms the stability and accuracy of the finite difference scheme presented in this article for computing the effective circumferential stress field.In the case of an impermeable reservoir, the seepage force term in Equation ( 1) is negligible, and the fluid primarily exerts a surface force on the rocks surrounding the wellbore.As depicted in Figure 4, for impermeable reservoirs, the effective circumferential stress resulting from fluid interaction around the wellbore is tensile stress.However, when fluid seepage enters the reservoir, the effective circumferential tensile stress induced by the fluid around the wellbore, due to the seepage force, increases significantly.A greater effective circumferential tensile stress increases the susceptibility of the well wall to tensile failure.This demonstrates that the seepage force significantly increases the possibility of tensile failure of the wellbore wall.Under identical injection conditions, an increasing reservoir fluid permeability amplifies the seepage effect, leading to a higher effective circumferential tensile stress around the wellbore.
As the reservoir's permeability gradually diminishes, the numerical solution begins to converge toward the analytical solution.When the reservoir's permeability approaches an impermeable state, the analytical and numerical solution curves essentially overlap.This confirms the stability and accuracy of the finite difference scheme presented in this article for computing the effective circumferential stress field.

Evolution Characteristics of Pore Pressure Field and Effective Stress Field
In this section, α is set to 1, resulting in the disappearance of the surface force term in Equation ( 1), with fluid exerting only seepage force on the rock framework.By combining the observations from Figures 5 and 6, it becomes evident that significant changes occur in the pore pressure field around the wellbore and the Effective Circumferential Stress (ECS) field formed by the action of seepage forces during the stable pressurization process of the wellbore.The pore pressure and ECS near the wellbore exhibit larger variations over time, while moving away from the wellbore results in a rapid decrease in pore pressure towards the initial pore pressure state.Under transient fluidsolid coupling conditions, the ECS formed by the action of seepage forces is maximized near the wellbore and gradually diminishes towards zero with an increasing distance from the wellbore.The magnitude of changes in the ECS due to seepage forces is more pronounced in the vicinity of the wellbore.By combining the observations from Figures 5 and 6, it becomes evident that significant changes occur in the pore pressure field around the wellbore and the Effective Circumferential Stress (ECS) field formed by the action of seepage forces during the stable pressurization process of the wellbore.The pore pressure and ECS near the wellbore exhibit larger variations over time, while moving away from the wellbore results in a rapid decrease in pore pressure towards the initial pore pressure state.Under transient fluid-solid coupling conditions, the ECS formed by the action of seepage forces is maximized near the wellbore and gradually diminishes towards zero with an increasing distance from the wellbore.The magnitude of changes in the ECS due to seepage forces is more pronounced in the vicinity of the wellbore.
pressurization process of the wellbore.The pore pressure and ECS near the wellbore exhibit larger variations over time, while moving away from the wellbore results in a rapid decrease in pore pressure towards the initial pore pressure state.Under transient fluidsolid coupling conditions, the ECS formed by the action of seepage forces is maximized near the wellbore and gradually diminishes towards zero with an increasing distance from the wellbore.The magnitude of changes in the ECS due to seepage forces is more pronounced in the vicinity of the wellbore.

Influence of Seepage Parameters
Using the established transient fluid-solid coupling model, an analysis of two key parameters, fluid viscosity and wellbore pressurization rate, was conducted.Figure 7

Influence of Seepage Parameters
Using the established transient coupling model, an analysis of two key parameters, fluid viscosity and wellbore pressurization rate, was conducted.Figure 7 presents pseudo-color plots of the pore pressure field around the wellbore and the corresponding ECS field formed by fluid action at the same time step but with varying fluid viscosities.Figure 8 displays pseudo-color plots of the pore pressure field and ECS field when the wellbore pressure is pressurized to the same level using different pressurization rates (α = 1; ν = 0.25, E = 5.0 × 10 4 MPa, K = 0.1 mD, φ = 0.1, and P 0 = 5 MPa).In Figure 7, it is evident that the fluid viscosity has a significant impact on the pore pressure field around the wellbore and the ECS field formed by seepage forces under transient conditions.At the same time step, when the fluid viscosity is high, the pore pressure rapidly decreases to the initial pore pressure state around the wellbore as fluid seeps into the reservoir.Conversely, as the fluid viscosity decreases, the pore pressure distribution In Figure 7, it is evident that the fluid viscosity has a significant impact on the pore pressure field around the wellbore and the ECS field formed by seepage forces under transient conditions.At the same time step, when the fluid viscosity is high, the pore pressure rapidly decreases to the initial pore pressure state around the wellbore as fluid seeps into the reservoir.Conversely, as the fluid viscosity decreases, the pore pressure distribution around the wellbore gradually becomes smoother.Additionally, with a decreasing fluid viscosity, the ECS field formed by seepage forces around the wellbore significantly strengthens, resulting in a notable increase in effective circumferential tensile stress values.Therefore, for low-permeability reservoirs, reducing fluid viscosity can significantly enhance the action of seepage forces around the wellbore.
Figure 8 demonstrates that, when pressurized to the same pressure level, the pressurization rate significantly affects the pore pressure field and ECS field.A smaller wellbore pressurization rate results in a longer time required to reach the same pressure level, but leads to higher pore pressure values around the wellbore.Conversely, a higher pressurization rate leads to a faster decline in the pore pressure around the wellbore.Additionally, as the pressurization rate decreases, the ECS values formed by seepage forces around the wellbore gradually increase.Therefore, a lower pressurization rate results in higher ECS values formed by seepage forces around the wellbore when pressurized to the same pressure level.Consequently, for reservoirs with a lower permeability, it may be advisable to consider lower fluid injection rates to enhance the action of seepage forces.

Influence of Reservoir Parameters
Using the transient fluid-solid coupling model established in this study, we can analyze the effects of key reservoir parameters, such as reservoir permeability and rock mechanics, on the pore pressure field around the wellbore and the ECS field formed by seepage forces (α = 1, µ = 1 mPa•s, φ = 0.1, P 0 = 5 MPa, and C = 0.25 MPa/s).
Figure 9 illustrates pseudo-color plots of the pore pressure field and ECS field under different reservoir permeability conditions, while Figure 10 shows the corresponding variations in the pore pressure and ECS along the radial direction.Combining the observations from Figures 8 and 9, it becomes apparent that, under low-viscosity fluid conditions, decreasing the reservoir permeability amplifies the magnitude of pore pressure variations.Near the wellbore, a lower reservoir permeability leads to a faster decay in pore pressure.
In the simulation case depicted in Figure 10, with low-viscosity fluid, once the reservoir permeability exceeds 1 mD, the pore pressure distribution tends toward a steady-state linear Darcy seepage distribution, and the influence of permeability on the pore pressure field and ECS field gradually diminishes.Consequently, under low-viscosity fluid conditions, a higher reservoir permeability results in elevated pore pressure values around the wellbore and increased ECS values formed by seepage forces.
Figure 9 illustrates pseudo-color plots of the pore pressure field and ECS field under different reservoir permeability conditions, while Figure 10 shows the corresponding variations in the pore pressure and ECS along the radial direction.Combining the observations from Figures 8 and 9, it becomes apparent that, under low-viscosity fluid conditions, decreasing the reservoir permeability amplifies the magnitude of pore pressure variations.Near the wellbore, a lower reservoir permeability leads to a faster decay in pore pressure.In the simulation case depicted in Figure 10, with low-viscosity fluid, once the reservoir permeability exceeds 1 mD, the pore pressure distribution tends toward a steady-state linear Darcy seepage distribution, and the influence of permeability on the pore pressure field and ECS field gradually diminishes.Consequently, under low-viscosity fluid conditions, a higher reservoir permeability results in elevated pore pressure values around the wellbore and increased ECS values formed by seepage forces.Combining the observations from Figures 11-14 reveals that, under fluid-solid coupling conditions, the elastic deformation of rock elements significantly affects the pore pressure field around the wellbore and the ECS field formed by seepage forces.Higher values of Young's modulus and Poisson's ratio result in elevated pore pressure values around the wellbore, corresponding to higher ECS values formed by seepage forces.Overall, Young's modulus (E) has a slightly greater influence on the pore pressure field and the stress field formed by seepage forces than Poisson's ratio (ν).However, compared to seepage parameters, rock mechanics parameters have a much weaker impact on the pore pressure field and the ECS field formed by seepage forces around the wellbore.In this paper, there is still room for improvement.Future research will delve into the analysis of well wall behavior, specifically in terms of tensile or shear failure, induced by fluid injection.This analysis will be based on the fluid-solid coupling model presented here, combined with rock failure criteria.Combining the observations from Figures 11-14 reveals that, under fluid-solid coupling conditions, the elastic deformation of rock elements significantly affects the pore pressure field around the wellbore and the ECS field formed by seepage forces.Higher values of Young's modulus and Poisson's ratio result in elevated pore pressure values around the wellbore, corresponding to higher ECS values formed by seepage forces.Overall, Young's modulus (E) has a slightly greater influence on the pore pressure field and the stress field formed by seepage forces than Poisson's ratio (ν).However, compared to seepage parameters, rock mechanics parameters have a much weaker impact on the pore pressure field and the ECS field formed by seepage forces around the wellbore.In this paper, there is still room for improvement.Future research will delve into the analysis of well wall behavior, specifically in terms of tensile or shear failure, induced by fluid injection.This analysis will be based on the fluid-solid coupling model presented here, combined with rock failure criteria.

Conclusions
Based on the Biot consolidation theory and finite difference method, this paper established a transient fluid-solid coupling model that can analyze the stress field and pore pressure field around an open-hole wellbore during the initiation of hydraulic fractures.(1) During the stable pressurization of the wellbore, significant changes were observed in the pore pressure field and effective stress field near the wellbore.As the distance from the wellbore increased to more than four times the radius, the fluid-solid coupling effect around the wellbore substantially diminished, returning to its initial state.
(2) With a decrease in fluid viscosity and wellbore pressurization rate, the pore pressure and ECS around the wellbore increased.Adjusting the fluid viscosity and wellbore pressurization rate allowed control over the action of seepage forces during wellbore fluid injection.
(3) Under conditions of low-viscosity fluid injection, a higher reservoir permeability yielded larger pore pressure values around the wellbore, accompanied by increased ECS values generated by seepage forces.
(4) When the Young's modulus and Poisson's ratio of the rock surrounding the wellbore were large, the pore pressure field and ECS were stronger.Notably, the Young's modulus exhibited a more pronounced impact on the fluid-solid coupling effect around the wellbore compared to Poisson's ratio.

17 Figure 5 .
Figure 5. Pseudo-color maps of pore pressure field and ECS field at different time steps.

Figure 5 .
Figure 5. Pseudo-color maps of pore pressure field and ECS field at different time steps.

Figure 6 .
Figure 6.Variation in radial pore pressure and ECS with radius at different time steps.

Figure 6 .
Figure 6.Variation in radial pore pressure and ECS with radius at different time steps.

Figure 7 .
Figure 7. Pseudo-color maps of pore pressure field and stress field under different fluid viscosities (C = 0.25 MPa/s).

Figure 7 .
Figure 7. Pseudo-color maps of pore pressure field and stress field under different fluid viscosities (C = 0.25 MPa/s).

Figure 7 .
Figure 7. Pseudo-color maps of pore pressure field and stress field under different fluid viscosities (C = 0.25 MPa/s).

Figure 8 .
Figure 8. Pseudo-color maps of pore pressure field and ECS field under different pressurization rates (µ = 10 mPa•s).

Figure 8 .
Figure 8. Pseudo-color maps of pore pressure field and ECS field under different pressurization rates (µ = 10 mPa•s).

Figure 10 .
Figure 10.Radial pore pressure and ECS change with radius under different reservoir permeability conditions.Figures 11 and 12 present pseudo-color plots illustrating the pore pressure field around the wellbore and the ECS field formed by seepage forces under varying conditions of Young's modulus and Poisson's ratio.Correspondingly, Figures 13 and 14 provide radial profiles of the pore pressure and ECS along the radius.

Figure 10 .
Figure 10.Radial pore pressure and ECS change with radius under different reservoir permeability conditions.Figures 11 and 12 present pseudo-color plots illustrating the pore pressure field around the wellbore and the ECS field formed by seepage forces under varying conditions of Young's modulus and Poisson's ratio.Correspondingly, Figures 13 and 14 provide radial profiles of the pore pressure and ECS along the radius.

Figure 12 .
Figure 12.Variation in pore pressure and ECS in the radial direction with radius under different Young's modulus conditions.Figure 12. Variation in pore pressure and ECS in the radial direction with radius under different Young's modulus conditions.

Figure 12 .
Figure 12.Variation in pore pressure and ECS in the radial direction with radius under different Young's modulus conditions.

Figure 14 .
Figure 14.Variation in pore pressure and ECS in the radial direction with radius at different Poisson's ratios.

Figure 14 .
Figure 14.Variation in pore pressure and ECS in the radial direction with radius at different Poisson's ratios.
Based on the Biot consolidation theory and finite difference method, this paper established a transient fluid-solid coupling model that can analyze the stress field and pore pressure field around an open-hole wellbore during the initiation of hydraulic fractures.The transient analytical solution of the pore pressure field was used to verify the accuracy of the model and the stability of the calculation method.Using this model, this paper analyzed the effects of reservoir parameters and seepage parameters on open-hole wellbore fluid-solid coupling under different injection conditions.The simulation results show that: