CN112305609B - Method for quantitatively characterizing fracture mechanical property based on structural style analysis - Google Patents
Method for quantitatively characterizing fracture mechanical property based on structural style analysis Download PDFInfo
- Publication number
- CN112305609B CN112305609B CN202011197303.9A CN202011197303A CN112305609B CN 112305609 B CN112305609 B CN 112305609B CN 202011197303 A CN202011197303 A CN 202011197303A CN 112305609 B CN112305609 B CN 112305609B
- Authority
- CN
- China
- Prior art keywords
- fracture
- construction
- degrees
- mechanical property
- types
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000004458 analytical method Methods 0.000 title claims abstract description 15
- 238000010276 construction Methods 0.000 claims abstract description 58
- 238000011158 quantitative evaluation Methods 0.000 claims abstract description 13
- 230000001427 coherent effect Effects 0.000 claims abstract description 12
- 238000013178 mathematical model Methods 0.000 claims abstract description 8
- 238000011156 evaluation Methods 0.000 claims abstract description 6
- 238000006243 chemical reaction Methods 0.000 claims abstract description 5
- 206010017076 Fracture Diseases 0.000 claims description 204
- 208000010392 Bone Fractures Diseases 0.000 claims description 197
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000011160 research Methods 0.000 claims description 10
- 230000000694 effects Effects 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000010586 diagram Methods 0.000 claims description 8
- 230000007480 spreading Effects 0.000 claims description 8
- 238000003892 spreading Methods 0.000 claims description 8
- 230000002159 abnormal effect Effects 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000013507 mapping Methods 0.000 claims description 4
- 241001415846 Procellariidae Species 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 2
- 238000005065 mining Methods 0.000 claims description 2
- 239000013598 vector Substances 0.000 claims description 2
- 238000001125 extrusion Methods 0.000 claims 2
- 230000035772 mutation Effects 0.000 claims 2
- 238000011161 development Methods 0.000 abstract description 14
- 230000018109 developmental process Effects 0.000 description 13
- 230000009471 action Effects 0.000 description 8
- 230000008859 change Effects 0.000 description 6
- 208000035126 Facies Diseases 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000010008 shearing Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 210000001503 joint Anatomy 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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/62—Physical property of subsurface
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
A method for quantitatively characterizing fracture mechanical properties based on structural style analysis comprises the following steps: reading fracture-related structural style characteristics based on the three-dimensional seismic data, the coherent data volume and the curvature data volume data; based on the difference characteristics of fracture under the background of regional stress with different properties, selecting five types of construction styles capable of reflecting fracture mechanical properties, and establishing a detailed construction style database; solving the weight coefficients and membership degrees of the five types of construction patterns by establishing membership functions and evaluation matrixes, and preliminarily realizing the conversion of geological features to quantitative data; establishing a quantitative evaluation mathematical model based on the weight coefficients of the five types of construction styles and the membership degree of the single construction style, and solving the numerical value of the mechanical property of each breakpoint; and drawing a plane development graph of the fracture mechanical property value of each trunk by using Petrel software based on the mechanical property value of each breakpoint to finish the division of the quantitative evaluation standard of the fracture mechanical property of the land-phase fault-trap basin. The method can efficiently, accurately and quantitatively depict the fracture property.
Description
Technical Field
The invention relates to a method for analyzing and quantitatively characterizing the mechanical property of a fracture zone. In particular to a method for quantitatively characterizing fracture mechanical properties based on structural style analysis.
Background
The development and evolution of the eastern China new continental facies collapse oily basin are wholly controlled by the double tectonic actions of horizontal stretching generated by upwelling of the valance and sliding shearing generated by the activity of a tancotta fracture zone (Maapricot wall, 1983; Xujiawei, 1995; Lacqueradex, 2014; and the like), so that the fracture in the basin is very developed and the fracture properties are complex and various, and the complex fracture control action is always a hotspot and a difficult point of research.
In recent years, with the continuous deepening of exploration and development research of Bohai sea oil fields, exploration practices show that fracture and reservoir control effects with different mechanical properties have obvious differences, and the specific expression is as follows: the fracture with strong sliding property has good lateral plugging effect, and even if the two disks are fractured and sand butt joint occurs, the trap can still be plugged and hidden; the fracture with strong stretching property has good oil-gas vertical transportation and conduction effects, but the lateral plugging property is poor, and the fracture two disks are difficult to plug and store once sand abutting and trapping occurs. Therefore, the fracture mechanical property is represented in a fine analysis and quantification mode, and the method has important significance for searching a favorable exploration target in an exploration stage and improving the water injection development efficiency in a development stage.
Fracture is a structure formed by that when the ground stress of a rock stratum reaches a certain degree, the continuous integrity of the rock stratum is damaged, fracture occurs, and obvious displacement occurs along a fracture surface (belt), and the mechanical property of fracture refers to the property of the ground stress in the fracture development process. The existing research method mainly aims at qualitative evaluation, and mainly has a rough analysis on fracture mechanical properties by analyzing the ground stress with different properties generated by plate motion in each property period and combining with the occurrence characteristics of fracture. However, the ground stress generated by the movement of the plate is large in scale, and the release of the stress of the target fracture specifically applied to exploration practice also has a certain difference, the ground stress is not tensile stress, torsional stress (shearing stress) or compressive stress, but is not single ground stress in the fracture development process, more is composite superposition between two or even three stresses, and the judgment of the fracture mechanical property is more complicated.
Taking the newly-generated fracture of the east China fractured-depressed basin as an example, based on the dual dynamic background of sliding and shearing, the mechanical properties of the fracture are qualitatively classified into three types of tension, tension (or torsion) and sliding according to fracture characteristics by the predecessors (Qijiafu, 2004; Yuyixin, 2012; Wuzhiping, 2013; and the like), but the classification scheme does not clearly define the classification basis and standard, and is easy to change due to the change of subjective recognition of the predecessors in the application process. The qualitative knowledge of the "ambiguous" cannot meet the increasingly refined exploration and development requirements, so that the quantitative representation of the qualitative geological knowledge of fracture mechanical properties is an important direction for the development of the subsequent exploration and development work.
Disclosure of Invention
The invention aims to solve the technical problem of providing a method for analyzing and quantitatively characterizing fracture mechanical properties based on a structural style, which can efficiently, accurately and quantitatively depict the fracture properties.
The technical scheme adopted by the invention is as follows: a method for quantitatively characterizing fracture mechanical properties based on structural style analysis comprises the following steps:
1) collecting a construction style: explaining the three-dimensional characteristics of the fracture based on the three-dimensional seismic data, the coherent data volume and the curvature data volume data, and reading the structural style characteristics related to the fracture;
2) preprocessing the construction style: based on the difference characteristics of fracture under the background of regional stress with different properties, selecting five types of construction styles capable of reflecting fracture mechanical properties, and establishing a detailed construction style database;
3) calculating a construction pattern, and solving the weight coefficient and the membership degree of the five types of construction patterns by establishing a membership function and a judgment matrix to preliminarily realize the conversion of geological features to quantitative data;
4) establishing a quantitative evaluation mathematical model based on the weight coefficient of the five types of construction styles and the membership degree of the single construction style obtained in the step 3), and obtaining the numerical value of the mechanical property of each breakpoint;
5) and drawing a plane layout diagram of the fracture mechanical property values of all the trunks by using Petrel software based on the solved mechanical property values of all the breakpoints, and completing the division of the quantitative evaluation standard of the fracture mechanical property of the land-phase fault-trap basin.
The method for quantitatively representing the fracture mechanical property based on the structural style analysis can efficiently, accurately and quantitatively depict the fracture property, thereby better serving research works such as fracture difference reservoir control action of different mechanical properties and the like. The invention has the beneficial effects that:
1. the method can be applied to quantitative depiction of mechanical properties of the developed fracture zone under the background of complex dynamics, and compared with the previous qualitative understanding of 'ambiguous' of the mechanical properties of the fracture zone, the method calculates the numerical values of the mechanical properties of the fracture zone by establishing a mathematical model and quantitatively divides the classification scheme of the mechanical properties of the fracture zone, so that the method can be more accurately and efficiently used for controlling oil-gas reservoir formation and related research work by the fracture difference of different properties;
2. the method not only can quantitatively depict the mechanical property of the target fracture zone integrally, but also can divide the mechanical property of the same fracture zone in a segmented manner based on the calculation of the point value of each fracture zone;
3. the mechanical properties of the fracture zone of the east continental facies fault-trap basin in China established by the method are unified and quantitatively evaluated, and after the mechanical properties of the fracture zone are quantitatively characterized, the comparative analysis of the mechanical properties of the fracture zone in different structural units can be realized, so that the effect of favorable target trap target optimization is achieved.
Drawings
FIG. 1 is a flow chart of a method for quantitatively characterizing fracture mechanical properties based on formation pattern analysis in accordance with the present invention;
FIG. 2 is a fine explanatory profile feature (ILN9800) of a break in the Bohai sea area, California 14-9 probe region;
FIG. 3 shows a Bohai sea area Jinzhou 14-9 exploration area plane fracture characteristic (T)3) And 5 backbone rupture zone divisions (F)1、F2、F3、F4、F5)；
FIG. 4 is a plot of the mechanical property values of 5 main fracture zones in the 14-9 probe region of the State based on Petrel software.
Detailed Description
The method for quantitatively characterizing fracture mechanical properties based on structural style analysis according to the present invention is described in detail below with reference to the following embodiments and the accompanying drawings.
As shown in FIG. 1, the method for quantitatively characterizing fracture mechanical properties based on structural style analysis of the invention comprises the following steps:
1) collecting a construction style: explaining the three-dimensional characteristics of the fracture based on the three-dimensional seismic data, the coherent data volume and the curvature data volume data, and reading the structural style characteristics related to the fracture; the method comprises the following steps: on the basis of a Petrel software platform, according to the fault relation of the homophase axes of three-dimensional seismic data in a research mining area, the fracture is explained, a coherent data body and a curvature data body are made by using a coherent algorithm and a curvature algorithm, coherent slices, curvature slices and fracture explanation results are mutually verified, the closure on a regional fracture explanation space is realized, fracture plane and fracture outline graphs of all layers of strata are realized, the trunk fracture of the research working area is combed according to fracture characteristics, and the construction style related to the fracture is read.
2) Preprocessing the construction style: based on the difference characteristics of fracture under the background of regional stress with different properties, selecting five types of construction styles capable of reflecting fracture mechanical properties, and establishing a detailed construction style database;
based on literature research, structural physical simulation experiments, numerical simulation and actual seismic interpretation schemes, the difference of fracture mechanical properties of fractures under different regional stress backgrounds shows that the five types of structural styles capable of reflecting the fracture mechanical properties are as follows: u1: fracture section combination pattern, U2: fracture plane spread morphology, U3: fault dip, U4: fracture surface cross-sectional morphology, U5: two disks of stratigraphic attitude. Wherein,
the fracture section combination mode comprises: a flower-like structure, a multi-stage Y-shaped structure and a Y-shaped structure; the fracture plane spreading form comprises: straight continuous, wild goose-row and saw-tooth configurations; the fault dip angle is divided into four ranges: 90 to 80 degrees, 80 to 60 degrees, 60 to 45 degrees and 45 degrees; the fracture surface profile comprises: plate, shovel, and ramp configurations; the two-disc stratum attitude comprises: crush-fold configuration, propensity-to-anomaly configuration, thickness-discontinuity configuration, yield-normal configuration.
The detailed construction pattern database is established by reading five types of construction patterns of each broken breakpoint at equal intervals based on a fracture system plan by taking the distance of 40 line numbers as intervals.
3) Calculating a construction pattern, namely solving the weight coefficient and the membership degree of the five types of construction patterns by establishing a membership function and a judgment matrix based on detailed statistics of the five types of construction patterns in the step 2), and preliminarily realizing the conversion of geological features to quantitative data; the method comprises the following steps:
(3.1) membership is Zadeh, L.A. the concept proposed in the "fuzzy sets" paper in 1965, which means that if there is a number F (x) e [0, 1] corresponding to any element x in the set U, F (x) is called membership of x to A. The method of the invention is based on the dynamic characteristics of the east China continental facies collapse basin, and defines: the walking and sliding effect is 1, the tension effect is 0, and on the basis of the standard, the membership value of each type of parameter is obtained by establishing a membership function;
the fracture section combination pattern is a combination characteristic of main branches and secondary fractures on a seismic section with vertical trunk fractures, and the flower-shaped structure, the flower-like structure, the multi-level Y-shaped structure and the Y-shaped structure represent that the sliding action of the fractures is weakened continuously and the tension action is strengthened continuously, so that the membership degrees of the flower-shaped structure, the flower-like structure, the multi-level Y-shaped structure and the Y-shaped structure are 1, 0.7, 0.5 and 0.3 respectively;
the fracture plane spreading form is a plane form which represents a main fracture section of a fracture system diagram of the same structural layer system, and the straight continuous structure, the wild goose-row structure and the sawtooth structure represent that the sliding action of fracture is continuously weakened, the tension action is continuously enhanced, and the sliding property fracture has small change from a tensile deep layer to a shallow layer; the horizontal fault-span of the tensile fault is greatly changed from the deep layer to the shallow layer, and the deep part is larger than the shallow part, so that the membership degrees of a straight continuous structure, an anser-type structure and a sawtooth-shaped structure are 1, 0.6 and 0.1 respectively;
the fault dip angle refers to an included angle between a fracture surface on a seismic section perpendicular to fracture and a horizontal plane, and the fracture dip angle range developed under stress environments with different properties has larger difference. The stratum in the Bohai sea area is broken under a stretching stress environment to form a fault, the dip angle of the positive fault is 60 degrees (45 degrees + theta/2) (theta is an internal friction angle of a rock, wherein the internal friction angle of the rock in the Bohai sea area is 30 degrees), and according to the specific breaking development characteristics of the Bohai sea area, the dip angle of the fault is divided into four ranges: 90-80 degrees, 80-60 degrees, 60-45 degrees and <45 degrees, and the slip property of the fracture is gradually enhanced along with the increase of the fault dip angle. Therefore, membership degrees of 90 degrees to >80 degrees, 80 degrees to >60 degrees, 60 degrees to >45 degrees and <45 degrees are 1, 0.8, 0.3 and 0.1, respectively;
the fracture section morphology is fracture track morphology on the seismic section perpendicular to the fracture trend of the trunk, and is divided into three categories according to the development characteristics of new generations of specific fractures of the fracture basin in the east of China: plate type, shovel type, sloping terrace type. The fracture section form is transited from a plate type to a sloping plateau type, which shows that the tensile property of the fracture is gradually enhanced, and the sliding property is gradually weakened. Thus, membership degrees of the plate type structure, the shovel type structure and the sloping plateau type structure are 1, 0.5 and 0.1 respectively;
the two-tray stratum attitude is a relative relation of tendency, thickness and deformation characteristics of two broken trays of stratums, according to the influence of different properties (sliding and stretching) of the Bohai sea area on the two broken trays of stratums of fracture development, the two-tray stratum attitude of the extended normal fault is taken as a standard, the two-tray stratum attitude of fracture caused by sliding is abnormal in a special displacement mode, wherein the abnormal stratum attitude caused by the sliding action in the Bohai sea area comprises: abnormal tendency, fold bending, abrupt thickness change. Therefore, the membership of the crush crease structure, the tendency abnormal structure, the thickness abrupt change structure and the occurrence normal structure is 0.9, 0.7, 0.5 and 0.2 respectively.
And (3.2) the weight coefficient is used for representing the importance of a certain index in the index system, and the size of the weight coefficient represents the size of the relevance of the index and the target importance degree. B is set as a fracture mechanical property quantitative characterization set; setting U as fuzzy subset of B to represent five types of construction style sets; setting V as a fuzzy subset of U to represent each single construction pattern factor set; firstly, establishing mapping from B to a function f (U) by comparing the importance of five types of construction styles pairwise, namely establishing a first-level hierarchical structure model judgment matrix; secondly, establishing mapping from U to function f (V) by pairwise comparing the importance of each factor in each type of construction style, namely establishing a secondary hierarchical structure model judgment matrix; finally, calculating and analyzing a structural hierarchical model by using yaahp hierarchical analysis software, setting the sum of the weight vectors of the five types of construction styles to be 1, and solving a single parameter weight coefficient; wherein,
(3.2.1) judging a matrix of the primary hierarchical structure model:
in the formula, aijRelative ratio representing the importance of two types of construction styles, where a11 ═U1/U1，a12＝U1/U2…；
(3.2.2) a secondary level model judgment matrix, wherein:
(a) fracture section combination pattern U1:
V-V1, V2, V3, V4, V5 flower-like, multi-stage Y-shaped, Y-shaped structures
In the formula, v11 represents the relative ratio of the importance of two factors in the fracture section combination pattern, wherein v11 is v1/v1, and a12 is v1/v2 …;
(b) fracture plane spread pattern U2:
v1, V2, V3 straight continuous structure, yankee structure and zigzag structure
In the formula, v11 represents the relative ratio of the importance of two factors in the fracture plane spreading form, wherein v11 is v1/v1, and a12 is v1/v2 …
(c) Fracture inclination angle U3:
V-V1, V2, V3, V4-90 to 80 degrees, 80 to 60 degrees, 60 to 45 degrees and 45 degrees
In the formula, v11 represents the relative ratio of the importance of two factors in the fracture inclination angle, wherein v11 is v1/v1, and a12 is v1/v2 …
(d) Fracture section morphology U4:
V-V1, V2, V3 are plate, shovel, or ramp structures
In the formula, v11 represents the relative ratio of the importance of two factors in fracture section morphology, wherein v11 is v1/v1, and a12 is v1/v2 …
(e) Two-disc formation attitude U5: 0.9, 0.7, 0.5, 0.2
V-V1, V2, V3, V4-crush-fold structure, tendency abnormal structure, thickness abrupt structure, and production normal structure
Where v11 represents the relative ratio of the importance of two factors in two-disc formation, where v 11-v 1/v1 and a 12-v 1/v2 …
Finally, weight coefficient standards of five types of construction styles are solved and established: u1: the weight coefficient of the fracture profile combination pattern was 0.4152, U2: the weight coefficient of the fracture plane spreading form is 0.2251, U3: the weight coefficient of the fault dip angle is 0.1578, U4: the weight coefficient of the fracture surface profile is 0.1163, U5: the two-disc formation attitude was weighted 0.0856.
4) Establishing a quantitative evaluation mathematical model based on the weight coefficient of the five types of construction styles and the membership degree of the single construction style obtained in the step 3), and obtaining the numerical value of the mechanical property of each breakpoint;
the quantitative evaluation mathematical model is as follows:
B＝U·R
[b1,b2,…,bn]＝[U1,U2,…,U5]·[rij]5×n
wherein, B ═ B1,b2,…,bn]Representing the quantitative value of the mechanical property of each breaking point; u ═ U1,U2,…,U5]Weight coefficients representing five types of construction styles; r ═ Rij]5×nAnd represents the membership of the jth breakpoint to the ith parameter evaluation.
The model can be used for solving the mechanical property value of each breaking point, so that the geological knowledge is really converted into quantitative data.
5) Based on the solved mechanical property value of each breakpoint, drawing a plane layout of the mechanical property value of each trunk fracture by using Petrel software to complete division of a quantitative evaluation standard of the mechanical property of the land-phase fracture-trap fracture; the method comprises the following steps: based on the solved mechanical property value of each breakpoint, the geodetic coordinates of the main fracture breakpoints are read by using double fox software, the geodetic coordinates of the breakpoints and the fracture mechanical property values are led into Petrel, a plane layout of the fracture mechanical property values of all the main fractures is drawn, the fracture mechanical property is analyzed more visually and accurately, and quantitative evaluation standards of the fracture mechanical property of the land-phase fracture-subsidence basin are divided.
Examples are given below:
aiming at a Bohai sea area Jinzhou 14-9 exploring area as an example, the specific implementation method comprises the following steps:
(1) the area of three-dimensional seismic data of 14-9 exploration areas of Jinzhou is 920km2Based on a Petrel software platform, fracture is finely explained according to the fault relation of the three-dimensional seismic data in-phase axis, a coherent data body and a curvature data body are made by using a coherent algorithm and a curvature algorithm and mutually verified with the fracture explanation result, and finally the iso-to closure in a regional fracture explanation space is realized, and a six-set-layer system fracture plane fracture class diagram (T) is implemented and drawn2、T3 U、T3 M、T3、T5、T8) And 5 trunk fracture zones (F) were marked1、F2、F3、F4、F5). Wherein FIG. 2 shows only the fracture characteristics (ILN9800) of one seismic line and a set of system plots (T) of fracture planes for layer series3)；
(2) Based on T3The fracture plane system diagram reads five types of structure styles of each fracture zone point of five fracture zones at equal intervals by taking 40 line number distances as intervals: u1 (fracture section combination style), U2 (fracture plane spread form), U3 (fault dip angle), U4 (fracture plane section form), U5 (two-disc stratum occurrence form), and FIG. 3 is a 5-class structure style feature table of each point of five main fracture zones;
(3) based on the obtained weight coefficients of the five types of parameters and the standard of single parameter membership, the weight coefficients and the membership are assigned to the 5 types of geological features read by the 5 trunk fracture zones to realize the preliminary conversion of the geological parameters to quantitative data, and FIG. 4 is a table of the weight coefficients and the membership values of the 5 types of construction patterns of each point of the five trunk fracture zones;
(4) and (4) solving the weight coefficients and membership values of various structural styles based on the step (3), and solving the mechanical property numerical values of the fracture points of the five fracture zones according to the established mathematical model. Take a certain point as an example
(5) Based on the mechanical property values of the 5 trunk fracture zones obtained in the step (4), a quantitative evaluation standard scheme for the fracture mechanical property of the east-China continental facies fracture-sink basin in China is established, and the scheme is as follows:
wherein B represents the fracture zone mechanical property value.
The geodetic coordinates of the trunk fracture zone points (F1, F2, F3 … fn) of the calculated mechanical property values are read by using double fox software, the mechanical property values and the geodetic coordinates are led into Petrel after being standardized, color icons are calibrated according to the evaluation standard scheme, a plane development diagram of the mechanical property values of all trunk fracture zones is drawn, the mechanical properties of the same fracture zone are found to have obvious segmentation from the color change of 5 trunk fracture zones in the diagram, and F1And F2The mechanical properties can be divided into three sections and four sections from south to north.
The above embodiments are only for illustrating the technical principle and practical application of the present invention, and the implementation steps of the method can be modified, and all changes or improvements based on the technical scheme of the present invention are within the protection scope of the present invention.
Claims (5)
1. A method for quantitatively characterizing fracture mechanical properties based on structural style analysis is characterized by comprising the following steps:
1) collecting a construction style: explaining the three-dimensional characteristics of the fracture based on the three-dimensional seismic data, the coherent data volume and the curvature data volume data, and reading the structural style characteristics related to the fracture;
2) preprocessing the construction style: based on the difference characteristics of fracture under the background of regional stress with different properties, selecting five types of construction styles capable of reflecting fracture mechanical properties, and establishing a detailed construction style database; the five types of structural styles capable of reflecting fracture mechanical properties are as follows:
u1: a fracture profile assembly pattern comprising: a flower-like structure, a multi-level Y-shaped structure and a single-level Y-shaped structure;
u2: fracture plane spreading form, including: straight continuous, wild goose-row and saw-tooth configurations;
u3: the fracture dip angle is divided into four ranges: 90 to 80 degrees, 80 to 60 degrees, 60 to 45 degrees and 45 degrees;
u4: fracture surface profile morphology comprising: plate, shovel, and tarmac constructions;
u5: two-plate stratigraphic attitude, comprising: an extrusion fold structure, a tendency abnormal structure, a thickness mutation structure and a production shape normal structure;
3) calculating a construction pattern, and solving the weight coefficient and the membership degree of the five types of construction patterns by establishing a membership function and a judgment matrix to preliminarily realize the conversion of geological features to quantitative data; the method comprises the following steps:
(3.1) definition: the walking and sliding effect is 1, the tension effect is 0, and the membership value of each type of parameter is obtained by establishing a membership function;
the fracture section combination pattern is a combination characteristic of a main branch and a secondary fracture on a seismic section which is vertical to a trunk fracture, and the membership degrees of a flower-shaped structure, a flower-like structure, a multi-level Y-shaped structure and a single-level Y-shaped structure are 1, 0.7, 0.5 and 0.3 respectively;
the fracture plane spreading form is a plane form representing a main fracture section of a fracture system diagram of the same structural layer series, and the membership degrees of a straight continuous structure, a goose-row structure and a sawtooth structure are 1, 0.6 and 0.1 respectively;
the fracture dip angle refers to an included angle between a fracture surface and a horizontal plane on a seismic section vertical to fracture, and the membership degrees of 90-80 degrees, 80-60 degrees, 60-45 degrees and 45 degrees are 1, 0.8, 0.3 and 0.1 respectively;
the fracture section form is a fracture track form on a seismic section perpendicular to the fracture trend of the trunk, and the membership degrees of the plate type structure, the shovel type structure and the slope level type structure are 1, 0.5 and 0.1 respectively;
the two-disc stratum attitude is the relative relation of the tendency, the thickness and the deformation characteristics of the fractured two-disc stratum, and the membership degrees of the extrusion fold structure, the tendency abnormal structure, the thickness mutation structure and the attitude normal structure are respectively 0.9, 0.7, 0.5 and 0.2;
(3.2) setting B as a name breaking point mechanical property quantitative value; setting U as fuzzy subset of B, representing weight coefficients of five types of construction styles; setting V as a fuzzy subset of U to represent each single construction pattern factor set; firstly, establishing mapping from B to a function f (U) by comparing the importance of five types of construction styles pairwise, namely establishing a first-level hierarchical structure model judgment matrix; secondly, establishing mapping from U to function f (V) by pairwise comparing the importance of each factor in each type of construction style, namely establishing a secondary hierarchical structure model judgment matrix; finally, calculating and analyzing a structural hierarchical model by using yaahp hierarchical analysis software, setting the sum of the weight vectors of the five types of construction styles to be 1, and solving a single parameter weight coefficient; wherein,
(3.2.1) judging a matrix of the primary hierarchical structure model:
in the formula, ai0j0A relative ratio representing the importance of the two types of construction styles;
(3.2.2) a secondary level model evaluation matrix, wherein:
(a) fracture section combination pattern U1:
V1＝v11、v12、v13、v14-flower-like structure, multi-stage Y-shaped structure, single-stage Y-shaped structure
In the formula, V1i1j1A relative ratio representing the importance of the two factors in the fracture profile combination pattern;
(b) fracture plane spread pattern U2:
V2＝v21、v22、v2straight continuous, wild goose-row and zigzag 3
In the formula, V2i2j2A relative ratio representing the importance of the two factors in the fracture plane spreading form;
(c) fracture inclination angle U3:
V3＝v31、v32、v33、v3e4 ═ 90 DEG E>Instant-wave of 80 DEG, 80 DEG>60, 60 degrees>45 degrees and<45 degrees (deg.)
In the formula, V3i3j3A relative ratio representing the importance of the two factors in fracture dip;
(d) fracture section morphology U4:
V4＝v41、v42、v43-plate structure, shovel structure, and ramp structure
In the formula, V4i4j4Representative fracture sectionThe relative ratio of the importance of the two factors in the surface morphology;
(e) two-disc formation attitude U5:
V5＝v51、v52、v53、v54-crush-fold structure, tendency-to-abnormal structure, thickness-abrupt structure, and occurrence-normal structure
In the formula, V5i5j5A relative ratio representing the importance of the two factors in the two-disc formation attitude;
finally, weight coefficient standards of five types of construction styles are solved and established: u1: the weight coefficient of the fracture profile combination pattern was 0.4152, U2: the weight coefficient of the fracture plane spreading form is 0.2251, U3: the weight coefficient of the fracture dip is 0.1578, U4: the weight coefficient of the fracture surface profile is 0.1163, U5: the weight coefficient of the two-disc formation attitude is 0.0856;
4) establishing a quantitative evaluation mathematical model based on the weight coefficient of the five types of construction styles and the membership degree of the single construction style obtained in the step 3), and obtaining the numerical value of the mechanical property of each breakpoint;
5) and drawing a plane layout diagram of the fracture mechanical property values of all the trunks by using Petrel software based on the solved mechanical property values of all the breakpoints, and completing the division of the quantitative evaluation standard of the fracture mechanical property of the land-phase fault-trap basin.
2. The method for quantitatively characterizing fracture mechanical properties based on the analysis of construction patterns according to claim 1, wherein step 1) comprises: on the basis of a Petrel software platform, according to the fault relation of the homophase axes of three-dimensional seismic data in a research mining area, the fracture is explained, a coherent data body and a curvature data body are made by using a coherent algorithm and a curvature algorithm, coherent slices, curvature slices and fracture explanation results are mutually verified, the closure on a regional fracture explanation space is realized, fracture plane and fracture outline graphs of all layers of strata are realized, the trunk fracture of the research working area is combed according to fracture characteristics, and the construction style related to the fracture is read.
3. The method for analyzing, quantitatively characterizing fracture mechanical properties based on construction styles according to claim 1, wherein the establishing of the detailed construction style database in step 2) is based on a fracture system plan, and the detailed construction style database is established by reading five types of construction styles of each fracture at equal intervals by taking 40 line number distances as intervals.
4. The method for quantitatively characterizing fracture mechanical properties based on structural style analysis according to claim 1, wherein the quantitative evaluation mathematical model of step 4) is as follows:
B＝U·R
[b1,b2,…,bn]＝[U1,U2,…,U5]·[rαβ]5×n
wherein, B ═ B1,b2,…,bn]Representing the quantitative value of the mechanical property of each breaking point; u ═ U1,U2,…,U5]Weight coefficients representing five types of construction styles; r ═ Rαβ]5×nAnd represents the degree of membership of the beta-th breakpoint to the alpha-th parameter evaluation.
5. The method for quantitatively characterizing fracture mechanical properties based on the analysis of construction patterns according to claim 1, wherein step 5) comprises: based on the solved mechanical property value of each breakpoint, the geodetic coordinates of the main fracture breakpoints are read by using double fox software, the geodetic coordinates of the breakpoints and the fracture mechanical property values are led into Petrel, a plane layout of the fracture mechanical property values of all the main fractures is drawn, the fracture mechanical property is analyzed more visually and accurately, and quantitative evaluation standards of the fracture mechanical property of the land-phase fracture-subsidence basin are divided.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011197303.9A CN112305609B (en) | 2020-10-31 | 2020-10-31 | Method for quantitatively characterizing fracture mechanical property based on structural style analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011197303.9A CN112305609B (en) | 2020-10-31 | 2020-10-31 | Method for quantitatively characterizing fracture mechanical property based on structural style analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112305609A CN112305609A (en) | 2021-02-02 |
| CN112305609B true CN112305609B (en) | 2022-06-07 |
Family
ID=74333361
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011197303.9A Active CN112305609B (en) | 2020-10-31 | 2020-10-31 | Method for quantitatively characterizing fracture mechanical property based on structural style analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112305609B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016037168A1 (en) * | 2014-09-05 | 2016-03-10 | Schlumberger Canada Limited | Well survivability in multidimensional geomechanical space |
| CN108680952A (en) * | 2018-04-03 | 2018-10-19 | 刘敬寿 | A kind of strike-slip fault structural evolution analytic method |
| CN110275205A (en) * | 2018-03-15 | 2019-09-24 | 中国石油化工股份有限公司 | A method of determining basin underground Small and Medium Sized Strike-slip faulted active times |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102982244A (en) * | 2012-11-30 | 2013-03-20 | 北京交通大学 | Fuzzy comprehensive evaluation method of geological strength index |
| US11231517B2 (en) * | 2018-02-27 | 2022-01-25 | Sanvean Technologies Llc | Azimuthal measurement for geosteering |
-
2020
- 2020-10-31 CN CN202011197303.9A patent/CN112305609B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016037168A1 (en) * | 2014-09-05 | 2016-03-10 | Schlumberger Canada Limited | Well survivability in multidimensional geomechanical space |
| CN110275205A (en) * | 2018-03-15 | 2019-09-24 | 中国石油化工股份有限公司 | A method of determining basin underground Small and Medium Sized Strike-slip faulted active times |
| CN108680952A (en) * | 2018-04-03 | 2018-10-19 | 刘敬寿 | A kind of strike-slip fault structural evolution analytic method |
Non-Patent Citations (6)
| Title |
|---|
| Cenozoic tectonic evolution of Liaodong dome, Northeast Liaodong Bay,Bohai, offshore China, constraints from seismic stratigraphy, vitrinite reflectance and apatite fission track data;Cheng Yanjun et al.;《Tectonophysics》;20150818;第152-165页 * |
| Fracture mechanics and mechanical fault detection by artificial intelligence methods: A review;Sara Nasiri et al.;《Engineering Failure Analysis》;20170708;第270-293页 * |
| 基于模糊数学的湘中凹陷圈闭评价;朱伟等;《特种油气藏》;20141231;第21卷(第06期);第1-5页 * |
| 渤海辽北地区新生代断裂特征及其演变过程;张江涛等;《现代地质》;20201013;第1-9页 * |
| 用模糊综合评判方法研究断层封闭性—以高邮凹陷陈堡地区为例;杨勇等;《新疆石油地质》;20050228;第26卷(第1期);第102-104页 * |
| 辽西凹陷北洼新生代断裂特征及控藏作用;张江涛等;《中国海上油气》;20171031;第29卷(第5期);第39-47页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112305609A (en) | 2021-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108680952B (en) | A kind of strike-slip fault structural evolution analytic method | |
| CN109102180B (en) | Comprehensive parameter evaluation method for double-dessert evaluation of tight sandstone reservoir | |
| Tariq et al. | A new artificial intelligence based empirical correlation to predict sonic travel time | |
| CN104614766B (en) | Seismic horizon skeletonization | |
| CN112419500B (en) | Three-dimensional geological model modeling method | |
| Ansari et al. | Clustering analysis of the seismic catalog of Iran | |
| CN102918423A (en) | Method for seismic hydrocarbon system analysis | |
| CN110007343B (en) | Method for evaluating complexity of regional geological structure | |
| CN101663667A (en) | Utilize comprehensive strata log diagram to generate the method for reservoir model | |
| CN106529667A (en) | Logging facies identification and analysis method based on fuzzy depth learning in big data environment | |
| CN101266299A (en) | Method for forecasting oil gas utilizing earthquake data object constructional features | |
| CN115166853A (en) | Shale gas reservoir natural fracture model establishing method and device, electronic equipment and medium | |
| CN110412662A (en) | Method for prediction of reservoirs of thin interbeded based on seismic multi-attribute deep learning | |
| CN112305609B (en) | Method for quantitatively characterizing fracture mechanical property based on structural style analysis | |
| CN117830561A (en) | Gold ore prospecting method based on three-dimensional geological modeling | |
| Chatterjee et al. | Ore grade estimation of a limestone deposit in India using an artificial neural network | |
| Sarkheil et al. | The fracture network modeling in naturally fractured reservoirs using artificial neural network based on image loges and core measurements | |
| CN106291715A (en) | A kind of low-grade fault law of development Forecasting Methodology based on fracture Self-similarity Theory | |
| Wong et al. | A neural‐network approach to knowledge‐based well interpolation: a case study OF a fluvial sandstone reservoir | |
| Wang et al. | Complex 3D geological modeling based on digital twin | |
| Li et al. | Slope deformation partitioning and monitoring points optimization based on cluster analysis | |
| Zellou et al. | Integrated fractured reservoir characterization using neural networks and fuzzy logic: Three case studies | |
| Li et al. | Construction and application of prediction methods for coal texture of CBM reservoirs at the block scale | |
| CN111414698A (en) | Corresponding analysis method for surface subsidence cause of subway tunnel excavation | |
| CN113240357A (en) | Rapid evaluation method for stability of highway side slope in red layer area |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |