A method for quantitatively analyzing internal structure of a fault-controlled reservoir

Abstract

The application discloses a kind of inside structure quantitative analysis method of faulted reservoir, it includes: step one, respectively calculating fracture, fracture-cave and the corresponding sensitive attribute body of connectivity;Step two, the fracture, fracture-cave, connectivity sensitive attribute body after calculation is carried out fracture-cave connectivity probability calculation of multi-source potential field diffusion conduction, obtain the connectivity intensity body of corresponding three sensitive attribute data, quantitatively characterize fracture-cave body inside structure characteristics;Step three, the dominant band component reflecting geological target is gradually stripped out, and the inside multi-element gaussian function attribute body of fracture-cave body is fused and quantitatively calculated, to obtain the overall connectivity intensity of fracture-cave body inside structure;Step four, carry out the quantitative calculation of reservoir inside longitudinal and transverse fracture contact relationship, further quantitatively analyze reservoir inside structure characteristics.The stress attribute of fracture-cave body in the reservoir, energy body attribute and coherent attribute are organically combined together, the connectivity between fracture-cave is finely described and quantitatively analyzed.

Description

Technical Field

[0001] This invention relates to a method for quantitative analysis of the internal structure of fault-controlled reservoirs, belonging to the field of seismic exploration technology. Background Technology

[0002] Deep carbonate oil and gas reservoirs are a key area of ​​China's oil and gas development strategy. The carbonate reservoirs in the Shunbei area of ​​the Tarim Basin are generally buried at depths below 7500m, with highly developed fractures, high reserves, and great potential, making it an important oil and gas field in China. However, the internal structure of these reservoirs is extremely complex, mainly characterized by strong lateral heterogeneity, distinct fracture segmentation and fracture-vuggy separation, and the development of small-scale fractures and fracture-vuggy bodies within the reservoir's fractured zones, with unclear connectivity. Due to limitations in seismic data resolution and the complexity of the reservoir's internal structure, reservoir evaluation faces significant challenges. Furthermore, with the increasing demands for evaluation accuracy, fault-reservoir mechanisms and fault-reservoir relationships are becoming increasingly complex. The relationship between segmentation and high-yield enrichment lacks clear universal patterns, highlighting the growing importance of detailed characterization of fault-dissolution reservoirs.

[0003] Within the Shunbei main fault zone, the differences between faults such as No. 1, No. 4, No. 5, and No. 8 are significant, and even within the same fault, different parts exhibit marked variations. It is difficult to find commonalities and uniform patterns, resulting in a situation where each segment is a unique entity, and each fault is distinct. Therefore, detailed description and connectivity of the reservoir's internal structure are crucial for characterization. Only by clarifying the internal structure and connectivity of the reservoir can multiple targets be identified within a single well, maximizing single-well productivity. Previous studies on the comprehensive characterization of fault-controlled reservoirs have yielded numerous analyses, leading to the development of multi-attribute fusion identification and characterization techniques, such as "amplitude gradient attributes + likelihood + chaotic structure reflections," which have achieved good results. However, these attributes can only qualitatively describe the internal structure of fractures and cavities, generally limited to qualitative analysis, lacking sufficient accuracy and clarity regarding connectivity. Quantitative analysis methods are also lacking, failing to meet well design requirements. Therefore, a new method capable of quantitatively characterizing the internal structure of fault-controlled reservoirs is urgently needed to improve the accuracy of fracture and cavity structure description and guide subsequent well location and trajectory design. Summary of the Invention

[0004] The purpose of this invention is to provide a quantitative analysis method for the internal structure of fault-controlled reservoirs. For different types of reservoirs with pores and fractures, based on the reinforcement of fractures and pores, the method analyzes the connectivity strength between pores and fractures mainly by stress properties. Then, by combining seismic sensitivity properties, the method organically combines the stress properties, energy properties, and coherence properties of the pores and fractures inside the reservoir, so as to achieve a detailed description and quantitative analysis of the connectivity between pores and fractures.

[0005] The technical solution adopted in this invention is: a method for quantitative analysis of the internal structure of a fault-controlled reservoir, which includes the following steps:

[0006] Step 1: Perform coherence volume calculation on the interpretive results data volume, and use the calculated azimuth and dip information as constraints to calculate the sensitive attribute volumes corresponding to fractures, cracks, and connectivity respectively.

[0007] Step 2: Calculate the connectivity probability of fractures, cavities, and connectivity-sensitive attribute bodies using the multi-source potential field diffusion and propagation method, and obtain the connectivity strength of the corresponding three sensitive attribute data to quantitatively characterize the internal structural features of fracture-controlled cavities.

[0008] Step 3: For the connectivity strength data volume corresponding to fractures, fissures, and connectivity, principal component analysis is used to gradually extract the dominant frequency band components reflecting the geological targets, and quantitative calculation is performed by fusing the multivariate Gaussian function attribute volume inside the fracture-controlled fissure body to obtain the overall connectivity strength of the internal structure of the fissure body.

[0009] Step 4: Optimize the overall connectivity strength threshold value by combining reliable fracture elements in key areas, and on this basis, perform quantitative calculations of the longitudinal and transverse fracture contact relationships within the reservoir to further quantitatively analyze the internal structural characteristics of the reservoir.

[0010] Furthermore, step one includes:

[0011] (1) Calculate the range of line number, track number and time constraint on the explanatory results data, and select the inclination calculation range according to the actual construction of the work area to perform coherence calculation, and obtain the line inclination angle, track inclination angle, line inclination amount, track inclination amount and coherence value;

[0012] (2) Based on the calculated line inclination angle, track inclination angle, line inclination amount, track inclination amount and coherence value information, the sensitive attributes for fine description of fracture enhancement and crack enhancement are selected.

[0013] (3) Under the constraints of the calculated azimuth and dip angle information, a three-dimensional structural stress field calculation and analysis is performed to obtain the maximum principal stress and crack net pressure data volume;

[0014] (4) The optimized fracture, crack, and connectivity sensitive attribute bodies are calculated by three-dimensional structural stress field, and the planar features of the sensitive attribute bodies are extracted. The distribution features of fracture and crack bodies on the planar section are compared and verified for optimization analysis.

[0015] Furthermore, in step (2), the preferred criterion is to conduct hot attribute testing and analysis on reservoirs of different types and locations based on geological analysis and seismic facies identification models, and to select sensitive attributes that finely describe the internal and external contours, internal information and connectivity of fractures and solutions based on the attribute information and the intrinsic relationship between the research objectives.

[0016] Furthermore, in step (3), the numerical simulation method of tectonic stress adopts the theory of elastic thin plates to simulate the direction of multi-scale tectonic stress fields and calculate the magnitude of stress, adapting to various complex models. This method first assumes that the strata under study are uniform, continuous, isotropic, and perfectly elastic, and considers that the strata are formed entirely by tectonic stress. In the rectangular coordinate system, the mid-surface of the thin plate is taken as the coordinate plane z=0. According to the right-hand rule, the X and Y coordinates are parallel to the earth coordinates, with upward as positive, and the displacements along the positive X, Y, and Z directions are ux, uy, u... z .

[0017] Furthermore, the specific steps are as follows:

[0018] a. Build a construction model using the information from the construction interpretation results;

[0019] b. Using the seismic inversion results and basic data of rock samples from the target layer, determine the internal physical properties of the target area;

[0020] c. Combine the regional geological background of the target area to conduct tectonic evolution analysis on the tectonic model and obtain the displacement information of the target area model;

[0021] d. Based on the geological background information of the target work area, determine the tectonic movement period experienced by the target layer;

[0022] e. Using stratigraphic spatial information, combined with the internal physical and mechanical parameters of the target layer, finite element simulation and calculation are performed to obtain the tectonic stress field of the target layer.

[0023] Furthermore, to align with geomechanical notation, a notation convention of positive for compressive stress and negative for tensile stress is adopted; a curvature less than zero indicates an upward bulge in the strata. After determining the stress along the coordinate axis at that point, its principal stresses and their directions can be calculated.

[0024] (1)

[0025] Angle with the X-axis , Angle with the X-axis :

[0026] (2)

[0027] In the formula For stress, max For the maximum principal stress, min For the minimum principal stress, x The horizontal principal stress is located in the X-axis direction. y The horizontal principal stress is in the y-axis direction. Shear force in the xy direction;

[0028] Therefore, by obtaining the perturbation equation of the ground plane or the curvature of points on its surface, the three-dimensional tectonic stress field on it can be estimated, and then the cracks caused by the stress can be calculated.

[0029] Furthermore, step two includes:

[0030] (1) The thermal diffusion coefficient of fractures, cracks, and connectivity-sensitive properties after the three-dimensional structural stress field calculation is calculated by introducing the thermal diffusion and conduction formula. The connectivity coefficient of cracks is obtained. If the thermal conductivity k is a constant, the thermal diffusion equation can be expressed as follows:

[0031] (3)

[0032] In the formula, x, y, and z are coordinate axes; T is the diffusion field; t is time; α is called the thermal diffusivity, which represents the ability of an object to achieve uniform temperature during non-steady-state heat conduction. That is, the larger the thermal conductivity, the faster the temperature tends to become uniform; q is the rate of heat generation per unit volume; the thermal diffusion equation shows that at any point in the medium, the net rate of heat conduction into a unit volume plus the rate of heat generation per unit volume must be equal to the rate of change of energy stored in the unit volume.

[0033] (2) Normalize the connection coefficient of the gap. Specifically, divide the scale value of the connection coefficient of the gap by the maximum value to obtain the connection coefficient value of the gap with a maximum value of 1.

[0034] Furthermore, step four includes:

[0035] (1) Using the clearly defined fracture elements in the key areas of the study area, the calculated overall connectivity strength is calibrated and verified by longitudinal profile and transverse plane connectivity strength, and the connectivity strength threshold value is determined;

[0036] (2) On the connected strength body after the threshold value is determined, the interval distance between fractures can be directly displayed on the flat section, including the width of the transverse fracture zone between fractures, the size of the interval between fracture contact points, and the size of the contact and misalignment distance of longitudinal fractures, which can intuitively and effectively guide the design of well location and well trajectory in the later stage.

[0037] The beneficial effects of this invention are as follows: This method addresses the current situation where the accuracy of describing the internal structure of fault-controlled reservoirs is insufficient, failing to meet the needs of later exploration and development. Starting with the net pressure and thermal diffusivity within the fractures, it quantitatively calculates the fracture-cavity connectivity probability, internal connectivity strength of the fracture-cavity body, and fracture contact relationships through multi-source potential field diffusion and transmission. This allows for the rapid and effective separate and detailed characterization and description of different types of pore-fracture reservoirs. Furthermore, the operating software is simple and highly scalable, making it of great significance. The improvement and promotion of this method can intuitively and effectively quantitatively characterize the internal structure of fault-controlled fracture-cavity bodies, effectively guiding later development and well location design, which is of great importance and has broad prospects.

[0038] Additional aspects and advantages of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. Attached Figure Description

[0039] Figure 1 The flowchart shows the quantitative analysis method for the internal structure of the control and storage group.

[0040] Figure 2 The diagrams are cross-sectional views of various properties of cracks and fractures; among them, (A) is the structural tensor cross-section (crack outline), (B) is the crack reinforcement cross-section (crack longitudinal characterization), (C) is the COR_thin cross-section, and (D) is the discontinuity cross-section.

[0041] Figure 3 The diagrams show the maximum principal stress and net pressure within the joint; where (A) is the maximum principal stress and (B) is the net pressure within the joint.

[0042] Figure 4 The diagram shows the sensitive properties of fracture, cavity, and net pressure within the cavity; where (A) represents Likelihood, (B) represents cavity reinforcement, and (C) represents net pressure within the cavity.

[0043] Figure 5 The diagram shows the seismic profile, fault, cavity, and net pressure connectivity within the fracture; where (A) is the seismic profile, (B) is the fracture-cavity enhanced connectivity, (C) is the net pressure connectivity within the fracture, and (D) is the fault-enhanced connectivity.

[0044] Figure 6 This is a three-dimensional display of the overall connectivity strength.

[0045] Figure 7 This is a schematic diagram for quantifying contact relationships.

[0046] Figure 8 This is a plan view of the overall connectivity strength of the HS region.

[0047] Figure 9A three-dimensional sculpted map of the overall connectivity strength of the HS region. Detailed Implementation

[0048] The present invention will now be described in detail with reference to the accompanying drawings and specific embodiments.

[0049] A quantitative analysis method for the internal structure of a fault-controlled reservoir, the method comprising:

[0050] Step 1: Perform coherence volume calculation on the interpretive results data volume, and use the calculated azimuth and dip information as constraints to calculate the sensitive attribute volumes corresponding to fractures, cracks, and connectivity, respectively.

[0051] (1) Calculate the range of line number, track number and time constraint on the explanatory results data, and select the inclination calculation range according to the actual construction of the work area to perform coherent body calculation, and obtain the line inclination angle, track inclination angle, line inclination amount, track inclination amount and coherent value.

[0052] (2) Combining the calculated line dip angle, trace dip angle, line inclination, trace inclination, and coherence value information, the sensitive attributes for fine description of fracture enhancement and fracture-cavity enhancement are selected. The selection criteria are as follows: based on geological analysis and seismic facies identification model, hot attribute tests and analyses are conducted on reservoirs of different types and locations of cavities and fractures. According to the attribute information and the intrinsic relationship of the research objectives, the sensitive attributes for fine description of the internal and external contours, internal information, and connectivity of fracture-cavities are selected. In terms of contour sensitive attribute selection, seismic attributes can clearly characterize the external boundary features of fracture-controlled fracture-cavities, obtain the relatively accurate location of the original fracture surface and the maximum abrupt change surface of fracture-cavities, and effectively identify the width and morphological characteristics of fracture-cavities. In terms of sensitive attribute selection of internal information of fracture-controlled fracture-cavities, seismic attributes can clearly characterize the degree of internal fragmentation of fracture-cavities, and the information is rich and reliable. In terms of sensitive attribute selection of internal connectivity of fracture-controlled fracture-cavities, seismic attributes can effectively divide the stretching and uplift areas of strike-slip fault zones, and can more accurately clarify the connectivity and connectivity size between reservoir fractures and fracture-cavities. Figure 2 ).

[0053] (3) Under the constraints of the calculated azimuth and dip angle information, a three-dimensional structural stress field calculation and analysis is performed to obtain the maximum principal stress and net crack pressure data volume ( Figure 3 ).

[0054] The numerical simulation method for tectonic stress, based on the theory of elastic thin plates, simulates the direction and magnitude of multi-scale tectonic stress fields (paleogeneous stress), adapting to various complex models. This method first assumes that the studied strata are homogeneous, continuous, isotropic, and perfectly elastic, and that the formation of the strata is entirely due to tectonic stress. In a rectangular coordinate system, the mid-surface of the thin plate is taken as the coordinate plane z=0. Following the right-hand rule, the X and Y coordinates are defined as parallel to the geodetic coordinates, with upward as positive. The displacements along the positive X, Y, and Z directions are ux, uy, and u, respectively. z The specific steps are as follows:

[0055] a. Build a construction model using the information from the construction interpretation results;

[0056] b. Using seismic inversion data and basic data such as rock samples from the target layer, determine the internal physical properties of the target area;

[0057] c. Combine the regional geological background of the target area to conduct tectonic evolution analysis on the tectonic model and obtain the displacement information of the target area model;

[0058] Based on the geological background information of the target work area, determine the tectonic movement period of the target layer;

[0059] e utilizes stratigraphic spatial information, combined with the internal physical and mechanical parameters of the target layer, to perform finite element simulation and calculation, thereby obtaining the tectonic stress field of the target layer.

[0060] To align with geomechanical notation, a notation convention of positive compressive stress and negative tensile stress is adopted; a curvature less than zero indicates an upward bulge in the strata. After determining the stress along the coordinate axis at that point, its principal stresses and their directions can be calculated.

[0061] (1)

[0062] Angle with the X-axis , Angle with the X-axis :

[0063] (2)

[0064] In the formula For stress, max For the maximum principal stress, min For the minimum principal stress, x The horizontal principal stress is located in the X-axis direction. y The horizontal principal stress is in the y-axis direction. This represents the shear force in the xy direction.

[0065] Therefore, if the perturbation equation of the stratum or the curvature of its points can be obtained, the three-dimensional tectonic stress field on it can be estimated, and then the fractures generated by the stress can be calculated. The quantitative calculation of tectonic stress has a good effect on fracture propagation law, reservoir internal connectivity and reservoir fine prediction.

[0066] (4) The optimized fracture, crack, and connectivity sensitive attribute bodies were calculated by three-dimensional structural stress field calculation, and the planar features of the sensitive attribute bodies were extracted. The analysis was compared and verified by combining the distribution features of fracture and crack bodies on the planar section. Figure 4 ).

[0067] Step 2: Using the multi-source potential field diffusion and propagation method, calculate the connection probability of fractures, fissures, and connectivity-sensitive attribute volumes after the three-dimensional structural stress field calculation. Obtain the corresponding connectivity strength volumes of the three sensitive attribute data, and quantitatively characterize the internal structural features of the fracture-controlled fissure bodies. Figure 5 ).

[0068] (1) The thermal diffusion coefficient of fractures, cracks, and connectivity-sensitive properties after the three-dimensional structural stress field calculation is calculated by introducing the thermal diffusion and conduction formula. The connectivity coefficient of cracks is obtained. If the thermal conductivity k is a constant, the thermal diffusion equation can be expressed as follows:

[0069] (3)

[0070] In the formula, x, y, and z are coordinate axes; T is the diffusion field; t is time; α is called the thermal diffusivity, which represents the ability of an object's internal temperature to become uniform during non-steady-state heat conduction. That is, the larger the thermal conductivity, the faster the temperature tends to become uniform; q is the rate of heat generation per unit volume; the thermal diffusion equation shows that at any point in the medium, the net rate of heat conduction into a unit volume plus the rate of heat generation per unit volume must be equal to the rate of change of energy stored in the unit volume.

[0071] (2) Normalize the connection coefficient of the gap. Specifically, divide the scale value of the connection coefficient of the gap by the maximum value to obtain the connection coefficient value of the gap with a maximum value of 1.

[0072] Step 3: Since pressure is interrelated with fractures and fissures, it has a multivariate normal effect on connectivity. For the connectivity strength data volume corresponding to fractures, fissures, and connectivity, principal component analysis is used to gradually extract the dominant frequency band components reflecting geological targets. A quantitative calculation of the multivariate Gaussian function attribute volume within the fracture-controlled fissure / cavity body is then performed to obtain the overall connectivity strength of the internal structure of the fissure / cavity body. Figure 6 ).

[0073] Step 4: Optimize the overall connectivity strength threshold value by combining reliable fracture elements in key areas, and on this basis, perform quantitative calculations of the longitudinal and transverse fracture contact relationships within the reservoir, further quantitatively analyzing the internal structural characteristics of the reservoir. Figure 7 ).

[0074] (1) Using the clearly defined fracture elements in the key areas of the study area, the calculated overall connectivity strength is calibrated and verified by longitudinal profile and transverse plane connectivity strength, and the connectivity strength threshold value is determined.

[0075] (2) On the connected strength body after the threshold value is determined, the interval distance between fractures can be directly displayed on the flat section, including the width of the transverse fracture zone between fractures, the size of the fracture contact point interval, and the size of the contact and misalignment distance of the longitudinal fracture, which can intuitively and effectively guide the design of well location and well trajectory in the later stage.

[0076] To verify the application effect of the new method for quantitative analysis of the internal structure of fault-controlled reservoirs proposed in this application, reservoir sensitivity attribute testing and sensitivity attribute optimization analysis were carried out in the HS study area. The sensitivity attributes of fracture, fracture-vuggy, and connectivity in this area were determined to be structural tensor, likelihood, and net pressure of fracture-vuggy, respectively. The connectivity strength was calculated using the three sensitivity attribute volumes, and the overall connectivity strength was finally obtained and displayed by 3D engraving. Figure 8 This is a planar distribution map of the overall connectivity strength in the area. It clearly shows the connectivity strength of the main fault zone and different parts within it, effectively displaying the extent of connectivity between well groups, which is a significant improvement over previous qualitative descriptions. Figure 9 This is a three-dimensional display of the overall connectivity intensity of the area. From the three-dimensional connectivity intensity diagram of the area, we can see the connectivity size of different directions and different parts of the fracture-controlled cavity. It provides intuitive and effective guidance for the design of well locations and well trajectories in the area, with remarkable results.

[0077] In summary, the method of this invention targets different types of reservoirs with pores and fractures. Based on fracture and fracture-cavity reinforcement, it analyzes the connectivity strength between fractures and cavities primarily using stress properties. Furthermore, by combining seismic sensitivity properties, it organically integrates the stress properties, energy properties, and coherence properties of the fracture-cavity body within the reservoir, thereby achieving a detailed description and quantitative analysis of the connectivity between fractures and cavities.

[0078] This invention further clarifies the internal structure and connectivity of fault-controlled fractured caverns in deep carbonate rocks of the Tarim Basin, guiding well location and trajectory design, significantly improving single-well productivity and recovery rate, and effectively enhancing economic benefits. The refinement and widespread application of this method enable intuitive and effective quantitative characterization of the internal structure of fault-controlled fractured caverns, effectively guiding subsequent development and well location design. It fills a gap in the quantitative calculation of the internal structure and connectivity of fault-controlled reservoirs both domestically and internationally, possessing significant importance and promising prospects.

[0079] The foregoing has shown and described the basic principles, main features, and advantages of the present invention. Those skilled in the art should understand that the above embodiments do not limit the scope of protection of the present invention in any way, and all technical solutions obtained by equivalent substitution or other means fall within the scope of protection of the present invention. Parts not covered in this invention are the same as or can be implemented using existing technology.

Claims

1. A method for quantitative analysis of the internal structure of a fault-controlled storage group, characterized in that... Includes the following steps: Step 1: Perform coherence volume calculations on the interpretive results data volume, and use the calculated azimuth and dip information as constraints to calculate the sensitive attribute volumes corresponding to fractures, fissures, and connectivity, respectively; specifically including: (1) Calculate the range of line number, track number and time constraint on the explanatory results data, and select the inclination calculation range according to the actual construction of the work area to perform coherence calculation, and obtain the line inclination angle, track inclination angle, line inclination amount, track inclination amount and coherence value; (2) Based on the calculated line inclination angle, track inclination angle, line inclination amount, track inclination amount and coherence value information, the sensitive attributes for fine description of fracture enhancement and crack enhancement are selected. (3) Under the constraints of the calculated azimuth and dip angle information, a three-dimensional structural stress field calculation and analysis is performed to obtain the maximum principal stress and crack net pressure data volume; (4) The optimized fracture, crack, and connectivity sensitive attribute bodies are calculated by three-dimensional structural stress field, and the planar features of the sensitive attribute bodies are extracted. The distribution features of fracture and crack bodies on the planar section are compared and verified for optimization analysis. Step 2: Utilizing the multi-source potential field diffusion and propagation method, calculate the connectivity probability of fractures, voids, and connectivity-sensitive attribute volumes using multi-source potential field diffusion and propagation. Obtain the connectivity strength volumes of the corresponding three sensitive attribute data, quantifying the internal structural characteristics of the fracture-controlled void body. Specifically, this includes: (1) The thermal diffusion coefficient of fractures, cracks, and interconnectivity-sensitive properties after the three-dimensional structural stress field calculation is calculated by introducing the thermal diffusion and conduction formula. The interconnectivity coefficient of cracks is obtained. If the thermal conductivity k is a constant, the thermal diffusion equation is expressed as follows: （3） In the formula, x, y, and z are coordinate axes; T is the diffusion field; t is time; α is called the thermal diffusivity, which represents the ability of an object to achieve uniform temperature during non-steady-state heat conduction. That is, the larger the thermal conductivity, the faster the temperature tends to become uniform; q is the rate of heat generation per unit volume; the thermal diffusion equation shows that at any point in the medium, the net rate of heat conduction into a unit volume plus the rate of heat generation per unit volume must be equal to the rate of change of energy stored in the unit volume. (2) Normalize the connection coefficient of the gap. Specifically, divide the scale value of the connection coefficient of the gap by the maximum value to obtain the connection coefficient value of the gap with a maximum value of 1. Step 3: For the connectivity strength volume corresponding to fractures, fissures, and connectivity, principal component analysis is used to gradually extract the dominant frequency band components reflecting the geological target, and quantitative calculation of the multivariate Gaussian function attribute volume inside the fracture-controlled fissure body is performed to obtain the overall connectivity strength of the internal structure of the fissure body. Step 4: Optimize the overall connectivity strength threshold value by combining reliable fracture elements in key areas, and on this basis, perform quantitative calculations of the longitudinal and transverse fracture contact relationships within the reservoir to further quantitatively analyze the internal structural characteristics of the reservoir.

2. The method for quantitative analysis of the internal structure of a fault-controlled reservoir according to claim 1, characterized in that, In step (2) of step one, the preferred criterion is to conduct hot attribute testing and analysis on reservoirs of different types and locations based on geological analysis and seismic facies identification models, and to select sensitive attributes that finely describe the internal and external contours, internal information and connectivity of fractures and solutions based on attribute information and the intrinsic relationship between research objectives.

3. The method for quantitative analysis of the internal structure of a fault-controlled reservoir according to claim 1, characterized in that, In step (3) of step one, the three-dimensional tectonic stress field calculation adopts the elastic thin plate theory to simulate the direction of multi-scale tectonic stress field and calculate the stress magnitude, adapting to various complex models. This method first assumes that the studied strata are uniform, continuous, isotropic, and perfectly elastic, and considers that the strata are formed entirely by tectonic stress. In the rectangular coordinate system, the mid-surface of the thin plate is set as the coordinate plane z=0. According to the right-hand rule, the X and Y coordinates are parallel to the earth coordinates, with upward as positive, and the displacements along the positive X, Y, and Z directions are u, respectively. x u y u z .

4. The method for quantitative analysis of the internal structure of a fault-controlled reservoir according to claim 3, characterized in that, The specific steps are as follows: a. Build a construction model using the information from the construction interpretation results; b. Using the seismic inversion results and basic data of rock samples from the target layer, determine the internal physical properties of the target area; c. Combine the regional geological background of the target area to conduct tectonic evolution analysis on the tectonic model and obtain the displacement information of the target area model; d. Based on the geological background information of the target work area, determine the tectonic movement period experienced by the target layer; e. Using stratigraphic spatial information, combined with the internal physical and mechanical parameters of the target layer, finite element simulation and calculation are performed to obtain the tectonic stress field of the target layer.

5. The method for quantitative analysis of the internal structure of a fault-controlled reservoir according to claim 4, characterized in that, To align with geomechanical notation, a notation convention of positive compressive stress and negative tensile stress is adopted; a curvature less than zero indicates an upward bulge in the strata. After determining the stress along the coordinate axis at that point, its principal stresses and their directions can be calculated. （1） Angle with the X-axis , Angle with the X-axis : （2） In the formula For stress, max For the maximum principal stress, min For the minimum principal stress, x The horizontal principal stress is located in the X-axis direction. y The horizontal principal stress is in the y-axis direction. Shear force in the xy direction; Therefore, by obtaining the perturbation equation of the ground plane or the curvature of points on its surface, the three-dimensional tectonic stress field on it can be estimated, and then the cracks caused by the stress can be calculated.

6. The method for quantitative analysis of the internal structure of a fault-controlled reservoir according to claim 1, characterized in that, Step four includes: (1) Using the clearly defined fracture elements in the key areas of the study area, the calculated overall connectivity strength is calibrated and verified by longitudinal profile and transverse plane connectivity strength, and the connectivity strength threshold value is determined; (2) On the connected strength body after the threshold value is determined, the interval distance between fractures can be directly displayed on the flat section, including the width of the transverse fracture zone between fractures, the size of the interval between fracture contact points, and the size of the contact and misalignment distance of longitudinal fractures, which can intuitively and effectively guide the design of well location and well trajectory in the later stage.

Images