Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Earth Sci., 09 September 2025

Sec. Georeservoirs

Volume 13 - 2025 | https://doi.org/10.3389/feart.2025.1659387

This article is part of the Research TopicIntelligent Artificial Lift and Multiphase Flow in the Wellbore in the Oil and Gas Production SystemsView all 3 articles

The research on intelligent evaluation methods for non-sealing faults in reservoirs based on well test analysis

Yiyi YangYiyi Yang1Boying Li
Boying Li2*Xudong WangXudong Wang1Mingming QiMingming Qi1Di ZhaiDi Zhai1
  • 1CNOOC International Limited, Beijing, China
  • 2College of Petroleum Engineering, China University of Petroleum, Beijing, China

In practical reservoirs, there exist a large number of non-sealing faults with flow capacity, which significantly influence oilfield development patterns and wellbore pressure transient behavior. The quantitative evaluation of fault sealing capacity is of great importance for characterizing remaining oil distribution and reconstructing flow fields. However, existing fault analysis methods are primarily qualitative, with limitations in the quantitative characterization of fault sealing. Traditional numerical well-test interpretation models do not account for fluid flow within faults, leading to significant deviations in well-test data interpretation, increased model-solving difficulties, and challenges in achieving quantitative analysis of reservoir sealing. Therefore, based on the fault, fluid, and reservoir property characteristics of the X reservoir, this study establishes a composite reservoir well-test interpretation mathematical model considering skin effects and solves the model using the Boltzmann transformation. By applying the “partial” mirror superposition principle, the dynamic response characteristics of typical curves under different fault boundary conditions are analyzed, and a quantitative sealing evaluation method suitable for non-sealing faults is developed. Furthermore, by integrating XGBoost multi-output regression and PSO algorithms, an intelligent hybrid inversion framework for identifying non-sealing faults in composite reservoirs is constructed: the XGBoost model predicts initial fault characteristic parameters, while the PSO algorithm performs global optimization to refine XGBoost parameters, ultimately inverting the fault connectivity coefficient (CfD) and effective connected thickness (hb). The results indicate the presence of a non-sealing fault F1 between Well B30Y and Well B1A, with inverted values of CfD=0.73 and hb=20.44, demonstrating strong fault connectivity. Additionally, the fitting trend of bottom-hole flowing pressure during shut-in periods in both wells verifies the validity and stability of the proposed model. The method presented in this study enables rapid, quantitative, and precise evaluation of non-sealing fault closure, providing robust technical support for subsequent remaining oil potential exploitation and development strategy optimization.

1 Introduction

Currently, research on sealing faults is relatively mature both domestically and internationally, whereas studies on non-sealing faults with certain conductivity are comparatively limited. In actual reservoirs, the majority of faults possess a certain degree of conductivity, necessitating consideration of factors such as fault sealing capacity, conductivity, and connectivity. Consequently, for well test data from composite reservoirs with non-sealing faults, the well test interpretation models designed for sealing faults are not applicable, and the results obtained from such interpretations are often unreasonable (Liu et al., 2016; Zeng et al., 2017; Xing et al., 2023; Ma, 2025).

The theories of mirror image and potential function is a classical theoretical method employed to address reservoir engineering challenges (Dake, 1978; Liu et al., 2017). A quintessential application of the mirror image in reservoir engineering is the determination of the distance from a test well to nearby fault utilizing pressure drawdown or buildup well test data. However, the traditional mirror image method has only investigated the scenario of sealed faults without fluid flow, where the slope of the well test curve in semi-logarithmic coordinates exhibits an integer or specific multiple change upon encountering a fault; the D. R. Horner semi-logarithmic curve method (Horner, 1951) is commonly used to identify sealed faults. Prasad (1975) extended the model of a single sealed fault to include multiple intersecting or parallel sealed faults.

However, in actual reservoirs, particularly for faults within the reservoir where the fault displacement is less than the thickness of the reservoir formation, the majority are non-sealing faults with a certain degree of conductivity (Fu et al., 2022; Zhang, 2025; Fu et al., 2024; Hu et al., 2024; Liu et al., 2024; Li et al., 2023; Wang and Zhang, 1993). These cannot be effectively interpreted using well test models designed for sealing faults. For non-sealing faults, Bixel et al. (1963) were the first to propose treating the fault as a non-sealing boundary, but the model they established only considered the scenario where the properties of the rock and fluid change abruptly in the plane; Kuchuk and Tarek (1997) further refined the solution method for this model. Building upon the research of Stewart et al. (1984) and Streltsova and McKinley (1984) utilized numerical simulation to investigate the impact of partially communicating faults on interference well testing. Yaxely (1987) derived an analytical model for an infinite homogeneous reservoir containing a partially communicating fault; Abdelaziz and Tiab (2004) expanded on the work of L. M. Yaxley to study the pressure transient behavior in a homogeneous reservoir with two intersecting fault boundaries. The existing physical models of non-sealing faults only consider the different fluid properties within the reservoir, while the flow of fluids within the fault itself has not been studied.

Meanwhile, most current intelligent well-test interpretation and inversion models fail to account for the presence of non-sealing faults in the reservoir. Li (2023) developed an automated well-test interpretation method based on convolutional neural networks (CNN), enabling the input of complete pressure derivative curves as variables and overcoming the limitations of traditional automatic matching approaches that rely on parameter optimization algorithms. Dong and Liao (2023) proposed an automated well-test interpretation method combining the Deep Deterministic Policy Gradient (DDPG) algorithm with a Long Short-Term Memory (LSTM) surrogate model, achieving efficient automatic curve matching under four distinct well-test model conditions. Mi et al. (2021) constructed a multi-modal stochastic analysis well-test interpretation model integrating model parameters and curves, employing ensemble Kalman filtering for well-test curve matching, thereby enhancing the interpretation of complex flow regimes and boundary effects. Li et al. (2020) introduced a CNN-based automated well-test interpretation method for radial composite reservoirs, realizing automatic parameter matching in such reservoirs.

Based on the fault, fluid, and reservoir properties of the X reservoir, this study establishes a mathematical model for well-test interpretation in composite reservoirs, which is solved using Boltzmann transformation. By applying the “partial” image superposition principle, we analyze the dynamic response characteristics of type curves under different fault boundary conditions. Furthermore, an inversion framework for non-sealing fault parameters in composite reservoirs is developed by integrating this model with eXtreme Gradient Boosting (XGBoost) multi-output regression and particle swarm optimization (PSO) algorithm, enabling efficient and accurate identification of fault characteristics.

2 Construction and intelligent solution of well test models for non-sealing faults

2.1 Establishment of well testing model for non-sealing faults

In an infinite homogeneous reservoir, the presence of a partially permeable fault with certain flow capacity is considered, where both sides of the fault exhibit partial conductivity. The physical model is illustrated in Figure 1. Well A is located at a distance d from fault F, and well A′ is a “mirror image” reflection well with respect to the fault plane F. Assuming the production rate of well A is qA, to simulate the scenario where fluid flows across the fault, the production rate of the “mirror image” well A′ is assumed to be: αqA. Here, α represents the degree of completeness of the mirror reflection. In an infinite homogeneous reservoir, the presence of a partially permeable fault with certain flow capacity is considered, where both sides of the fault exhibit partial conductivity. The physical model is illustrated in Figure 1. Well A is located at a distance d from fault F, and well A′ is a “mirror image” reflection well with respect to the fault plane F. Assuming the production rate of well A is qA, to simulate the scenario where fluid flows across the fault, the production rate of the “mirror image” well A′ is assumed to be: αqA. Here, α represents the degree of completeness of the mirror reflection. As shown in Figure 1, the reservoir is divided into an inner zone (Region 1) and an outer zone (Region 2), with distinct properties such as permeability, porosity, and fluid characteristics (i.e., a radial composite reservoir model). The outer zone contains a non-sealing boundary (i.e., a conductive fault), modeled as a linear conductive boundary. For simplification, this boundary is assumed to be located in the outer zone at a distance L from well A (representing a 1D radial flow boundary response).

Figure 1
Diagram illustrating fluid flow through a reservoir. A circle labeled

Figure 1. Schematic diagram of mirror image reflection in a homogeneous reservoir.

Similar to homogeneous and dual-porosity reservoirs, when interpreting well test data from composite oil and gas reservoirs with faults, the pressure derivative curve of the test data can become highly complex due to the presence of faults, affecting the mid- and late-time behavior of the curve. Therefore, the influence of faults on the characteristic well test model response of composite reservoirs can be analyzed by applying the partial image method and the superposition principle.

The model assumptions are as follows:

1. Both the inner and outer zones exhibit homogeneous, isothermal, single-phase fluid flow, with gravity and capillary forces neglected.

2. The well is located at the center of the inner zone (r = 0), producing at a constant rate (q). Wellbore storage effects (storage coefficient: C) and skin effects (skin factor: S) are considered.

3. The non-sealing boundary satisfies the conductive boundary condition: pressure continuity across the boundary; Flow rate proportional to the pressure difference between inner and outer zones (representing cross-boundary fluid transfer capacity)

Based on the assumptions, applying the theory of fluid flow mechanics, the model for well testing in a homogeneous reservoir is established:

Based on flow theory in porous media, the following mathematical model can be established:

1.Inner zone (0<r<rf) Equations:

1rrrk1p1r=ϕ1μ1Ct1p13.6t(1)

Where, k1: inner zone permeability, md; ϕ1: inner zone porosity; Ct1: inner zone total compressibility, MPa-1; p1: inner zone pressure, MPa; μ1: inner zone fluid viscosity, mPa·s.

2. Outer zone (r>rf) Equations:

1rrrk2P2r=ϕ2μ2Ct2P23.6t(2)

Where, k2: outer zone permeability, md; ϕ2: outer zone porosity; Ct2: outer zone total compressibility, MPa-1; p2: outer zone pressure, MPa; μ2: outer zone fluid viscosity, mPa·s.

3. Inner-outer zone interface conditions:

Pressure and flow rate continuity across the interface are governed by:

P1rf,t=P2rf,t(3)
k1μ1P1r|r=rf=k2μ2P2r|r=rf(4)

4. Non-sealing boundary conditions (outer zone boundary: r=d):

k2hμ2P2rr=d=kbh´bμ2P2Pi(5)

Where, kbh´bμ2: conductivity of the non-sealing boundary (CF); Pi: initial reservoir pressure, MPa.

5. Wellbore conditions and initial conditions

Accounting for wellbore storage and skin effects, the wellbore boundary condition is expressed as:

Cdpwdt=qB172.8πrwk1hp1r|r=rw(6)
Pwt=P1rw,t+172.8πqμ1Bk1hS(7)

Where, h: net pay thickness, m; B: formation volume factor; q: surface production rate, m3/d.

The initial reservoir condition is defined as:

P1r,0=P2r,0=Pi(8)

2.2 Solution methodology for non-sealing fault well testing model (laplace transform approach)

1. Dimensionless transformation

The dimensionless parameters are defined as shown in Equations 917.

PD1=k1h172.8πqμBPiP1(9)
PD2=k2h172.8πqμBPiP2(10)
tD=3.6k1ϕ1μCt1rw2t(11)
rD=rrw(12)
rfD=rfrw(13)
dD=drw(14)
CD=2πϕ1μCt1rw2C(15)
ω=k1ϕ2μ2Ct2k2ϕ1μ1Ct1(16)
CfD=kbhbrwk2hLb(17)

Where, CD: dimensionless storage coefficient; C: storage coefficient, m3/MPa; CfD: fault connectivity coefficient; kb: fault zone permeability, mD; hb: effective connected thickness of the fault: m; ω: storativity ratio between inner and outer zones; Ct1, Ct2: total compressibility of the inner zone and outer zone, MPa-1.

Equation 16 characterizes the flow capacity contrast between the two zones, and Equation 17 represents the dimensionless conductivity of the boundary. By applying dimensionless transformation to Equations 18. We obtain Equations 1824:

1rDrDrDP1DrD=P1DtD(18)
1rDrDrDP2DrD=ωP2DtD(19)
P1DrfD,tD=P2DrfD,tD(20)
P1DrDrD=rfD=P2DrDrD=rfD(21)
P2DrDrD=dD=CfDP2D(22)
CDdPWDdtD=1P1DrDrD=1(23)
PWDtD=P1DSP1DrD|rD=1(24)

2. Solution in Laplace Domain

By applying the Laplace transform to Equations 1822, the solutions for the inner and outer zones are obtained.

Inner zone solution:

P1D¯S=1SK0SrD+SS+CDSSSP1DSP1D¯rD|rD=1(25)

Outer zone solution:

P2D¯S=P1D¯S|rD=rfD=K0ωSrD/rfD+CfDωSK1ωSrD/rfD1+CfDωSK1ωSrD/rfDK0ωS(26)

2.3 Discussion on the solution of the non-sealing fault well test model

According to the law of images in seepage mechanics, the partial “image” reflection method is applied. Let α represent the degree of completeness of the image reflection. When α = 0, it reflects the scenario where the fault conductivity is entirely equal to the reservoir conductivity; when α = 1, it reflects the scenario where the fault is completely sealing; and when 0 < α < 1, it reflects the scenario of partial fault conductivity. After reflection, the fault disappears, and the pressure decline behavior of the original well is then solved using the superposition principle.

As shown in Figure 2, Well A is located in a two-zone composite reservoir, where the inner zone (Region 1) is near the wellbore, and the outer zone (Region 2) extends farther away. Well A is positioned at a distance d from Fault F. Well A′ is a “mirror image” reflection well generated with respect to the fault plane F. Assuming Well A produces at a rate q, to simulate fluid flow across the fault, the production rate of the “mirror image” well A′ is set as q' = αq, where α represents the degree of completeness of the mirror reflection.

Figure 2
Two circles labeled

Figure 2. Composite reservoir mirror image method schematic diagram.

Let η1=k1φμCt1, η2=k2φμCt2, M12=k/μ1k/μ2,the bottomhole pressure drop of the well A at time t after opening the well is:

ΔP=PiPwf=qμ1B345.6πk1hEirw214.4η1t+Eirf214.4η1tM12Eirf214.4η1texprf214.4η1η2tη2η1+aqμ2B345.6πk2hEi2d214.4η2texprf214.4η1η2tη2η1(27)

In Equation 27, the first term on the right-hand side represents the pressure drawdown at the wellbore induced by fluid flow from the near-well formation (Region 1), the second term corresponds to the pressure response from the distant formation (Region 2), while the third term characterizes the pressure interference caused by the image well (A′). During the early-time period before the pressure transient reaches the boundary of region 1 (t < t1), the contributions of both the second and third terms become negligible. By applying the logarithmic approximation to the exponential integral function -Ei(-x) and incorporating skin effects, Equation 27 can be simplified to Equation 28:

Pwf=Pi2.12×103qμ1Bk1hlgt+lgη1rw2+0.9077+0.8686S=Pim1lgt+lgη1rw2+0.9077+0.8686S(28)

Where, Pwf: bottom hole Pressure, MPa.

As shown in Figure 3, plotting the semi-logarithmic relationship between Pwf and lgt yields a straight-line segment with a slope of m1, where the expression for m1 is shown in Equation 29:

m1=2.12×103qμ1Bk1h(29)

Figure 3
Graph depicting bottom hole pressure (MPa) against log-transformed time. Segments labeled:

Figure 3. Schematic diagram of semi-logarithmic well test curve showing bottomhole pressure variation over time.

When the pressure drawdown from the producing well has propagated to the boundary of region 1 but before the pressure disturbance from the image well reaches the producing well (i.e., during the time period t < tF), the third term in Equation 27 becomes negligible. By approximating the exponential integral function -Ei(-x) with its logarithmic equivalent and incorporating skin effects, Equation 27 can be simplified to Equation 30:

Pwf=Pi2.12×103qμ2Bk2hlgt+2M12lgrfrw+lgη2rf2+0.9077+0.8686S=Pim2lgt+2M12lgrfrw+lgη2rf2+0.9077+0.8686S(30)

As shown in Figure 3, plotting the semi-logarithmic relationship between Pwf and lgt yields a straight-line segment with a slope of m2. The expression for m2 can be represented by Equation 31:

m2=2.12×103qμ2Bk2h(31)

As illustrated in Figure 3, m1, m2, and m3 correspond to the slopes of the straight-line segments in the semi-log plot, reflecting the pressure depletion rates during different flow regimes. These slope values quantitatively characterize the flow capacity of fluids passing through distinct reservoir regions: m1: near-wellbore zone; m2: outer reservoir region; m3: fault boundary zone. In this study, we focus on evaluating the flow conductivity near the fault based on this methodology. Thus, special emphasis is placed on the analysis of m2 and m3, which are critical for assessing. The bottom-hole flowing pressure is shown in Equation 32:

Pwf=Pi2.12×103qμ2Bk2hlgt+2M12lgrfrw+lgη2rf2+0.9077+0.8686S2.12×103αqμ2Bk2hlgt+2M12lgdrf+lgη2d2+0.3056=Pim3lgtm22M12lgrfrw+lgη2rf2+0.9077+0.8686Sαm22M12lgdrf+lgη2d2+0.3056(32)

As shown in Figure 3, plotting the semi-logarithmic relationship between Pwf and lgt yields a straight-line segment with a slope of m3. Where the expression for m3 can be represented by Equation 33:

m3=2.12×1031+αqμ2Bk2h=1+αm2(33)

Where the expression for a can be represented by Equation 34:

α=m3m21(34)

When α = 1, m3 = 2m2, indicating a completely sealed fault condition where no fluid flows across the fault; when α = 0, m3 = m2, reflecting that the fault’s flow conductivity is exactly equal to the reservoir’s flow conductivity, meaning fluid can completely pass through the fault; when 0 < α < 1, the relationship between m3 and m2 is non-integer, representing partial flow conductivity of the fault, which shows that the fault partially blocks fluid flow.

3 Analysis of pressure response characteristics in well testing models for non-sealing faults

The Stehfest numerical inversion algorithm was applied to Equations 25, 26 to obtain the time-domain solution of dimensionless bottomhole pressure, from which the dimensionless pressure derivative was calculated. Equation 35 presents the Stehfest Laplace numerical inversion algorithm. Subsequently, the characteristic response curves of the non-sealing fault well-testing model were generated, as illustrated in Figures 4, 5.

pwDt=ln2ti=1Mλip¯wDuiui=ln2/tiλi=1i+M/2k=i+1/2mini.M/2kM/22k!/M/2k!×k!k1!ik!2ki!(35)

Figure 4
Log-log plot showing different curves of PwD’ against tD/CD across a range of CfD values from 1.0 to 0.0, indicated by colors from orange to blue. The x-axis is labeled log(tD/CD) and the y-axis is labeled log(PWD) and log(PwD’). The plot highlights four sections marked by green dashed lines, labeled from one to four.

Figure 4. Characteristic well test response curves of composite reservoirs under various non-sealing fault conditions (outer zone mobility < inner zone mobility).

Figure 5
Log-log plot illustrating PwD’ values against log(tD/CD) with different CfD values. The legend shows CfD values ranging from 0.0 to 1.0. The graph contains colored lines representing each CfD value with sections labeled by numbers one to four, separated by green dashed vertical lines. Y-axis is labeled

Figure 5. Characteristic well test response curves of composite reservoirs under various non-sealing fault conditions (outer zone mobility > inner zone mobility).

Where, pwD: the dimensionless bottomhole pressure in Laplace space, MPa; z: the Laplace variable.

The well-testing curves (pressure and pressure derivative) for non-sealing faults can be divided into four characteristic stages:

1. Early-time wellbore storage and skin-dominated flow regime: The dimensionless pressure derivative curve initially exhibits a unit-slope straight line (indicating wellbore storage dominance); Subsequently transitions to an upward-convex curve (reflecting skin effects)

2. Mid-time radial flow regime: the dimensionless pressure derivative curve stabilizes into a horizontal line with zero slope (approaching 0.5); Demonstrates formation fluid flowing radially toward the wellbore

3. Transitional Response Phase: this stage reflects the flow transition from the inner zone to the outer zone in the composite reservoir system. When outer zone mobility is lower than inner zone mobility: the dimensionless pressure derivative curve initially exhibits an upward deflection, subsequently entering a stabilized pseudo-radial flow regime (Figure 4). When outer zone mobility exceeds inner zone mobility: the derivative curve first shows a “dip”, then stabilizes into pseudo-radial flow (Figure 5).

4. In the late outer boundary response stage, the characteristic curve of the dimensionless pressure derivative is affected by the mobility contrast between the outer and inner zone and the fault sealing parameter (CfD). When the mobility in the outer zone is smaller than that in the inner zone, as shown in Figure 4: when the outer boundary is a fully connected fault (CfD=1.0), the characteristic curve is a horizontal line with the value of 0.5k1μ2k2μ1; the curve of the fully sealing fault (CfD=0.0) is upwardly curved first, and finally stabilized at the horizontal line with the value of k1μ2k2μ1; the curve of the non-sealing fault (0.0<CfD<1.0) has a similar shape, which is upwardly curved first and then stabilized with the stabilized value between the two previous ones. On the contrary, as shown in Figure 5, when the flow in the outer zone is larger than that in the inner zone, the curve characteristic law is consistent with the above case. The characteristic curves of unconfined faults, which are ‘upturned and then stabilized, and the stabilization value is in the middle’, are the core markers that distinguish them from the confined faults.

Figures 6, 7 display the characteristic response curves of the non-sealing fault well-test model under varying boundary distances. Analytical results demonstrate that for both cases where outer zone mobility is either lower or higher than inner zone mobility, the onset time of dimensionless pressure derivative increase shows a significant positive correlation with boundary distance. Specifically, as the boundary distance progressively increases, the initiation point of the dimensionless pressure derivative rise exhibits a clearly delayed trend.

Figure 6
Log-log plot showing multiple curves of log(PwD) and log(PwD’) against log(tD/CD), with dD values from 1000.0 to 11000.0 represented by different colors. Green vertical dashed lines at significant points labeled 1 to 4.

Figure 6. Characteristic well test response curves of composite reservoirs under varying boundary distances (outer zone mobility < inner zone mobility).

Figure 7
Graph showing multiple lines representing different PwD’ values plotted against log(tD/CD) on a log-log scale. Each colored line corresponds to a specific dD value, ranging from 1000.0 to 11000.0. The y-axis is labeled log(PwD) and log(PWD’). Vertical green dashed lines divide the x-axis into sections labeled (1), (2), (3), and (4).

Figure 7. Characteristic well test response curves of composite reservoirs under varying boundary distances (outer zone mobility > inner zone mobility).

4 Intelligent evaluation methodology for fault sealing capacity

As previously discussed, for reservoirs incorporating an image well within a Cartesian coordinate system (Figure 8), the formation pressure P2(t) at any arbitrary point M (x, y) when t > tf can be mathematically represented by Equation 36:

P2t=Pi2.12×103qμ1Bk1hlg14.4η1tγx+d2+y2lg14.4η1tγrf2+M12lg14.4η1tγrf2exprf214.4η1η2tη2η2+0.8686Sα2.12×103qμ2Bk2hlg14.4η2tγxd2+y2exprf214.4η1η2tη2η2(36)

Figure 8
Diagram of a coordinate plane with point M marked at (x, y). Two circles centered at points A and A’ with radii r_A and r_A’ are shown. A is at (-d, 0) with blue lines; A’ is at (d, 0) with red lines. The diagram illustrates geometric relationships or forces.

Figure 8. Schematic diagram of composite reservoir Configuration and well locations in cartesian coordinates.

The derivative of Equation 35 in the x-direction can be expressed by Equation 37:

P2tx=1.842×103qμ1Bk1hx+dx+d2+y2+αM12xdxd2+y2exprf214.4η1η2tη2η1(37)

Thus, the pressure gradient at the location of the fault plane, i.e., on the y-axis, is expressed by Equation 38.

P2txx=0=1.842×103qμ1Bk1h×dd2+y2αM12dd2+y2exprf214.4η1η2tη2η1(38)

The flow rate across a unit length fault at any instant is expressed by Equation 39:

qF=k2hμ2P2txx=0=1.842×103qμ2BCfdk2hd2+y2×1αM12exprf214.4η1η2tη2η1(39)

Where, Cf is the conductivity of the reservoir, i.e., is the conduction coefficient or is the flow coefficient (k1h/μ1).

The conductivity of the fault, denoted as CF, satisfies Equation 40:

CF=1αM12exprf214.4η1η2tη2η1Cf(40)

Substituting Equation 33 into Equation 40 yields Equation 41:

CF=2m3m2Cf(41)

There β can be expressed by Equation 42.

β=m3m2(42)
CF=2βCf(43)

By determining the slope change magnification factor β from semi-log well test curves and substituting it into Equation 43, the fault conductivity or transmissibility across the fault can be obtained.

To quantitatively describe the observable connectivity between reservoir units achieved by fault conductivity, we define the effective fault thickness (h) as shown in Equation 44:

hb=CFμkF=2βCfμk=2βh(44)

In order to express the degree of connectivity of the fault to the reservoir on both sides, the fault connectivity coefficient is defined by the following Equation 45:

CfD=CFCf=2β(45)

To achieve rapid and accurate inversion of key fault parameters under complex non-sealing fault well-testing conditions, this paper proposes an intelligent hybrid modeling workflow. Based on the mathematical model established in Section 2, a theoretical type curve database was first generated. After feature extraction and normalization preprocessing, the XGBoost algorithm was employed for multi-output regression prediction. The PSO algorithm was then introduced to perform global search and optimal configuration of hyperparameters. Building upon the characteristic response features of multi-stage flow regimes, K-fold cross-validation was applied to thoroughly validate and optimize the model performance, thereby enhancing the robustness and generalization capability of fault parameter identification.

1. Dataset establishment

Based on the previously developed composite non-sealing fault well-test interpretation mathematical model, a numerical simulation method was employed to generate type curve data samples to support machine learning training and validation.

Sample Size: a total of 1,000 synthetic pressure and derivative curves for composite reservoirs were generated, covering various combinations of fault and reservoir properties, including: inner and outer zone permeability: 100–1,000 mD; Porosity: 0.15–0.30; Mobility ratio: 0.5–2; Fault boundary distance: 1,000–10,000 m.

2. Data processing and feature extraction

To enhance the efficiency and stability of model training, preprocessing was performed on the simulated curves and reservoir properties:

a. Feature denoising: wavelet transform was applied to the simulated pressure derivative curves to eliminate numerical fluctuations and improve feature extraction accuracy.

b. Input Feature Extraction: for each type curve, the following input features were extracted:

① reservoir property parameters: inner and outer zone permeability, porosity, mobility ratio and fault boundary distance.

② Flow regime characteristics: based on pressure derivative response, the curves were divided into: wellbore storage effect, radial flow regime, transitional flow regime, outer boundary effect regime, critical time points (m2, m3) for regime transitions and duration of each flow stage.

To eliminate the influence of differing dimensions and numerical scales across parameters, a normalization approach was implemented. Permeability values first underwent logarithmic transformation followed by standardization (Z-score normalization) to normalize their distribution. The mobility ratio and fault boundary distance were scaled to the [0,1] interval using min-max normalization to maintain consistent value ranges. For time-dependent features including characteristic transition points (m2, m3) and flow regime durations, global min-max normalization was applied based on the extreme values observed across the entire dataset.

3. Dataset partitioning and cross-validation

To rigorously evaluate the model’s generalization capability and prevent overfitting, this study employed 5-fold cross-validation (k = 5) for dataset partitioning and validation. The implementation process was as follows:

Random Partitioning: The type curve database was randomly divided into 5 equally sized subsets while maintaining data distribution consistency.

Iterative Training-Validation: in each iteration, 4 subsets were used as the training set, while the remaining 1 subset served as the validation set. This process was repeated until all 5 subsets had been used as the validation set exactly once. After completing all iterations, the average performance metrics across the 5 validation rounds were calculated to provide a comprehensive evaluation of the model’s robustness.

4. XGBoost model and hyperparameter configuration

This study employed XGBoost as a multi-output regression model to establish the nonlinear mapping relationship between type curve characteristics and fault-related parameters (m2, m3). The Mean Squared Error (MSE) was adopted as the loss function to quantify and minimize the discrepancy between the model’s predicted values and the actual target parameters (m2, m3).

5. PSO for hyperparameter tuning

To further enhance the generalization capability of the XGBoost model, this study employed the PSO algorithm for global search and automated optimization of key hyperparameters. The PSO implementation adopted the following configuration: A swarm of 30–50 particles was randomly initialized to explore the predefined parameter space. During each iteration, particle positions were dynamically updated based on both individual best and global best solutions, with the optimization process continuing until either fitness convergence or reaching the maximum iteration limit (50–100 cycles). The algorithm incorporated an inertia weight that linearly decayed from 0.9 to 0.4 and uniform learning factors (c1 = c2 = 1.5). The search space encompassed critical XGBoost parameters including: learning rate (0.01–0.2), maximum tree depth (4–8), and regularization parameters (L1: 0–1, L2: 0–2). This systematic approach balanced exploration and exploitation to identify optimal hyperparameter combinations that maximize model performance.

6. Optimal XGBoost model training and multi-target prediction

Upon obtaining the optimized hyperparameters, the XGBoost model was retrained using the complete training dataset while monitoring validation error to ensure optimal performance. The finalized model was then employed to conduct multi-target predictions on an independent test set, simultaneously outputting the critical parameters m2 and m3.

5 Application and analysis of fault seal evaluation in typical reservoirs

Faults are highly complex geological phenomena, and their conductivity is influenced by factors such as fracture stress, faulting mechanisms, fault orientation, types of minerals on the fault plane, and the distribution patterns of these minerals. No fault is absolutely sealing; each possesses a certain degree of conductivity, albeit varying in magnitude. While the location of a fault can be determined through seismic or other geological analysis methods and confirmed via well test analysis, the conductivity of a fault can only be determined through dynamic testing and analysis methods.

The X reservoir is located on the slope belt of the Grampian Arch, within a graben controlled by two NE-SW trending faults to the north and south. The oilfield features well-developed internal faults, including a set of NE-SW trending faults parallel to the graben and a contemporaneous set of NW-SE trending faults, forming a conjugate fault pattern. Planar distribution divides the oilfield into six fault blocks based on fault development characteristics from south to north. As shown in Figure 9: This displays a cross-section of the 3D fault-architecture model surrounding Well B1A, where Well B1A and Well B30Y are separated by Fault F1. The fault edges on both sides of F1 exhibit sandstone-to-sandstone juxtaposition, with a throw smaller than the reservoir thickness (28 m), indicating partial connectivity across the fault zone.

Figure 9
(a) A diagram shows a map with labeled zones connected by red lines. Black and white circles are scattered throughout, marked with letters and numbers. (b) A geological cross-section features colored layers representing different geological structures, such as sand-filled channels and mudstone. A legend on the side explains the color coding, including sand-filled channels in red, mudstone in grey, and more.

Figure 9. Spatial distribution map of well locations and fault characteristics. (a) Well locations and fault plane distribution. (b) Cross-sectional view of 3D fault-architecture model.

To rapidly quantify the transmissibility of the non-sealing fault F1, a pre-established intelligent hybrid inversion framework was implemented to analyze the composite reservoir between wells B1A and B30Y.

Permeability and mobility data from wells B1A and B30Y, along with the estimated fault throw, were first standardized or normalized to enhance the model’s stability and predictive accuracy. These processed parameters were then fed into a pretrained XGBoost multi-output regression model to obtain initial estimates of the key flow regime transition points, m2 and m3.

The initial predictions of m2 and m3 from the XGBoost model were subsequently used as the starting search centers for the Particle Swarm Optimization (PSO) algorithm. Taking the misfit between the measured pressure derivative curve and the simulated response as the fitness function, PSO iteratively updated particle positions to refine m2 and m3 estimates, ensuring convergence toward values more representative of the actual reservoir conditions.

To balance exploration efficiency and solution accuracy, the PSO algorithm employed a linearly decreasing inertia weight (from 0.9 to 0.4), with the maximum number of iterations set between 50 and 100 generations. During the training phase, five-fold cross-validation was incorporated to evaluate model generalization, and PSO was simultaneously utilized to perform automated hyperparameter tuning of the XGBoost model. This integrated modeling approach significantly enhanced the robustness and generalizability of the inversion results, demonstrating its effectiveness and practical value under complex fault-controlled reservoir conditions. As shown in Figure 10 are the inverted parameters m2 and m3.

Figure 10
Graph depicting pressure in megapascals (MPa) versus superposition time. The curve shows an upward trend, segmented into three slopes marked as \( m_1 \), \( m_2 = 0.225 \), and \( m_3 = 0.286 \), increasing toward \( t_F \).

Figure 10. Semi-Log Well Test Curve Schematic of Bottomhole Flowing Pressure vs. Time for Well B1A.

As evidenced by the analysis of Figure 10, the semi-log curve of Well B1A exhibits characteristics consistent with an infinite-acting composite reservoir system. The initial slope m1) corresponds to the formation properties of the inner zone, while the subsequent slope (m2) reflects the characteristics of the outer zone. The third slope (m3) clearly indicates the presence of a non-sealing fault boundary. Key reservoir parameters interpreted using our model include: inner zone permeability: 200 mD; composite radius: 220 m; outer zone permeability: 480 mD.

Figure 10 presents the pressure buildup test results from Well B1A conducted on 18 May 2013, where bottomhole flowing pressure was continuously monitored in real-time via a permanently installed downhole pressure gauge. Prior to shut-in, the well maintained an average production rate of 350 m3/d. The semi-log pressure derivative curve demonstrates characteristic behavior during the outer boundary response phase, exhibiting an initial upward deflection followed by stabilization - a diagnostic pattern indicative of non-sealing fault influence. This confirms the presence of a fault boundary south of Well B1A. The inversion framework quantified two distinct pressure depletion rates: a second-stage slope (m2) of 0.225 MPa representing the reservoir-dominated flow regime, and a third-stage slope (m3) of 0.286 MPa reflecting fault-affected flow. These inverted parameters (m2 = 0.225, m3 = 0.286) were subsequently input into Equations 4649, yielding the fault’s conductivity coefficient of 0.73 and effective connected thickness of 20.44 m.

β=m3m2=1.27(46)
CF=2βCf=0.73Cf(47)
hb=CFμkF=2βh=20.44(48)
CfD=CFCf=2β=0.73(49)

Previous geological architecture studies reveal that Fault F1 connects thicker, high-quality sandstone units on its northern side with thinner but equally permeable sandstones to the south, exhibiting partial hydraulic communication. Analysis indicates that Fault F1 has a geometric connectivity thickness of 0 m but an effective connected thickness of 24; Quantitative analysis demonstrates the fault’s connectivity coefficient (CfD) ranges between 0–1, with values trending closer to 1. This suggests that the fault is partially connected and has a high degree of connectivity, with a high degree of conductivity and a weak shading effect.

Based on this interpretation, development well B31Y was optimally positioned on the southern flank of Fault F1. After commissioning, the well achieved exceptional productivity of 2,000 cubic meters per day (m3/d), demonstrating highly effective reservoir development.

Figure 11 displays the bottom hole pressure profiles for Wells B1A and B30Y throughout 2014, with particular emphasis on the pressure buildup test initiated on 29 July 2014. Prior to shut-in, the wells operated at distinct production rates: B30Y maintained a robust output averaging 2,400 m3/d, while B1A produced at a more moderate rate of 480 m3/d. Continuous BHP monitoring was facilitated by permanently installed downhole pressure gauges in both wells, providing real-time, high-resolution pressure data throughout the test sequence.

Figure 11
Line graph showing bottom hole pressure (in megapascals) over time in days from January to December 2014. The blue line represents Wellbore B30Y and the red line represents Wellbore B1A. Pressure ranges from sixteen to twenty-seven megapascals, with notable fluctuations for B1A, particularly between April and September.

Figure 11. B1A and B30Y wellbore flow pressure variation curves.

It can be seen from the figure that the bottomhole flow pressure of the shut-in wells B1A and B30Y has a similar change trend but is not exactly the same, which indicates that there is pressure conduction between the two wells, but there is a certain difference in the change trend due to the shielding effect of fault F1. It further proves that fault F1 is a non-closed fault with strong pressure conduction ability, which verifies the reasonableness of the model calculation in this paper.

6 Conclusion

1. This study establishes a novel well-test interpretation model for composite reservoirs to characterize the limited conductivity of non-sealing faults. By incorporating the fluid flow mechanisms within fault zones into the analytical framework, the model significantly improves reservoir evaluation accuracy in the presence of conductive fault boundaries. The proposed approach provides new theoretical guidance and engineering methodologies for similar complex structural reservoirs.

2. A systematic analysis was conducted on the pressure and pressure derivative responses of composite reservoirs under varying fault conductivity and boundary distance conditions. The results demonstrate that non-sealing faults exhibit a diagnostic “upward deflection followed by stabilization” trend on the pressure derivative curve, serving as a key indicator for fault sealing evaluation. Additionally, the dimensionless derivative curve shows a strong positive correlation between the time of upward deflection and the fault boundary distance, providing a robust physical basis for quantitative inversion of fault location and sealing capacity.

3. An intelligent hybrid inversion framework (XGBoost-PSO) was developed and rigorously validated via a field case study in the X reservoir. Leveraging multi-stage pressure derivative characteristics and reservoir fluid properties, the inversion results quantitatively demonstrate that Fault F1 functions as a partially conductive boundary, with its conductivity coefficient and effective connected thickness confirming significant fault connectivity. Bottom-hole pressure history matching during the shut-in periods for Wells B1A and B30Y exhibits strong alignment with field monitoring data, validating the robustness and practical applicability of the proposed model and inversion framework. Future research should prioritize enhancing the model’s generalizability to ensure broader deployment across diverse reservoir scenarios.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

YY: Methodology, Writing – review and editing. BL: Writing – original draft. XW: Formal Analysis, Writing – review and editing. MQ: Formal Analysis, Writing – review and editing. DZ: Data curation, Formal Analysis, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded by the National Science and Technology Major Project of China, grant number 2017ZX05009004-005 and National Natural Science Foundation of China, grant numbers 52074322, U1762210.

Conflict of interest

Authors YY, XW, MQ, and DZ were employed by CNOOC International Limited.

The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Abdelaziz, B., and Tiab, D. (2004). “Pressure behavior of a well between two intersecting leaky faults,” in Paper presented at the Canadian International Petroleum Conference, Calgary, Alberta. doi:10.2118/2004-209

CrossRef Full Text | Google Scholar

Bixel, H. C., Larkin, B. K., and Van, H. K. (1963). Effect of linear discontinuities on pressure build-up and drawdown behavior. J. Pet. Technol. 15 (06), 885–895. doi:10.2118/611-PA

CrossRef Full Text | Google Scholar

Dake, L. P. (1978). Fumdamentals of reservoir engineering. Elsevier Scientific Publishing Company.

Google Scholar

Dong, P., and Liao, X. (2023). An efficient approach for automatic well-testing interpretation based on surrogate model and deep reinforcement learning. Chin. J. Comput. Phys. 40 (01), 67–80. doi:10.19596/j.cnki.1001-246x.8533

CrossRef Full Text | Google Scholar

Fu, X., Lv, D., and Huang, J. (2022). Fault-caprock coupling sealing mechanism and fault trap hydrocarbon accumulation model. Nat. Gas. Ind. 42 (03), 21–28. doi:10.3787/j.issn.1000-0976.2022.03.003

CrossRef Full Text | Google Scholar

Fu, X., Wang, H., and Song, X. (2024). Quantitative evaluation method and its application for vertical fault seal of thick caprock intervals. Petroleum Geol. and Oilfield Dev. Daqing 43 (04), 95–104. doi:10.19597/J.ISSN.1000-3754.202405013

CrossRef Full Text | Google Scholar

Horner, D. R. (1951). Pressure build-up in Wells 3rd world petroleum conference. Netherlands: Society of Petroleum Engineers, 503–521.

Google Scholar

Hu, X., Sandouhaxi, R., Liu, Y., and Lv, Y. (2024). Differences of lateral sealing ability between consequent fault and antithetic fault in Wen'an slope. Chin. J. Geol. 59 (4), 1057–1069. doi:10.12017/dzkx.2024.075

CrossRef Full Text | Google Scholar

Kuchuk, F. J., and Tarek, H. (1997). Pressure behavior of laterally composite reservoir. SPE Form. Eval. 12 (1), 47–56. doi:10.2118/24678-PA

CrossRef Full Text | Google Scholar

Li, X. (2023). Research on automatic well test analysis method based on deep learning. Anhui: Hefei University of Technology. doi:10.27101/d.cnki.ghfgu.2023.000182

CrossRef Full Text | Google Scholar

Li, D., Liu, X., and Cha, W. (2020). Automatic well test interpretation based on convolutional neural network for a radial composite reservoir. Petroleum Explor. Dev. 47 (03), 583–591. doi:10.11698/PED.2020.03.14

CrossRef Full Text | Google Scholar

Li, Y., Luo, L., and Wang, J. (2023). A geochemical evaluation method of fault sealing evolution and its controlling effect on hydrocarbon accumulation: a case study of the hongche fault zone in the northwest margin of the junggar basin. Nat. Gas. Ind. 43 (08), 12–25. doi:10.3787/j.issn.1000-0976.2023.08.002

CrossRef Full Text | Google Scholar

Liu, H., Wang, G., and Wu, S. (2016). A study on pressure transmission of unsteady fluid flow in banded fault block reservoir. J. Southwest Petroleum Univ. Sci. and Technol. Ed. 38 (05), 135–142. doi:10.11885/issn-1674-5086.2014.08.26.02

CrossRef Full Text | Google Scholar

Liu, Q., Jin, J., and Gui, F. (2017). Characteristics of well test curves corresponding to a combined reservoir with linear faults. J. Southwest Petroleum Univ. Sci. and Technol. Ed. 39 (04), 113–118. doi:10.11885/j.issn.1674-5086.2016.03.02.01

CrossRef Full Text | Google Scholar

Liu, Y., Zhao, S., and Zhou, W. (2024). Dynamic evaluation model of fault sealing in water drive oilfield. Fault-Block Oil and Gas Field 31 (05), 877–884. doi:10.6056/dkyqt202405017

CrossRef Full Text | Google Scholar

Ma, Q. (2025). Well test model of fractured well in linear shale oil reservoir with straight-line fault. Complex Hydrocarb. Reserv. 18 (01), 24–30. doi:10.16181/j.cnki.fzyqc.2025.01.004

CrossRef Full Text | Google Scholar

Mi, L., Gu, S., and Xu, L. (2021). An intelligent well test interpretation method based on data driven technology: a case study from one water-producing gas well of water-bearing gas reservoirs. Nat. Gas. Ind. 41 (02), 119–124. doi:10.3787/j.issn.1000-0976.2021.02.014

CrossRef Full Text | Google Scholar

Prasad, R. K. (1975). Pressure transient analysis in the presence of two intersecting boundaries. J. Pet. Technol. 27 (1), 89–96. doi:10.2118/4560-PA

CrossRef Full Text | Google Scholar

Stewart, G., Gupta, A., and Westaway, P. (1984). “The interpretation of interference tests in a reservoir with a sealing and a partially communicating faults,” in Presented at the 1984 european petroleum conference. London: Society of Petroleum Engineers, 411–429.

Google Scholar

Streltsova, T. D., and McKinley, R. M. (1984). Effect of flow time duration on buildup pattern for reservoirs with heterogeneous properties. SPEJ 24 (3), 294–306. doi:10.2118/11140-PA

CrossRef Full Text | Google Scholar

Wang, W., and Zhang, W. (1993). Application of semi logarithmic typical curve in well testing interpretation. Well Test. 2 (04), 42–46+50.

Google Scholar

Xing, C., Yin, H., and Zhang, Y. (2023). Well test analysis model for ribbon-shaped composite reservoirs with partially permeable faults. Petroleum Reserv. Eval. Dev. 13 (02), 215–222. doi:10.13809/j.cnki.cn32-1825/te.2023.02.010

CrossRef Full Text | Google Scholar

Yaxely, L. M. (1987). Effect of partially communicating fault on transient pressure behavior. SPE Form. Eval. 2 (4), 590–598. doi:10.2118/14311-PA

CrossRef Full Text | Google Scholar

Zeng, Y., Wang, B., and Nie, R. (2017). The porous flow model of multi-stage fractured horizontal well in linear composite oil reservoirs. Acta Pet. Sin. 38 (06), 687–695. doi:10.7623/syxb201706008

CrossRef Full Text | Google Scholar

Zhang, S. (2025). Application of pressure buildup test data in fine characterization of reservoirs in hade oilfield. Well Test. 34 (02), 73–78. doi:10.19680/j.cnki.1004-4388.2025.02.012

CrossRef Full Text | Google Scholar

Keywords: non-sealing fault, well test analysis, intelligent evaluation, Boltzmann transformation, “partial” image method

Citation: Yang Y, Li B, Wang X, Qi M and Zhai D (2025) The research on intelligent evaluation methods for non-sealing faults in reservoirs based on well test analysis. Front. Earth Sci. 13:1659387. doi: 10.3389/feart.2025.1659387

Received: 04 July 2025; Accepted: 19 August 2025;
Published: 09 September 2025.

Edited by:

Leilei Zhang, Rice University, United States

Reviewed by:

Zhengzheng Cao, Henan Polytechnic University, China
Xiaojun Wu, Changzhou University, China
Kun Qian, Changzhou University, China
Zongxiao Ren, Xi’an Shiyou University, China
Dongfeng Zhao, School of Food and Pharmacy, Zhejiang Ocean University, China

Copyright © 2025 Yang, Li, Wang, Qi and Zhai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Boying Li, MjAyMjMxMDE1OUBzdHVkZW50LmN1cC5lZHUuY24=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.