WO2021008630A1 - 圈闭断层的封闭性定量分析方法、装置及*** - Google Patents
圈闭断层的封闭性定量分析方法、装置及*** Download PDFInfo
- Publication number
- WO2021008630A1 WO2021008630A1 PCT/CN2020/111223 CN2020111223W WO2021008630A1 WO 2021008630 A1 WO2021008630 A1 WO 2021008630A1 CN 2020111223 W CN2020111223 W CN 2020111223W WO 2021008630 A1 WO2021008630 A1 WO 2021008630A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- fault
- mudstone
- smear
- core
- data volume
- Prior art date
Links
- 238000007789 sealing Methods 0.000 title claims abstract description 135
- 238000000034 method Methods 0.000 title claims abstract description 113
- 238000004445 quantitative analysis Methods 0.000 title claims abstract description 43
- 239000011148 porous material Substances 0.000 claims abstract description 200
- 238000005070 sampling Methods 0.000 claims abstract description 150
- 238000004088 simulation Methods 0.000 claims description 90
- 238000004458 analytical method Methods 0.000 claims description 58
- 239000007787 solid Substances 0.000 claims description 38
- 238000003384 imaging method Methods 0.000 claims description 19
- 238000005553 drilling Methods 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 13
- 239000000463 material Substances 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000004891 communication Methods 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims 1
- 239000011435 rock Substances 0.000 abstract description 16
- 230000008569 process Effects 0.000 description 40
- 238000010586 diagram Methods 0.000 description 32
- 238000012545 processing Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 8
- 238000005755 formation reaction Methods 0.000 description 7
- 238000006073 displacement reaction Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 239000012530 fluid Substances 0.000 description 6
- 238000005325 percolation Methods 0.000 description 6
- 238000004451 qualitative analysis Methods 0.000 description 6
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 239000003208 petroleum Substances 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 3
- 230000035699 permeability Effects 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- 239000006004 Quartz sand Substances 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000004927 clay Substances 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000011158 quantitative evaluation Methods 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 230000032258 transport Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 244000035744 Hura crepitans Species 0.000 description 1
- 208000002599 Smear Layer Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003032 molecular docking Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Definitions
- This application relates to the technical field of oil and gas exploration and development, and in particular to a method, device and system for quantitative analysis of the sealing of trap faults.
- fault sealing is an important factor in the formation of oil and gas reservoirs and controlling the scale of oil and gas reservoirs. Its research is of great significance for understanding faults in oil and gas accumulation and petroleum exploration and development. Therefore, the study of fault sealing is very important. At the same time, there are many factors that affect the sealability of faults, and the analysis of these factors is also very necessary. What role faults play in the process of oil and gas migration and accumulation depends on their sealing properties.
- fault sealability refers to the ability of fault planes or fault zones to seal formation fluids and prevent fluid seepage.
- trap fault is the fault of stratum. After the stratum is fractured, a fault is formed by moving it again. There is only a fracture, no movement, that is a crack. The extension range of fractures is small, and the extension range of trap faults is large. Fractures usually act as seepage channels, while trap faults have more complicated functions, sometimes sealing oil and gas, and sometimes transporting oil and gas. Faults are dislocated, which is the result of shearing. Shearing can cause rocks on both sides to be broken. Some clastics fall off and fill the cracks, and then evolve into wider or narrow rock layers. Therefore, the trap fault is a clastic rock layer formed after being filled. Clastic rocks have porosity and permeability.
- the physical properties of the fault clastic rocks are also very different. They are heterogeneous layers. Some parts are low porosity and low permeability, and some parts are high porosity and high permeability.
- the fault rock facing the reservoir is clastic rock with relatively poor physical properties
- the fault plays a role of sealing oil and gas.
- the fault rock facing the reservoir is clastic rock with relatively good physical properties
- the fault plays a role in transporting oil and gas.
- faults cannot seal water, but can only seal oil and gas.
- the mechanism for faults to seal oil and gas is exactly the same as that of caprocks.
- a fault transports or seals oil and gas has nothing to do with the occurrence of the fault, has nothing to do with the positive and negative properties, and has nothing to do with the tension and compression.
- Underground faults are all closed and cannot be opened. Closed faults can seal oil and gas, and can also transport oil and gas.
- Tectonic movement may destroy the fault again, but it only changes the properties of the fault rock. For example, it may produce cracks in the fault rock, but it is impossible to change the closure properties of the fault.
- the sealability of trap faults has always been one of the core issues of petroleum geology research, and has attracted much attention from domestic and foreign petroleum geologists. 80% of fault block oil and gas reservoirs are controlled by fault sealability, and the trap fault rock mudstone smear type is Mainly closed type.
- the analysis of the sealability of the trap fault is directly related to the ability to accurately judge whether the trap controlled by the fault can accumulate and the hydrocarbon height. In the practice of petroleum industry exploration and production, for fault block traps, whether it is possible to accurately predict the oil and gas height and the size of the trap area is critical to improving the efficiency of oil and gas exploration, and will directly affect exploration decision deployment and investment .
- the embodiment of the present invention proposes a method for quantitatively analyzing the sealing properties of trap faults, which is used to quantitatively analyze the sealing properties of trap faults to be drilled with high accuracy.
- the method includes:
- each sampling point analyzes the sealability of each sampling point of the mudstone smear type fault to be drilled.
- the embodiment of the present invention provides a quantitative analysis device for the sealability of trap faults, which is used to quantitatively analyze the sealability of the trap fault to be drilled with high accuracy.
- the device includes:
- the first module is used to obtain a three-dimensional data volume of the mudstone smear-type fault in the trap to be drilled, and the three-dimensional data volume is obtained by scanning and reconstructing the fault simulation model of the mudstone smear-type fault in the trap to be drilled;
- the second module is used to obtain the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points according to the three-dimensional data volume;
- the third module is used to fit the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points, and determine the fitting coefficient of the mudstone smear thickness and fault gouge ratio of the fault zone;
- the fourth module is used to determine the mudstone smear sealing factor of the fracture zone of multiple sampling points according to the fitting coefficient and the fault gouge ratio of multiple sampling points of the mudstone smear type fault to be drilled;
- the fifth module is used to analyze the sealability of each sampling point of the mudstone smear-type fault to be drilled according to the mudstone smear sealing factor of each sampling point.
- the embodiment of the present invention provides a quantitative analysis system for the sealability of trap faults, which is used to quantitatively analyze the sealability of the trap fault to be drilled with high accuracy.
- the system includes: a fault simulation model structural unit, a scanning unit and the above-mentioned traps
- a quantitative analysis device for the sealing of faults in which,
- the fault simulation model structural unit is used to construct the fault simulation model of the mudstone smear type fault to be drilled;
- the scanning unit is used to scan the fault simulation model to generate the surface data of the mudstone smear type fault to be drilled and send it to the sealing quantitative analysis device of the trap fault;
- a quantitative analysis device for the sealing of trap faults is used to reconstruct the received surface data to obtain a three-dimensional data volume; according to the three-dimensional data volume, obtain the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points;
- the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points are fitted to determine the fitting coefficient of the mudstone smear thickness and the fault gouge ratio of the fault zone; according to the fitting coefficient and the mudstone smear type fault to be drilled
- the ratio of fault mud at multiple sampling points determines the mudstone smear sealing factor of the fracture zone at multiple sampling points; according to the mudstone smear sealing factor of the fracture zone at each sampling point, analyze each sampling point of the mudstone smear type fault to be drilled The closure.
- the embodiment of the present invention proposes a fault sealing analysis method, which is used to quantitatively analyze the fault sealing performance with high accuracy.
- the method includes:
- the fault sealing is analyzed.
- the embodiment of the present invention provides a fault sealing analysis device, which is used to quantitatively analyze the fault sealing performance with high accuracy.
- the device includes:
- the first data obtaining module is used to obtain a first core data volume containing pores, and the first core data volume containing pores is obtained by scanning a core model of a target interval containing faults;
- the second data acquisition module is used to perform three-dimensional reconstruction of the first core data volume containing pores to obtain a solid second core data volume;
- the third data acquisition module is used to make a difference between the second solid core data volume and the first core data volume containing pores to obtain the pore data volume;
- the core fracture zone seepage field acquisition module is used to obtain the core fracture zone seepage field according to the pore data volume
- Calculation module used to calculate the strength of the core fracture zone seepage field at the fracture zone
- the analysis module is used to analyze the fault sealability according to the strength of the core fracture zone's seepage field at the fracture zone.
- the embodiment of the present invention provides a fault sealing analysis system for quantitative analysis of fault sealing with high accuracy.
- the system includes: the above-mentioned fault sealing analysis device and a scanning unit, wherein:
- the scanning unit is used for:
- the embodiment of the present invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and running on the processor, and the processor realizes the closure of the trap fault when the processor executes the computer program. Quantitative analysis method or fault sealing analysis method.
- the embodiment of the present invention also provides a computer-readable storage medium that stores a computer program that executes the above-mentioned trapped fault sealing quantitative analysis method or fault sealing analysis method.
- a first data volume of the mudstone smear-type fault in the trap to be drilled is obtained, and the first data volume is obtained by scanning and reconstruction of the fault simulation model of the mudstone smear-type fault in the trap to be drilled.
- the fault simulation model contains pores, and the first data volume is a data volume containing pores; the first data volume is three-dimensionally reconstructed to obtain a second data volume, and the second data volume is a solid data volume; The second data volume is compared with the first data volume to obtain the pore data volume; according to the pore data volume, the core fracture zone seepage field is obtained, and the strength of the core fracture zone seepage field at the fracture zone is calculated; according to the first data volume, Obtain the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points; fit the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points, and determine the fitting coefficient of the mudstone smear thickness and fault gouge ratio of the fault zone; Determine the mudstone smear sealing factor of the fracture zone of multiple sampling points according to the fitting coefficient and the fault mud ratio of the multiple sampling points of the mudstone smear type fault to be drilled; according to the mudstone smear seal of the
- the strength of the seepage field of the core fracture zone at the fracture zone and the mudstone smear sealing factor of the fracture zone at multiple sampling points were calculated at the same time, so as to analyze the mudstone smear type fault to be drilled for each sampling point. Closeness, compared with qualitative analysis, has high accuracy.
- a first core data volume containing pores is obtained, and the first core data volume containing pores is obtained by scanning a core model of a target interval containing faults; and the first core data volume containing pores Perform three-dimensional reconstruction of the volume to obtain a solid second core data volume; make a difference between the solid second core data volume and the first core data volume containing pores to obtain the pore data volume; obtain the core fracture zone seepage field according to the pore data volume ; Calculate the strength of the seepage field of the core fracture zone at the fracture zone; analyze the fault sealability according to the strength of the seepage field of the core fracture zone at the fracture zone.
- the process of analyzing the sealability of a fault is the process of quantitatively analyzing the sealability of a fault. Compared with qualitative analysis, the accuracy is high.
- Fig. 1 is a flowchart of a method for quantitative analysis of the sealability of a trap fault in an embodiment of the present invention
- FIG. 2 is a detailed flow chart of a method for quantitative analysis of the sealability of trap faults according to an embodiment of the present invention
- FIG. 3 is a front view of a physical simulation model of a trap to be drilled constructed in an embodiment of the present invention
- FIG. 4 is a schematic diagram of a cross-sectional view of a fault zone obtained by scanning in an embodiment of the present invention
- Figure 5 is a schematic diagram of a three-dimensional data volume in an embodiment of the present invention.
- Fig. 6 is a schematic diagram of a relationship chart in an embodiment of the present invention.
- FIG. 7 is a schematic diagram of a sealing quantitative analysis device for a trap fault in an embodiment of the present invention.
- Fig. 8 is a schematic diagram of a quantitative analysis system for the sealing of a trap fault in an embodiment of the present invention.
- FIG. 9 is a flowchart of a method for analyzing the sealability of an interrupted layer according to an embodiment of the present invention.
- FIG. 10 is a detailed flowchart of a fault sealing analysis method proposed by an embodiment of the present invention.
- Figure 11 is a schematic diagram of a first core data volume containing pores in an embodiment of the present invention.
- FIG. 12 is a schematic diagram of a solid second core data volume in an embodiment of the present invention.
- Figure 13 is a schematic diagram of a pore data volume in an embodiment of the present invention.
- Fig. 14 is a schematic diagram of a percolation communication channel in an embodiment of the present invention.
- Figure 15 is a schematic diagram of the seepage line in the embodiment of the present invention.
- 16 is a schematic diagram of the strength of the seepage field of the core fracture zone at the fracture zone in the embodiment of the present invention.
- FIG. 17 is a schematic diagram of a device for analyzing the sealability of an interrupted layer according to an embodiment of the present invention.
- FIG. 18 is a schematic diagram of a system for analyzing the sealability of an interrupted layer according to an embodiment of the present invention.
- Fig. 1 is a flowchart of a method for quantitative analysis of the sealability of a trap fault in an embodiment of the present invention. As shown in Fig. 1, the method includes:
- Step 101 Obtain a three-dimensional data volume of a mudstone smear-type fault in a trap to be drilled, and the three-dimensional data volume is obtained by scanning and reconstructing a fault simulation model of the mudstone smear-type fault in a trap to be drilled;
- Step 102 Obtain the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points according to the three-dimensional data volume;
- Step 103 Fitting the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points, and determine the fitting coefficient of the mudstone smear thickness of the fault zone and the fault gouge ratio;
- Step 104 Determine the mudstone smear sealing factor of the fault zone of the multiple sampling points according to the fitting coefficient and the fault gouge ratio of the multiple sampling points of the mudstone smear type fault to be drilled;
- Step 105 Analyze the sealability of each sampling point of the mudstone smear-type fault to be trapped by the trap according to the mudstone smear sealing factor of each sampling point.
- the sealability of the trap fault to be drilled can be analyzed, and in the analysis process, the fitting coefficient of the mudstone smear thickness of the fault zone and the fault mud ratio is determined, and the fitting coefficient can indicate the effective smearing degree of the mudstone. Therefore, the mudstone smear sealing factor of the fault zone is determined more accurately, which makes the final analysis of the sealability of each sampling point of the mudstone smear fault to be drilled more accurate.
- obtaining the three-dimensional data volume of the mudstone smear-type fault in the trap to be drilled includes:
- fault data and horizon data of the trap to be drilled determine the structural parameters and fault deformation stress of the physical simulation model of the trap to be drilled;
- the fault simulation model is obtained by loading the fault deformation stress on the physical simulation model of the trap to be drilled, and the physical simulation model of the trap to be drilled is constructed based on the structural parameters.
- the method before determining the structural parameters of the physical simulation model of the trap to be drilled, the method further includes:
- the fault data and horizon data of the trap to be drilled are obtained.
- the fault data and horizon data of the trap to be drilled can be obtained by interpretation in the seismic interpretation system or other related software (such as Geoeast software, Landmark software or Geoframe software) based on the post-stack seismic data. Obtain the section distribution law and distance information. According to the results of interpretation, sand and mudstone formations can be obtained, and the fault data and horizon data of the traps to be drilled provided by the interpretation provide accurate parameters and basis for physical simulation, ensuring the accuracy and science of the physical simulation model of drilling traps expected later. Sex.
- the structural parameters and fault deformation stress of the physical simulation model of the trap to be drilled can be determined. This process is also called the physical simulation of the trap to be drilled. Schematic design of the model.
- the structural parameters of the physical simulation model of the trap to be drilled include one or any combination of the model boundary, the physical simulation model similarity ratio, the physical simulation duration similarity ratio, the simulated formation material and the simulated mudstone material.
- the steps for constructing a physical simulation model of the trap to be drilled are as follows:
- S2 According to the principle of material similarity, determine the simulated formation material and the simulated mudstone layer material; for example, quartz sand is used to simulate sandstone layer and clay is used to simulate mudstone layer. These two types of materials are similar, relatively high and stable, and are currently commonly used.
- S3 Construct a physical simulation model of the trap to be drilled according to the structural parameters of the physical simulation model of the trap to be drilled.
- the physical simulation model of the trap to be drilled is loaded with the fault deformation stress to obtain the fault simulation model.
- the loaded fault deformation stress can be used in the deformed sandbox to cause the strata to move and form a fault.
- the mudstone layer descends through the hanging wall and drags smearing to form a mudstone smear layer, and finally a mudstone smear type fault.
- the mudstone thickness at the location where the mudstone smear loses continuity can determine the minimum mudstone thickness H of the fault seal by analyzing the mudstone smear disconnection, which is also called the fault zone mudstone smear thickness H.
- the fault deformation stress can be stopped, the structural deformation is completed, and the fault simulation model is obtained.
- a three-dimensional data volume of the mudstone smear-type fault in the trap to be drilled is obtained, and the three-dimensional data volume is obtained by scanning and reconstructing the fault simulation model of the mudstone smear-type fault in the trap to be drilled.
- Scanning methods such as industrial CT can be used to scan the tomographic simulation model.
- the dynamic monitoring scan position, scan frequency, scan interval, etc. can be set according to the accuracy requirements to ensure the scan accuracy.
- scanning the fault simulation model it is mainly to scan the fault zone in the fault simulation model to obtain the surface data of the mudstone smear type fault to be drilled. The more surface data obtained, the higher the accuracy of the final fault zone.
- the specific scanning process can be as follows:
- S1 Set different collection intervals according to the research accuracy requirements. Place the position to be scanned for the tomographic simulation model under the industrial CT, and scan the surface data.
- a difference method or a fitting method is used to perform three-dimensional reconstruction of the surface data to obtain a three-dimensional data volume of the mudstone smear type fault to be drilled.
- the three-dimensional data volume obtained by the three-dimensional reconstruction method has high accuracy.
- the three-dimensional data volume facilitates the comprehensive identification of fault simulation models.
- Reconstruct the surface data in three dimensions that is, construct unknown surface data (surface data between two sets of adjacent surface data) through the known surface data obtained by scanning, and finally recombine all surface data to form a three-dimensional Data body.
- the reconstruction method such as the difference can be realized with the help of software or other computer equipment.
- the difference method is used to reconstruct the surface data in three dimensions, that is, between the two sets of adjacent surface data (for example, two gray-scale scanned images), mathematical methods such as Gaussian difference are used to calculate the difference between the two sets of adjacent surface data.
- the 3D reconstruction of the surface data is carried out using the fitting method, that is, the process of fitting the shape of other surface data according to the change law of the surface data through a certain known surface data to obtain the unknown surface data.
- obtaining the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points according to the three-dimensional data volume includes:
- the fault gouge ratio of multiple sampling points is calculated.
- the attribute information of the fault zone includes the mudstone smear thickness of the fault zone at multiple sampling points, the fault distance of the fault, and the mud content of the section, while the fault distance of the fault and the mud content of the section at multiple sampling points.
- SGR is the ratio of fault gouge at each sampling point
- Z is the shale content of each sampling point
- in cm is the distance between each sampling point in cm.
- extracting the fault zone attribute information of the mudstone smear type fault from the three-dimensional data volume includes:
- the gray-scale characteristic value of the fault zone is generally obtained from the analysis of simulated materials.
- the gray-scale characteristic value of the fault zone is generally about 430.
- the fault zone of the mudstone smear-type fault can be identified .
- the identification of the fracture zone can be realized by means of gray-scale processing software such as VG, etc., of course, other methods can also be used, and the relevant changes should fall within the protection scope of the present invention.
- the gray-scale characteristic value of sandstone and mudstone is 436, that is, more than 436 is mudstone, and less than 436 is sandstone.
- the minimum mudstone thickness H of the fault closure can be determined by analyzing the mudstone smear disconnection, which is also called the fault zone mudstone smear thickness H. Plot the fracture zone mudstone smear thickness H and SGR on a chart to obtain the fault zone mudstone smear Thickness and SGR relationship chart, based on the relationship chart, determine the fitting coefficient between the mudstone smear thickness of the fault zone and the fault gouge ratio.
- the mudstone smear thickness and the fault gouge ratio of the fault zone at the multiple sampling points have a linear relationship.
- the fitting formula of the linear relationship may be as follows:
- H is the mudstone smear thickness of the fault zone
- A is the fitting coefficient.
- B is a constant value parameter.
- the mudstone smear thickness of the fault zone has a linear relationship with the fault gouge ratio, the larger the mudstone smear thickness of the fault zone, the larger the fault gouge ratio and the better the sealing.
- the corresponding SGR is 18%, that is, mudstone less than 0.1cm thick loses continuity is invalid smearing, and the corresponding SGR is less than 18% can not play an effective plugging effect. Therefore, using this method, the extreme point of continuous mudstone smearing can be accurately defined on the relationship chart, that is, the SGR lower limit of fault sealing.
- the following formula is used to determine the mudstone smear sealing factor of the fracture zone at multiple sampling points based on the fitting coefficient and the fault gouge ratios at multiple sampling points of the mudstone smear-type fault to be drilled:
- SGRN is the mudstone smear sealing factor in the fault zone of each sampling point
- A is the fitting coefficient
- C is a constant value parameter
- SGR is the fault gouge ratio at each sampling point.
- the physical meaning of the fitting coefficient is to indicate the effective degree of mudstone smearing, also known as the effective continuous mudstone smearing contribution rate, which is related to regional geological characteristics, and its value range is 0-1, more specifically,
- the resultant coefficient A is between 0.62 and 0.67.
- the mudstone smear sealing factor of the fault zone obtained after the calculation of the above formula takes into account the effective degree of mudstone smear and the effectiveness and heterogeneity of mudstone smear, rather than a rough simple Calculate the mudstone of all formations.
- the present invention can more accurately analyze the sealability of the mudstone smear-type fault to be drilled.
- the effective smearing degree of mudstone comes from quantitative physical simulation and analysis, which is suitable for the real situation under the set geological conditions, with high reliability, strong scientificity, and closer to the real underground core characteristics. In fact, according to different geological conditions in different blocks, the value of the fitting coefficient may be different.
- analyzing the sealing properties of each sampling point of the mudstone smear-type fault to be drilled in the trap to be drilled according to the mudstone smear sealing factor of the fracture zone of each sampling point including:
- the sampling point For each sampling point, if the mudstone smear sealing factor of the fault zone of the sampling point is greater than the threshold value, it is determined that the sampling point is closed; otherwise, the sampling point is not closed.
- the method further includes: determining whether to drill the trap to be drilled according to the sealing property of the mudstone smear type fault.
- the process of calculating the fitting coefficients according to this method can also be obtained by sampling from field geological outcrops, performing industrial CT scanning and then following similar steps. It should also belong to the protection scope of this application example.
- FIG. 2 is a detailed flow chart of the method for quantitative analysis of the sealing of trap faults according to an embodiment of the present invention.
- the detailed process of the method for quantitative analysis of the sealability of trap faults includes:
- Step 201 Obtain fault data and horizon data of the trap to be drilled according to the post-stack seismic data of the trap to be drilled;
- Step 202 Determine the structural parameters and fault deformation stress of the physical simulation model of the trap to be drilled according to the geological background information, fault data and horizon data of the trap to be drilled;
- Step 203 Obtain the surface data of the mudstone smear-type fault in the trap to be drilled, and perform three-dimensional reconstruction on the surface data to obtain a three-dimensional data volume of the mudstone smear-type fault in the trap to be drilled.
- the surface data is a simulation of the fault
- the model is obtained by scanning, and the fault simulation model is obtained by loading the fault deformation stress on the physical simulation model of the trap to be drilled;
- Step 204 Identify the fracture zone from the three-dimensional data volume according to the gray feature value of the fracture zone
- Step 205 Extract the attribute information of the fault zone body, the attribute information includes the mudstone smear thickness of the fault zone at multiple sampling points, the fault distance of the fault, and the shale content of the section;
- Step 206 Calculate the fault gouge ratio of the multiple sampling points according to the fault distance of the multiple sampling points and the mud content of the cross section;
- Step 207 Fit the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points, and determine the fitting coefficient of the mudstone smear thickness and fault gouge ratio of the fault zone;
- Step 208 Determine the mudstone smear sealing factor of the fault zone of the multiple sampling points according to the fitting coefficient and the fault gouge ratio of the multiple sampling points of the mudstone smear type fault to be drilled;
- Step 209 Determine the threshold value of the mudstone smear sealing factor of the fault zone to be drilled;
- Step 210 For each sampling point, if the mudstone smear sealing factor of the fracture zone of the sampling point is greater than the threshold value, it is determined that the sampling point is closed; otherwise, the sampling point is not closed.
- S11 Load the post-stack seismic data of the trap to be drilled into the seismic interpretation system or other related software, such as Geoeast software, and interpret the fault data and horizon data of the trap to be drilled.
- the size of the physical simulation model of the trap to be drilled is 48cm ⁇ 24cm ⁇ 26cm; the mudstone layer has 3 layers, The thickness is 1.5cm; the simulated displacement and deformation of the physical simulation model of the trap to be drilled is determined according to the fault distance value interpreted by the post-stack seismic data volume, that is, the section distance is the maximum displacement, which is 9cm in this example;
- quartz sand simulates a sandstone layer, and clay simulates a mudstone layer.
- Fig. 3 is a front view of a physical simulation model of a trap to be drilled constructed in an embodiment of the present invention.
- FIG. 4 is a schematic diagram of the fracture zone section obtained by scanning in the embodiment of the present invention, forming surface data. Here, 30 fracture zone sections are scanned and stored in DICOM format.
- the graphics workstation performs visualization processing and three-dimensional reconstruction.
- the professional three-dimensional reconstruction software Mimics is used to establish a 3D model by interpolation and perform editing processing to obtain a three-dimensional data volume.
- FIG. 5 is a schematic diagram of the three-dimensional data volume in an embodiment of the present invention.
- S15 Identify the fault zone from the three-dimensional data volume according to the gray feature value of the fault zone.
- the sand and mudstone are distinguished according to the gray feature value in the Mimics software.
- the gray characteristic value is 436, the ones larger than 436 are mudstone, and those smaller than 436 are sandstone.
- the attribute information of the fault zone body is extracted, and the attribute information includes the mudstone smear thickness of the fault zone at multiple sampling points, the distance of the fault, and the mud content of the section.
- the fault gouge ratio SGR of multiple sampling points is calculated; the obtained mudstone smear thickness value H and SGR are generated into a scatter plot to form the mudstone smear thickness value H and SGR relationship chart. Based on the relationship chart, the fitting coefficient of the mudstone smear thickness of the fault zone and the fault gouge ratio is determined through fitting.
- the fitting formula can be as follows:
- 0.6248 is the fitting coefficient, which is also the slope of the curve in the relationship chart.
- Figure 6 is a schematic diagram of the relationship chart in the embodiment of the present invention.
- the fitting coefficient 0.6248 represents the effective mudstone smearing degree, which is also called effective continuous mudstone smearing contribution rate.
- S17 Determine the threshold of the mudstone smear sealing factor of the fault zone to be drilled.
- the lower limit of the SGR for the fault sealability is 18%.
- the mudstone of the fault zone to be drilled is smeared with the threshold of the seal factor It is determined to be 18%.
- the threshold of the seal factor It is determined to be 18%.
- For each sampling point if the mudstone smear sealing factor of the fault zone of the sampling point is greater than the threshold, it is determined that the sampling point is closed; otherwise, the sampling point is not closed. So as to find out the weak points of sealing and avoid drilling risks.
- a three-dimensional data volume of the mudstone smeared fault to be drilled is obtained, and the three-dimensional data volume is a fault simulation model of the mudstone smeared fault to be drilled.
- Obtained by scanning reconstruction according to the three-dimensional data volume, obtain the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points; fit the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points to determine the fault Fitting coefficient of mudstone smear thickness and fault gouge ratio; according to the fitting coefficient and the fault gouge ratios of multiple sampling points of the mudstone smear type fault to be drilled, determine the mudstone smear closure of multiple sampling points in the fracture zone Factor; According to the mudstone smear sealing factor of each sampling point, analyze the sealability of each sampling point of the mudstone smear-type fault to be drilled.
- the present invention can analyze the sealability of the trap fault to be drilled, and in the analysis process, the fitting coefficient of the mudstone smear thickness of the fault zone and the fault mud ratio is determined, and the fitting coefficient can indicate the effective smearing degree of the mudstone. Therefore, the mudstone smear sealing factor of the fault zone is determined more accurately, which makes the final analysis of the sealability of each sampling point of the mudstone smear fault to be drilled more accurate.
- the present invention uses post-stack seismic data to improve the accuracy of the attribute information of mudstone smear-type faults, and uses advanced methods of physical simulation and industrial CT data quantitative acquisition and analysis to obtain three-dimensional data volume with high accuracy and quantitatively solve existing faults.
- the idealization of the calculation model the incomplete consideration of the factors, the scientificity and the low accuracy of the problem, provide a scientific basis for the closure evaluation, accurately evaluate the closure and screen out the leakage points, avoid the risk of drilling investment, and have a good
- the technical application prospects and economic benefits have improved the success rate of drilling and achieved the technical effect of evaluating fault block traps with higher accuracy based on fault sealing.
- an embodiment of the present invention also provides a sealing quantitative analysis device for trap faults, as described in the following embodiments. Since these problem-solving principles are similar to the quantitative analysis method of the sealing of the trap fault, the implementation of the device can refer to the implementation of the method, and the repetition will not be repeated.
- Fig. 7 is a schematic diagram of a device for quantitative analysis of the sealability of a trap fault in an embodiment of the present invention. As shown in Fig. 7, the device includes:
- the first module 701 is used to obtain a three-dimensional data volume of a mudstone smear-type fault in a trap to be drilled, and the three-dimensional data volume is obtained by scanning and reconstructing a fault simulation model of the mudstone smear-type fault in a trap to be drilled;
- the second module 702 is used to obtain the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points according to the three-dimensional data volume;
- the third module 703 is used to fit the mudstone smear thickness and the fault gouge ratio of the fault zone at multiple sampling points, and determine the fitting coefficient of the mudstone smear thickness and the fault gouge ratio of the fault zone;
- the fourth module 704 is used to determine the mudstone smear sealing factor of the fault zone of multiple sampling points according to the fitting coefficient and the fault gouge ratio of multiple sampling points of the mudstone smear type fault to be drilled;
- the fifth module 705 is used to analyze the sealability of each sampling point of the mudstone smear-type fault to be drilled according to the mudstone smear sealing factor of the fracture zone of each sampling point.
- the third module 703 is specifically used for:
- the second module 702 is specifically used for:
- the fault gouge ratio of multiple sampling points is calculated.
- the second module 702 is specifically used for:
- the first module 701 is specifically used for:
- fault data and horizon data of the trap to be drilled determine the structural parameters and fault deformation stress of the physical simulation model of the trap to be drilled;
- the fault simulation model is obtained by loading the fault deformation stress on the physical simulation model of the trap to be drilled, and the physical simulation model of the trap to be drilled is constructed based on the structural parameters.
- the first module 701 is specifically used for:
- the fault data and horizon data of the trap to be drilled are obtained.
- the structural parameters of the physical simulation model of the trap to be drilled include one or any combination of the model boundary, the physical simulation model similarity ratio, the physical simulation duration similarity ratio, the simulated formation material and the simulated mudstone material.
- the first module 701 is specifically used for:
- a difference method or a fitting method is used to perform three-dimensional reconstruction on the surface data to obtain a three-dimensional data volume of the mudstone smear type fault to be drilled.
- the fifth module 705 is specifically used for:
- the sampling point For each sampling point, if the mudstone smear sealing factor of the fault zone of the sampling point is greater than the threshold value, it is determined that the sampling point is closed; otherwise, the sampling point is not closed.
- a three-dimensional data volume of mudstone smeared faults to be drilled is obtained, and the three-dimensional data volume is a fault simulation model of the mudstone smeared faults to be drilled.
- Obtained by scanning reconstruction according to the three-dimensional data volume, obtain the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points; fit the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points to determine the fault Fitting coefficient of mudstone smear thickness and fault gouge ratio; according to the fitting coefficient and the fault gouge ratios of multiple sampling points of the mudstone smear type fault to be drilled, determine the mudstone smear closure of multiple sampling points in the fracture zone Factor; According to the mudstone smear sealing factor of each sampling point, analyze the sealability of each sampling point of the mudstone smear-type fault to be drilled.
- the present invention can analyze the sealability of the trap fault to be drilled, and in the analysis process, the fitting coefficient of the mudstone smear thickness of the fault zone and the fault mud ratio is determined, and the fitting coefficient can indicate the effective smearing degree of the mudstone. Therefore, the mudstone smear sealing factor of the fault zone is determined more accurately, which makes the final analysis of the sealability of each sampling point of the mudstone smear fault to be drilled more accurate.
- the present invention uses post-stack seismic data to improve the accuracy of the attribute information of mudstone smear-type faults, and uses advanced methods of physical simulation and industrial CT data quantitative acquisition and analysis to obtain three-dimensional data volume with high accuracy and quantitatively solve existing faults.
- the idealization of the calculation model the incomplete consideration of the factors, the scientificity and the low accuracy of the problem, provide a scientific basis for the closure evaluation, accurately evaluate the closure and screen out the leakage points, avoid the risk of drilling investment, and have a good
- the technical application prospects and economic benefits have improved the success rate of drilling and achieved the technical effect of evaluating fault block traps with higher accuracy based on fault sealing.
- FIG. 8 is a schematic diagram of the closedness quantitative analysis system for trap faults in an embodiment of the present invention, including: a fault simulation model construction unit 801 and a scanning unit 802 And the aforementioned trap and fault sealing quantitative analysis device 803, in which,
- the fault simulation model structure unit 801 is used to construct the fault simulation model of the mudstone smear type fault to be drilled;
- the scanning unit 802 is used to scan the fault simulation model to generate surface data of the mudstone smear type fault to be drilled and send to the sealing quantitative analysis device of the trap fault;
- the sealing quantitative analysis device 803 of the trap fault is used to reconstruct the received surface data to obtain a three-dimensional data volume; according to the three-dimensional data volume, obtain the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points; Fit the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points, and determine the fitting coefficient of the mudstone smear thickness and the fault gouge ratio in the fault zone; according to the fitting coefficient and the mudstone smear type fault to be drilled Determine the mudstone smear sealing factor of the fault zone at multiple sampling points, and analyze the mudstone smear seal factor of the fracture zone mudstone at each sampling point to analyze each sample of the mudstone smear type fault to be drilled Point of closure.
- a three-dimensional data volume of mudstone smeared faults to be drilled is obtained, and the three-dimensional data volume is a fault simulation model of the mudstone smeared faults to be drilled.
- Obtained by scanning reconstruction according to the three-dimensional data volume, obtain the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points; fit the mudstone smear thickness and fault gouge ratio of the fault zone at multiple sampling points to determine the fault Fitting coefficient of mudstone smear thickness and fault gouge ratio; according to the fitting coefficient and the fault gouge ratios of multiple sampling points of the mudstone smear type fault to be drilled, determine the mudstone smear closure of multiple sampling points in the fracture zone Factor; According to the mudstone smear sealing factor of each sampling point, analyze the sealability of each sampling point of the mudstone smear-type fault to be drilled.
- the present invention can analyze the sealability of the trap fault to be drilled, and in the analysis process, the fitting coefficient of the mudstone smear thickness of the fault zone and the fault mud ratio is determined, and the fitting coefficient can indicate the effective smearing degree of the mudstone. Therefore, the mudstone smear sealing factor of the fault zone is determined more accurately, which makes the final analysis of the sealability of each sampling point of the mudstone smear fault to be drilled more accurate.
- the present invention uses post-stack seismic data to improve the accuracy of the attribute information of mudstone smear-type faults, and uses advanced methods of physical simulation and industrial CT data quantitative acquisition and analysis to obtain three-dimensional data volume with high accuracy and quantitatively solve existing faults.
- the idealization of the calculation model the incomplete consideration of the factors, the scientificity and the low accuracy of the problem, provide a scientific basis for the closure evaluation, accurately evaluate the closure and screen out the leakage points, avoid the risk of drilling investment, and have a good
- the technical application prospects and economic benefits have improved the success rate of drilling and achieved the technical effect of evaluating fault block traps with higher accuracy based on fault sealing.
- FIG. 9 is a flowchart of the method for analyzing the sealability of a fault in an embodiment of the present invention. As shown in FIG. 9, the method includes:
- Step 901 Obtain a first core data volume containing pores, where the first core data volume containing pores is obtained by scanning a core model of a target interval containing faults;
- Step 902 Perform a three-dimensional reconstruction on the first core data volume containing pores to obtain a solid second core data volume
- Step 903 Make a difference between the second solid core data volume and the first core data volume containing pores to obtain a pore data volume
- Step 904 Obtain the seepage field of the core fracture zone according to the pore data volume
- Step 905 Calculate the strength of the core fracture zone seepage field at the fracture zone
- Step 906 According to the strength of the core fracture zone seepage field at the fracture zone, the fault sealing is analyzed.
- a first core data volume containing pores is obtained, and the first core data volume containing pores is obtained by scanning a core model of a target interval containing faults; and the first core data volume containing pores Perform three-dimensional reconstruction of the volume to obtain a solid second core data volume; make a difference between the solid second core data volume and the first core data volume containing pores to obtain the pore data volume; obtain the core fracture zone seepage field according to the pore data volume ; Calculate the strength of the seepage field of the core fracture zone at the fracture zone; analyze the fault sealability according to the strength of the seepage field of the core fracture zone at the fracture zone.
- the process of analyzing the sealability of a fault is the process of quantitatively analyzing the sealability of a fault. Compared with qualitative analysis, the accuracy is high.
- the target interval core model containing faults includes mudstone smear type faults, sand-shale butt-type faults and other types.
- the first core data volume containing pores is obtained by scanning the core model of the target interval containing faults. There can be many scanning methods. For example, high-precision industrial CT scanning can be used, and reasonable scanning steps and ranges can be set according to specific needs. , The scanning range must include the structure of the fracture zone to ensure a clear scan near the fracture zone. Of course, other three-dimensional scanning methods can also be used for scanning.
- the second solid core data volume does not contain voids.
- the difference between the second solid core data volume and the first core data volume containing pores is to obtain the pores in the core data volume, which is called the void data volume.
- the pore data volume the seepage field of the core fracture zone is obtained; the intensity of the seepage field of the core fracture zone at the fracture zone is calculated.
- the intensity is a specific quantitative value. According to the intensity, the sealing of the fault can be quantitatively analyzed.
- the first core data volume containing pores includes core skeleton volume data and/or pore volume data.
- performing three-dimensional reconstruction on the first core data volume containing pores to obtain a solid second core data volume includes:
- the pores in the first core data volume containing pores are reconstructed and filled to obtain a solid second core data volume.
- the first core data volume containing pores is a three-dimensional structure
- the three-dimensional reconstruction is to reconstruct and fill the pores to obtain a solid second core data volume.
- obtaining the seepage field of the core fracture zone according to the pore data volume includes:
- the connected pores and the pores that can be connected under the set pressure are formed to form multiple seepage connection channels perpendicular to the fault direction;
- the seepage field of the core fracture zone is constructed.
- the pore three-dimensional imaging data can be obtained by loading the pore data volume into professional imaging software, for example, VG or Simpleware software.
- the starting point of the seepage connection channel is generally on the side of the fault, and the end point of the seepage connection channel is the fault.
- N1, N2, N3... easy to analyze in subsequent steps.
- obtaining connected pores and pores capable of being connected under a set pressure according to the three-dimensional imaging data of the pores includes:
- the connected pores and the pores that can be connected under the set pressure are obtained, including:
- the pore size and the distance between the pores From the pore three-dimensional imaging data, find the connected pores and the pores that can be connected under the set pressure.
- the size of the gap and the distance between the pores can be set according to the principle of gap connection, and from the pore three-dimensional imaging data, the connected pores and the pores that can be connected under the set pressure can be found.
- calculating the strength of the core fracture zone seepage field at the fracture zone includes:
- the strength of the seepage field of the core fracture zone at the fracture zone is obtained.
- the percolation lines L1, L2, L3... of the percolation connected channel are obtained, wherein the density of the percolation line Indicates the strength of the seepage field.
- the direction of the seepage line from the starting point through the fault to the end point is the direction of the calibrated seepage field.
- the following formula is used to obtain the strength of the core fracture zone's seepage field at the fracture zone according to the number of seepage lines and the cross-sectional area of the fault:
- W is the strength of the core fracture zone's seepage field at the fracture zone
- L is the number of seepage lines
- S is the cross-sectional area of the fault.
- the core fracture zone seepage field is near the fracture zone where the strength coefficient W is stronger, the fault sealing is weaker, and the core fracture zone seepage field is faulted at the place where the fracture zone strength coefficient W is weaker or where there is no seepage line. The stronger the closure.
- the analysis of the fault sealability according to the strength of the core fracture zone's seepage field at the fracture zone includes:
- the intensity threshold is determined according to the drilling data.
- the strength W of the seepage field of the core fracture zone at the fracture zone is compared with the intensity threshold W1. Where W is greater than W1, it indicates that the seepage ability is strong, the fault is not strong against fluid, and the fluid is easy to pass, which confirms the core fracture.
- the seepage field of the zone is not closed at the fault zone.
- the point where W is less than W1 indicates that the seepage capacity is weak, and the fault has a strong ability to hinder fluids, and fluid is not easy to pass. It is determined that the seepage field of the core fracture zone is closed at the fault zone.
- it may be determined whether to drill the fault block trap under the control of the target fault according to the sealing properties of various places in the spatial section of the target fault.
- FIG. 10 is a detailed flowchart of the fault sealing analysis method proposed by the embodiment of the present invention. As shown in FIG. 10, a In the embodiment, the detailed process of the fault sealing analysis method includes:
- Step 1001 Obtain a first core data volume containing pores, where the first core data volume containing pores is obtained by scanning a core model of a target interval containing faults;
- Step 1002 reconstruct and fill the pores in the first core data volume containing pores to obtain a solid second core data volume
- Step 1003 making a difference between the second solid core data volume and the first core data volume containing pores to obtain a pore data volume
- Step 1004 Obtain pore three-dimensional imaging data of the pore data volume
- Step 1005 determine the pore size and the distance between the pores
- Step 1006 according to the pore size and the distance between the pores, from the pore three-dimensional imaging data, search for connected pores and pores that can be connected under a set pressure;
- Step 1007 Determine the start point and end point of the seepage connection channel
- Step 1008 according to the start point and the end point of the seepage communication channel, connect the pores and the pores that can be connected under the set pressure to form a plurality of seepage communication channels perpendicular to the fault direction;
- Step 1009 constructing the seepage field of the core fracture zone according to the multiple seepage communication channels
- Step 1010 Obtain the seepage lines corresponding to multiple seepage communication channels in the seepage field of the core fracture zone;
- Step 1011 Obtain the strength of the core fracture zone's seepage field at the fracture zone according to the number of seepage lines and the cross-sectional area of the fault;
- Step 1012 Compare the intensity of the seepage field at the core fracture zone with the intensity threshold of the core fracture zone seepage field. If the intensity of the core fracture zone seepage field at the fracture zone is greater than the intensity threshold, then the core fracture zone seepage The field is not closed at the fracture zone, otherwise, the core fracture zone seepage field is closed at the fracture zone; the intensity threshold is determined according to the drilling data.
- the target interval core containing the target fault A and target fault B with a length of 80mm and a core radius of 79mm.
- the core of the target interval is processed as necessary, and the length is 80mm.
- a columnar body with a width of 40mm and a height of 40mm forms a core model of the target interval, in which the external excess sandstone is cut off, so that the scanning rays can more easily penetrate the external surface to reach the internal fault.
- Preliminary observations have identified that the fault in the core of the target interval is a reverse fault, and the lithological docking method is sand and mud interconnection.
- the first core data volume with pores obtained by scanning contains the core skeleton volume and pore volume.
- the data near the fault zone is characterized by tight mudstone Floor.
- Figure 11 is a schematic diagram of the first core data volume with pores in the embodiment of the present invention.
- the first core data volume with pores is reconstructed by professional three-dimensional reconstruction VG software to obtain a solid second core data volume.
- Figure 12 shows A schematic diagram of the second solid core data volume in the embodiment of the present invention.
- the second solid core data volume and the first core data volume containing pores are compared to obtain the pore data volume.
- Figure 13 shows the pore data in the embodiment of the present invention. Schematic diagram of the body.
- FIG. 14 is a schematic diagram of the seepage connection channel in the embodiment of the present invention.
- Figure 15 is a schematic diagram of the seepage line in the embodiment of the present invention.
- the density of the seepage line indicates the strength of the seepage field, the direction of the seepage line from the starting point through the fault to the end point is the direction of the calibrated seepage field, and the density of the seepage line represents The strength of the seepage field, the direction of the seepage line passing through the fault from the starting point to the end point is the direction of the calibrated seepage field.
- Figure 16 It is a schematic diagram of the strength of the seepage field of the core fracture zone at the fracture zone in the embodiment of the present invention.
- the strength threshold of the core fracture zone seepage field is determined to be 2.8, and the intensity of the core fracture zone seepage field at the fracture zone is compared with the intensity threshold. If the core fracture zone seepage field strength at the fracture zone is greater than all If the strength threshold is stated, the core fracture zone seepage field is not closed at the fracture zone. Otherwise, the core fracture zone seepage field is closed at the fracture zone.
- Table 1 shows the analysis of the core fracture zone seepage field closed at the fracture zone in the example of the present invention. As a result, it can be seen from Table 1 that the seepage field of the core fracture zone is closed at 31 fracture zones.
- a first core data volume containing pores is obtained, and the first core data volume containing pores is obtained by scanning a core model of a target interval containing faults; Perform a three-dimensional reconstruction of a core data volume to obtain a solid second core data volume; make a difference between the solid second core data volume and the first core data volume containing pores to obtain a pore data volume; obtain a core fracture according to the pore data volume Zone seepage field; calculate the strength of the core fracture zone seepage field at the fracture zone; analyze the fault sealing according to the strength of the core fracture zone seepage field at the fracture zone.
- the process of analyzing the sealability of a fault is the process of quantitatively analyzing the sealability of a fault. Compared with qualitative analysis, the accuracy is high.
- an embodiment of the present invention also provides a fault sealing analysis device, as described in the following embodiment. Since these problem-solving principles are similar to the fault sealing analysis method, the implementation of the device can refer to the implementation of the method, and the repetition will not be repeated.
- FIG. 17 is a schematic diagram of a device for analyzing the sealability of an interrupted layer according to an embodiment of the present invention. As shown in FIG. 17, the device includes:
- the first data obtaining module 1701 is used to obtain a first core data volume containing pores, and the first core data volume containing pores is obtained by scanning a core model of a target interval containing faults;
- the second data obtaining module 1702 is used to perform three-dimensional reconstruction of the first core data volume containing pores to obtain a solid second core data volume;
- the third data obtaining module 1703 is used to make a difference between the second solid core data volume and the first core data volume containing pores to obtain a pore data volume;
- the core fracture zone seepage field acquisition module 1704 is used to obtain the core fracture zone seepage field according to the pore data volume;
- Calculation module 1705 used to calculate the strength of the core fracture zone seepage field at the fracture zone
- the analysis module 1706 is used to analyze the fault sealability according to the strength of the core fracture zone's seepage field at the fracture zone.
- the second data obtaining module 1702 is specifically configured to:
- the pores in the first core data volume containing pores are reconstructed and filled to obtain a solid second core data volume.
- the core fracture zone seepage field obtaining module 1704 is specifically used for:
- the connected pores and the pores that can be connected under the set pressure are formed to form multiple seepage connection channels perpendicular to the fault direction;
- the seepage field of the core fracture zone is constructed.
- the core fracture zone seepage field obtaining module 1704 is specifically used for:
- the pore size and the distance between the pores From the pore three-dimensional imaging data, find the connected pores and the pores that can be connected under the set pressure.
- calculation module 1705 is specifically configured to:
- the strength of the seepage field of the core fracture zone at the fracture zone is obtained.
- the analysis module 1706 is specifically configured to:
- the core fracture zone seepage field is not closed at the fracture zone. Otherwise, the seepage field of the core fracture zone is closed at the fracture zone, and the intensity threshold is determined according to the drilling data.
- calculation module 1705 is specifically configured to:
- W is the strength of the core fracture zone's seepage field at the fracture zone
- L is the number of seepage lines
- S is the cross-sectional area of the fault.
- the first core data volume containing pores includes core skeleton volume data and/or pore volume data.
- a first core data volume containing pores is obtained, and the first core data volume containing pores is obtained by scanning a core model of a target interval containing faults; Perform a three-dimensional reconstruction of a core data volume to obtain a solid second core data volume; make a difference between the solid second core data volume and the first core data volume containing pores to obtain a pore data volume; obtain a core fracture according to the pore data volume Zone seepage field; calculate the strength of the core fracture zone seepage field at the fracture zone; analyze the fault sealing according to the strength of the core fracture zone seepage field at the fracture zone.
- the process of analyzing the sealability of a fault is the process of quantitatively analyzing the sealability of a fault. Compared with qualitative analysis, the accuracy is high.
- FIG. 18 is a schematic diagram of the fault sealing analysis system according to the embodiment of the present invention.
- the system includes:
- the second scanning unit 1802 is configured to:
- the first core data volume containing pores is sent to the first data obtaining module 1801.
- a first core data volume containing pores is obtained, and the first core data volume containing pores is obtained by scanning the core model of the target interval containing faults; Perform a three-dimensional reconstruction of a core data volume to obtain a solid second core data volume; make a difference between the solid second core data volume and the first core data volume containing pores to obtain a pore data volume; obtain a core fracture according to the pore data volume Zone seepage field; calculate the strength of the core fracture zone seepage field at the fracture zone; analyze the fault sealing according to the strength of the core fracture zone seepage field at the fracture zone.
- the process of analyzing the sealability of a fault is the process of quantitatively analyzing the sealability of a fault. Compared with qualitative analysis, the accuracy is high.
- the embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.
- a computer-usable storage media including but not limited to disk storage, CD-ROM, optical storage, etc.
- These computer program instructions can also be stored in a computer-readable memory that can guide a computer or other programmable data processing equipment to work in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture including the instruction device.
- the device implements the functions specified in one process or multiple processes in the flowchart and/or one block or multiple blocks in the block diagram.
- These computer program instructions can also be loaded on a computer or other programmable data processing equipment, so that a series of operation steps are executed on the computer or other programmable equipment to produce computer-implemented processing, so as to execute on the computer or other programmable equipment.
- the instructions provide steps for implementing functions specified in a flow or multiple flows in the flowchart and/or a block or multiple blocks in the block diagram.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
Description
Claims (22)
- 一种圈闭断层的封闭性定量分析方法,其特征在于,包括:获得待钻圈闭的泥岩涂抹型断层的三维数据体,所述三维数据体是对待钻圈闭的泥岩涂抹型断层的断层模拟模型进行扫描重构获得的;根据所述三维数据体,获得多个采样点的断裂带泥岩涂抹厚度和断层泥比率;对多个采样点的断裂带泥岩涂抹厚度和断层泥比率进行拟合,确定断裂带泥岩涂抹厚度与断层泥比率的拟合系数;根据所述拟合系数和待钻圈闭的泥岩涂抹型断层的多个采样点的断层泥比率,确定多个采样点的断裂带泥岩涂抹封闭因子;根据每个采样点的断裂带泥岩涂抹封闭因子,分析待钻圈闭的泥岩涂抹型断层的每个采样点的封闭性。
- 如权利要求1所述的圈闭断层的封闭性定量分析方法,其特征在于,还包括:在进行多个采样点的断裂带泥岩涂抹厚度和断层泥比率的拟合时,多个采样点的断裂带泥岩涂抹厚度和断层泥比率为线性关系。
- 如权利要求1所述的圈闭断层的封闭性定量分析方法,其特征在于,根据所述三维数据体,获得多个采样点的断裂带泥岩涂抹厚度和断层泥比率,包括:从所述三维数据体中,提取泥岩涂抹型断层的断裂带属性信息,所述属性信息包括多个采样点的断裂带泥岩涂抹厚度、断层的断距和断面的泥质含量;根据多个采样点的断层的断距和断面的泥质含量,计算多个采样点的断层泥比率。
- 如权利要求3所述的圈闭断层的封闭性定量分析方法,其特征在于,从所述三维数据体中,提取泥岩涂抹型断层的断裂带属性信息,包括:根据断裂带体的灰度特征值,从所述三维数据体中,识别出断裂带体;提取出断裂带体的属性信息。
- 如权利要求1所述的圈闭断层的封闭性定量分析方法,其特征在于,获得待钻圈闭的泥岩涂抹型断层的三维数据体,包括:根据待钻圈闭的地质背景信息、断层数据和层位数据,确定待钻圈闭物理模拟模型的构造参数、断层变形应力;获得待钻圈闭的泥岩涂抹型断层的面数据,对所述面数据进行三维重构,获得待钻圈闭的泥岩涂抹型断层的三维数据体,所述面数据是对断层模拟模型进行扫描获得的, 所述断层模拟模型是对待钻圈闭物理模拟模型加载断层变形应力获得的,所述待钻圈闭物理模拟模型是基于所述构造参数构造的。
- 如权利要求5所述的圈闭断层的封闭性定量分析方法,其特征在于,待钻圈闭物理模拟模型的构造参数包括模型边界、物理模拟模型相似比、物理模拟时长相似比、模拟地层材料和模拟泥岩层材料中的其中一种或任意组合。
- 如权利要求5所述的圈闭断层的封闭性定量分析方法,其特征在于,采用差值法或拟合法对所述面数据进行三维重构,获得待钻圈闭的泥岩涂抹型断层的三维数据体。
- 如权利要求1所述的圈闭断层的封闭性定量分析方法,其特征在于,根据每个采样点的断裂带泥岩涂抹封闭因子,分析待钻圈闭的泥岩涂抹型断层的每个采样点的封闭性,包括:确定待钻圈闭的断裂带泥岩涂抹封闭因子的阈值;对每个采样点,若该采样点的断裂带泥岩涂抹封闭因子大于所述阈值,确定该采样点的封闭;否则,该采样点不封闭。
- 一种圈闭断层的封闭性定量分析装置,其特征在于,包括:第一模块,用于获得待钻圈闭的泥岩涂抹型断层的三维数据体,所述三维数据体是对待钻圈闭的泥岩涂抹型断层的断层模拟模型进行扫描重构获得的;第二模块,用于根据所述三维数据体,获得多个采样点的断裂带泥岩涂抹厚度和断层泥比率;第三模块,用于对多个采样点的断裂带泥岩涂抹厚度和断层泥比率进行拟合,确定断裂带泥岩涂抹厚度与断层泥比率的拟合系数;第四模块,用于根据所述拟合系数和待钻圈闭的泥岩涂抹型断层的多个采样点的断层泥比率,确定多个采样点的断裂带泥岩涂抹封闭因子;第五模块,用于根据每个采样点的断裂带泥岩涂抹封闭因子,分析待钻圈闭的泥岩涂抹型断层的每个采样点的封闭性。
- 一种圈闭断层的封闭性定量分析***,其特征在于,包括:断层模拟模型构造单元、扫描单元和权利要求9所述的圈闭断层的封闭性定量分析装置,其中,断层模拟模型构造单元,用于构造对待钻圈闭的泥岩涂抹型断层的断层模拟模型;扫描单元,用于对断层模拟模型进行扫描,生成待钻圈闭的泥岩涂抹型断层的面数据,并发送至圈闭断层的封闭性定量分析装置;圈闭断层的封闭性定量分析装置,用于对接收的面数据进行重构,获得三维数据体;根据所述三维数据体,获得多个采样点的断裂带泥岩涂抹厚度和断层泥比率;对多个采样点的断裂带泥岩涂抹厚度和断层泥比率进行拟合,确定断裂带泥岩涂抹厚度与断层泥比率的拟合系数;根据所述拟合系数和待钻圈闭的泥岩涂抹型断层的多个采样点的断层泥比率,确定多个采样点的断裂带泥岩涂抹封闭因子;根据每个采样点的断裂带泥岩涂抹封闭因子,分析待钻圈闭的泥岩涂抹型断层的每个采样点的封闭性。
- 一种断层封闭性分析方法,其特征在于,包括:获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体;对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;根据孔隙数据体,获得岩心断裂带渗流场;计算岩心断裂带渗流场在断裂带处的强度;根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性。
- 如权利要求11所述的断层封闭性分析方法,其特征在于,对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体,包括:对含有孔隙的第一岩心数据体中的孔隙进行重建填充,获得实心的第二岩心数据体。
- 如权利要求11所述的断层封闭性分析方法,其特征在于,根据孔隙数据体,获得岩心断裂带渗流场,包括:获得孔隙数据体的孔隙三维成像数据;根据孔隙三维成像数据,获得联通孔隙和能够在设定压力下形成联通的孔隙;确定渗流联通通道的起点和终点;根据渗流联通通道的起点和终点,联通孔隙和能够在设定压力下形成联通的孔隙,形成多个垂直通过断层方向的渗流联通通道;根据所述多个渗流联通通道,构造岩心断裂带渗流场。
- 如权利要求13所述的断层封闭性分析方法,其特征在于,根据孔隙三维成像数据,获得联通孔隙和能够在设定压力下形成联通的孔隙,包括:确定孔隙尺寸和孔隙间的距离;根据孔隙尺寸和孔隙间的距离,从孔隙三维成像数据中,查找联通孔隙和能够在设定压力下形成联通的孔隙。
- 如权利要求13所述的断层封闭性分析方法,其特征在于,计算岩心断裂带渗流场在断裂带处的强度,包括:获得岩心断裂带渗流场中多个渗流联通通道对应的渗流线;根据渗流线的数量和断层的断面面积,获得岩心断裂带渗流场在断裂带处的强度。
- 如权利要求15所述的断层封闭性分析方法,其特征在于,采用如下公式,根据渗流线的数量和断层的断面面积,获得岩心断裂带渗流场在断裂带处的强度:W=L/SW为岩心断裂带渗流场在断裂带处的强度;L为渗流线的数量;S为断层的断面面积。
- 如权利要求11所述的断层封闭性分析方法,其特征在于,根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性,包括:将岩心断裂带渗流场在断裂带处的强度与岩心断裂带渗流场的强度阈值进行比较,若岩心断裂带渗流场在断裂带处的强度大于所述强度阈值,则岩心断裂带渗流场在断裂带处不封闭,否则,岩心断裂带渗流场在断裂带处封闭;所述强度阈值是根据钻井数据确定的。
- 如权利要求11所述的断层封闭性分析方法,其特征在于,含有孔隙的第一岩心数据体包括岩心骨架体积数据和/或孔隙体积数据。
- 一种断层封闭性分析装置,其特征在于,包括:第一数据获得模块,用于获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;第二数据获得模块,用于对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体;第三数据获得模块,用于对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;岩心断裂带渗流场获得模块,用于根据孔隙数据体,获得岩心断裂带渗流场;计算模块,用于计算岩心断裂带渗流场在断裂带处的强度;分析模块,用于根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性。
- 一种断层封闭性分析***,其特征在于,包括:权利要求19所述的断层封闭性分析装置,第二扫描单元,其中,所述第二扫描单元,用于:对含有断层的目标层段岩心模型进行扫描,获得含有孔隙的第一岩心数据体;将含有孔隙的第一岩心数据体发送至第一数据获得模块。
- 一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至8、权利要求11至18任一项所述方法。
- 一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至8、权利要求11至18任一项所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GB2112593.5A GB2595809B (en) | 2019-07-18 | 2020-08-26 | Method, apparatus, and system for quantitative analysis of seal of trap faults |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910649401.2 | 2019-07-18 | ||
CN201910649401.2A CN112329187B (zh) | 2019-07-18 | 2019-07-18 | 断层封闭性分析方法、装置及*** |
CN202010303957.9A CN113534248A (zh) | 2020-04-17 | 2020-04-17 | 圈闭断层的封闭性定量分析方法、装置及*** |
CN202010303957.9 | 2020-04-17 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2021008630A1 true WO2021008630A1 (zh) | 2021-01-21 |
Family
ID=74210189
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2020/111223 WO2021008630A1 (zh) | 2019-07-18 | 2020-08-26 | 圈闭断层的封闭性定量分析方法、装置及*** |
Country Status (2)
Country | Link |
---|---|
GB (1) | GB2595809B (zh) |
WO (1) | WO2021008630A1 (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112946245A (zh) * | 2021-03-11 | 2021-06-11 | 中国电建集团成都勘测设计研究院有限公司 | 红层地区软弱夹层的强度参数确定方法及装置 |
CN112965108A (zh) * | 2021-02-05 | 2021-06-15 | 中国石油天然气股份有限公司 | 圈闭盖层垂向封闭性确定方法及*** |
CN113466931A (zh) * | 2021-05-31 | 2021-10-01 | 中国石油天然气股份有限公司 | 一种地质时期断层封闭性评价方法及*** |
CN113687441A (zh) * | 2021-07-28 | 2021-11-23 | 中国海洋石油集团有限公司 | 一种陆相断陷盆地断裂带的断层封闭能力定量评价方法 |
CN113870676A (zh) * | 2021-08-31 | 2021-12-31 | 中国石油大学(北京) | 一种断层模拟装置及模拟方法 |
CN114004104A (zh) * | 2021-11-09 | 2022-02-01 | 中铁第一勘察设计院集团有限公司 | 基于棋盘格测试的cors站点选取方法 |
CN114994282A (zh) * | 2022-06-01 | 2022-09-02 | 国家石油天然气管网集团有限公司 | 一种面向交变应力环境的断层滑移判别方法 |
CN116146202A (zh) * | 2023-03-30 | 2023-05-23 | 中国石油大学(华东) | 一种礁滩碳酸盐岩的圈闭划分方法和装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7840914B1 (en) * | 2005-05-13 | 2010-11-23 | Massachusetts Institute Of Technology | Distributing computations in a parallel processing environment |
CN106772675A (zh) * | 2016-12-05 | 2017-05-31 | 中国石油化工股份有限公司 | 基于断层断裂结构评价断层启闭性的方法 |
CN106918839A (zh) * | 2015-12-25 | 2017-07-04 | 中国石油天然气股份有限公司 | 构造圈闭模型的定量评价方法及装置 |
CN107292435A (zh) * | 2017-06-16 | 2017-10-24 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种基于含油高度定量预测的断块圈闭评价方法 |
CN109239778A (zh) * | 2018-09-14 | 2019-01-18 | 中国石油大学(华东) | 一种断层侧向封闭性定量评价方法 |
-
2020
- 2020-08-26 WO PCT/CN2020/111223 patent/WO2021008630A1/zh active Application Filing
- 2020-08-26 GB GB2112593.5A patent/GB2595809B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7840914B1 (en) * | 2005-05-13 | 2010-11-23 | Massachusetts Institute Of Technology | Distributing computations in a parallel processing environment |
CN106918839A (zh) * | 2015-12-25 | 2017-07-04 | 中国石油天然气股份有限公司 | 构造圈闭模型的定量评价方法及装置 |
CN106772675A (zh) * | 2016-12-05 | 2017-05-31 | 中国石油化工股份有限公司 | 基于断层断裂结构评价断层启闭性的方法 |
CN107292435A (zh) * | 2017-06-16 | 2017-10-24 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种基于含油高度定量预测的断块圈闭评价方法 |
CN109239778A (zh) * | 2018-09-14 | 2019-01-18 | 中国石油大学(华东) | 一种断层侧向封闭性定量评价方法 |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112965108B (zh) * | 2021-02-05 | 2023-08-22 | 中国石油天然气股份有限公司 | 圈闭盖层垂向封闭性确定方法及*** |
CN112965108A (zh) * | 2021-02-05 | 2021-06-15 | 中国石油天然气股份有限公司 | 圈闭盖层垂向封闭性确定方法及*** |
CN112946245A (zh) * | 2021-03-11 | 2021-06-11 | 中国电建集团成都勘测设计研究院有限公司 | 红层地区软弱夹层的强度参数确定方法及装置 |
CN113466931A (zh) * | 2021-05-31 | 2021-10-01 | 中国石油天然气股份有限公司 | 一种地质时期断层封闭性评价方法及*** |
CN113466931B (zh) * | 2021-05-31 | 2023-08-22 | 中国石油天然气股份有限公司 | 一种地质时期断层封闭性评价方法及*** |
CN113687441A (zh) * | 2021-07-28 | 2021-11-23 | 中国海洋石油集团有限公司 | 一种陆相断陷盆地断裂带的断层封闭能力定量评价方法 |
CN113687441B (zh) * | 2021-07-28 | 2022-11-08 | 中国海洋石油集团有限公司 | 一种陆相断陷盆地断裂带的断层封闭能力定量评价方法 |
CN113870676A (zh) * | 2021-08-31 | 2021-12-31 | 中国石油大学(北京) | 一种断层模拟装置及模拟方法 |
CN113870676B (zh) * | 2021-08-31 | 2022-11-11 | 中国石油大学(北京) | 一种断层模拟装置及模拟方法 |
CN114004104A (zh) * | 2021-11-09 | 2022-02-01 | 中铁第一勘察设计院集团有限公司 | 基于棋盘格测试的cors站点选取方法 |
CN114004104B (zh) * | 2021-11-09 | 2024-06-07 | 中铁第一勘察设计院集团有限公司 | 基于棋盘格测试的cors站点选取方法 |
CN114994282A (zh) * | 2022-06-01 | 2022-09-02 | 国家石油天然气管网集团有限公司 | 一种面向交变应力环境的断层滑移判别方法 |
CN116146202B (zh) * | 2023-03-30 | 2023-10-13 | 中国石油大学(华东) | 一种礁滩碳酸盐岩的圈闭划分方法和装置 |
CN116146202A (zh) * | 2023-03-30 | 2023-05-23 | 中国石油大学(华东) | 一种礁滩碳酸盐岩的圈闭划分方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
GB202112593D0 (en) | 2021-10-20 |
GB2595809B (en) | 2023-05-24 |
GB2595809A (en) | 2021-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2021008630A1 (zh) | 圈闭断层的封闭性定量分析方法、装置及*** | |
RU2651719C1 (ru) | Способ усовершенствованного планирования высоты трещины гидроразрыва в слоистой породе подземного пласта | |
Iding et al. | Evaluating the impact of fractures on the long-term performance of the In Salah CO2 storage site | |
CN103258091B (zh) | 非常规油气藏水平井段三维岩体力学模型建立的方法及装置 | |
Espinoza et al. | Natural and induced fractures in coal cores imaged through X-ray computed microtomography—Impact on desorption time | |
Al-Marzouqi | Digital rock physics: Using CT scans to compute rock properties | |
CN104948176B (zh) | 一种基于渗透增大率识别碳酸盐岩储层裂缝的方法 | |
RU2601733C2 (ru) | Способ построения геолого-гидродинамических моделей двойной среды залежей баженовской свиты | |
CN107558999B (zh) | 一种孤立缝洞体剩余油定量预测方法 | |
CN107795320B (zh) | 一种水平井碳酸盐岩储层参数的计算方法 | |
CN104977617A (zh) | 储层裂缝识别方法和成像测井储层裂缝识别方法 | |
CN110632657B (zh) | 泥岩涂抹型断层封闭性分析方法及装置 | |
CN110632658B (zh) | 断层的侧向封闭性分析方法及装置 | |
Wennberg et al. | A brief introduction to the use of X-ray computed tomography (CT) for analysis of natural deformation structures in reservoir rocks | |
CN110632659A (zh) | 断层空间封闭性分析方法、装置及*** | |
CN107831540A (zh) | 储层物性参数直接提取新方法 | |
Huang et al. | Modeling well interference and optimal well spacing in unconventional reservoirs using the fast marching method | |
CN114690243A (zh) | 一种碎屑岩断层古封闭性三维数值模拟恢复方法及装置 | |
Yang et al. | Study on the stimulation effectiveness evaluation of large-scale hydraulic fracturing simulation experiment based on optical scanning technology | |
Hadavand et al. | A practical methodology for integration of 4D seismic in steam-assisted-gravity-drainage reservoir characterization | |
He et al. | Self-supporting conductivity in shallow and ultra shallow shale reservoirs: Under hydraulic, CO2 and Sc-CO2 fracturing | |
CN107942381B (zh) | 一种致密油储层层理缝定量预测方法 | |
CN113534248A (zh) | 圈闭断层的封闭性定量分析方法、装置及*** | |
CN114114411B (zh) | 三维断层疏导性的定量确定方法及装置 | |
Fossen | Structural core observations in a siliciclastic reservoir-scale framework |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 20840334 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 202112593 Country of ref document: GB Kind code of ref document: A Free format text: PCT FILING DATE = 20200826 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 20840334 Country of ref document: EP Kind code of ref document: A1 |