CN113589372B - Quantitative identification method applicable to broken basin torsion Zhang Duanceng - Google Patents
Quantitative identification method applicable to broken basin torsion Zhang Duanceng Download PDFInfo
- Publication number
- CN113589372B CN113589372B CN202010370976.3A CN202010370976A CN113589372B CN 113589372 B CN113589372 B CN 113589372B CN 202010370976 A CN202010370976 A CN 202010370976A CN 113589372 B CN113589372 B CN 113589372B
- Authority
- CN
- China
- Prior art keywords
- fault
- calculating
- distance
- torsion
- data
- 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 25
- 238000006073 displacement reaction Methods 0.000 claims abstract description 21
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000011160 research Methods 0.000 claims description 3
- 230000001427 coherent effect Effects 0.000 claims description 2
- 238000012512 characterization method Methods 0.000 abstract description 7
- 208000010392 Bone Fractures Diseases 0.000 description 10
- 238000010276 construction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 238000007499 fusion processing Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009189 diving Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000001228 spectrum Methods 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
- 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/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
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
The invention provides a quantitative identification method suitable for torsion and tension faults in a land phase fracture basin, which comprises the following steps: extracting a two-dimensional seismic section based on the three-dimensional seismic data, the well data and the seismic coherence data volume; calculating a section true dip angle alpha based on a two-dimensional seismic section perpendicular to the fault strike; determining the horizontal movement distance L of a fault based on a fault cutting relation method, and calculating the vertical movement distance H of the fault and the displacement distance D along the fault tendency; calculating the total displacement distance M of the fault; calculating the proportion R of the horizontal fault distance L in the total fault displacement M, wherein R is defined as the slip ratio; and quantitatively identifying faults according to the magnitude of the slip ratio R. The quantitative identification method suitable for torsion and tension faults in the land fracture basin is simple and convenient to apply, comprehensive in consideration factors, reasonable in method, and wide in popularization and application value in exploration of complex structural bands of the torsion and tension basin, and quantitative characterization is realized.
Description
Technical Field
The invention relates to the technical field of oil field exploration and development, in particular to a quantitative identification method suitable for a fractured basin torsion Zhang Duanceng.
Background
In the context of extended formations, stick-slip effects are becoming increasingly important in oil field formation analysis. Under the influence of the diving of the Pacific plate and the Indian plate, more twisting factors are injected into the oil-gas-containing basin in China, and a series of twisting faults are formed by overlapping the stretching-collapsing actions of the basin. The study of the structure oil and gas reservoir controlled by the discovered torsion fault is carried out by students at home and abroad, and the result shows that the torsion fault influences the aspects of the material source channel, the deposition condition, the sedimentation rate, the trap type and the like of the oil and gas basin, has obvious control effect on the oil and gas conditions such as oil and gas generation, accumulation, migration, trap, preservation and the like, and is a favorable zone for the migration and accumulation of the oil and gas. However, the current identification of torsion and tension faults has some difficulties, and qualitative analysis can be performed only through the combined forms of some faults, and the fault is still in a blank state in the aspect of quantitative identification.
In application number: 201710247321.5, a fault identification method and device, the method comprises: acquiring seismic data of a target interval; performing spectrum decomposition processing on the seismic data based on a preset frequency interval and the frequency range of the seismic data to obtain a plurality of seismic frequency division data of the target interval; preprocessing the plurality of seismic frequency division data to obtain preprocessed seismic frequency division data; performing color fusion processing on the seismic data of the target interval based on the preprocessed seismic frequency division data; and identifying the position of the fault in the target interval based on the seismic data after the color fusion processing. The application is mainly based on frequency division data of seismic data, and is used for identifying faults, wherein the importance is that whether faults exist or not is identified, and the nature of the faults cannot be judged.
In application number: 201811128028.8, a method and a device for identifying faults and structural anomalies, the method comprises the following steps: acquiring a propagation rule of diffracted waves in a spherical coordinate system; acquiring three-dimensional forward seismic data; obtaining diffracted waves in the forward seismic data; finding an energy maximum point of the diffracted wave, and positioning the energy maximum point as a diffraction point of the diffracted wave; imaging the diffraction point, and identifying faults and structural abnormal bodies. The application mainly relates to the field of seismic data imaging, and by the method, the seismic data can reflect the fault morphology more clearly. The patent is also directed to imaging and identifying faults, and cannot be used for identifying different fault types such as walk, stretch, reverse fault and the like.
Therefore, the invention provides a novel quantitative characterization method suitable for identifying fault properties in a torsion structure, and solves the technical problems.
Disclosure of Invention
The invention aims to provide a quantitative characterization method which is simple and convenient to apply, considers the horizontal movement influence of faults and the vertical movement of the faults, comprehensively considers the activities of the faults in different directions and is suitable for the fault properties in a torsion and tension structure.
The aim of the invention can be achieved by the following technical measures: the quantitative characterization method suitable for the fault properties of the torsion structure comprises the following steps: step 1, extracting a two-dimensional seismic section based on three-dimensional seismic data, well data and a seismic coherent data volume; step 2, calculating a section true dip angle alpha based on a two-dimensional seismic section perpendicular to the fault strike; step 3, determining the horizontal movement distance L of the fault based on a fault cutting relation method, and calculating the vertical movement distance H of the fault and the displacement distance D along the fault tendency; step 4, calculating the total displacement distance M of the fault; step 5, calculating the proportion R of the horizontal fault distance L in the total fault displacement M, and defining R as the slip ratio; and 6, quantitatively identifying faults according to the magnitude of the sliding ratio R.
The aim of the invention can be achieved by the following technical measures:
in the step 1, three-dimensional seismic data, speed data and well data of a research area are comprehensively analyzed, fault coherence attribute analysis is carried out on the three-dimensional seismic data, and objective layer structure interpretation is carried out according to coherence data to obtain an objective layer fracture distribution map; and selecting a fracture zone to be analyzed, and extracting a two-dimensional seismic section perpendicular to the trend direction of the fracture zone.
In step 2, the two-dimensional seismic profile is converted into a two-dimensional geological profile according to the velocity data, and the true dip angle alpha of the profile is calculated.
In the step 3, determining the horizontal movement distance L of the fault along the fault trend based on a fault cutting relation method; along fault tendency, calculating the vertical movement distance H of the fault by using the drop of two discs of the fault, and calculating the movement distance D of the fault along tendency by using the vertical movement distance H and combining the inclination angle alpha of the fault, wherein D=H/Sin (alpha).
In step 4, the section inclination angle alpha is comprehensively considered, and the total displacement distance M of the fault is calculated based on the fault horizontal breaking distance L and the inclination displacement D, wherein
In step 5, the ratio R of the fault horizontal distance L in the fault total displacement M is calculated, R being defined as the slip ratio, where r=l/M.
In step 6, quantitatively identifying faults according to the magnitude of the slip ratio R; wherein when
When < 1, the fracture is a torsion fracture.
The quantitative characterization method suitable for the fault properties in the torsion structure is simple and convenient to apply, comprehensively considers the moving distance of the fault in the horizontal direction and the vertical direction, and establishes quantitative identification parameters of the torsion fault based on the proportion of the moving distance of the fault in the horizontal direction. The method has comprehensive consideration factors, realizes quantitative characterization, is reasonable, has wider popularization and application value in exploration of complex construction zones of the torsion basin, and provides favorable support for complex construction trap evaluation and oil and gas exploration.
Drawings
FIG. 1 is a flow chart of a quantitative identification method for a torsion fault according to an embodiment of the present invention;
FIG. 2 is a schematic diagram showing a method for calculating an interruption tilt angle α, a vertical interruption distance H and a tendency interruption distance D according to an embodiment of the present invention;
FIG. 3 is a schematic diagram illustrating a method for calculating an interruption layer trend displacement L according to an embodiment of the present invention;
FIG. 4 is a schematic diagram showing the calculation of the total displacement M of the interrupt layer and the quantitative identification parameters of the torsion fault according to an embodiment of the present invention;
fig. 5 is a plan view of a torsion Zhang Duanceng in accordance with an embodiment of the present invention.
Detailed Description
It should be noted that the following detailed description is exemplary and is intended to provide further explanation of the invention. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the present invention. As used herein, the singular forms also are intended to include the plural forms unless the context clearly indicates otherwise, and furthermore, it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, and/or combinations thereof.
As shown in fig. 1, fig. 1 is a flowchart of a quantitative identification method of a torsion fault according to the present invention.
In step 101, three-dimensional seismic data, velocity data and well data of a research area are comprehensively analyzed, fault coherence attribute analysis is conducted on the three-dimensional seismic data, and target layer structure interpretation is conducted according to the coherence data, so that a target layer fracture distribution map is obtained. Selecting a fracture zone to be analyzed, extracting a two-dimensional seismic section perpendicular to the direction of the fracture zone, as shown in fig. 2, and the flow proceeds to step 102;
in step 102, converting the two-dimensional seismic section into a two-dimensional geological section according to the velocity data, and calculating the true inclination angle α of the section, as shown in fig. 2, the flow proceeds to step 103;
in step 103, determining a horizontal movement distance L of the fault along the fault trend based on a fault cutting relationship method, as shown in fig. 3; along the fault trend, calculating the vertical movement distance H of the fault by using the drop of two discs of the fault, and calculating the movement distance D of the fault along the trend by using the vertical movement distance H and combining the inclination angle alpha of the fault, wherein D=H/Sin (alpha), as shown in fig. 2. The flow proceeds to step 104;
in step 104, the section inclination angle alpha is comprehensively considered and is based onThe total displacement distance M of the fault is calculated by the horizontal fault distance L and the inclined displacement D, wherein
As in fig. 4. The flow proceeds to step 105;
in step 105, calculating a ratio R of the fault horizontal distance L in the fault total displacement M, defining R as a walk-slip ratio, where r=l/M, and the flow proceeds to step 106;
in step 106, faults are quantitatively identified according to the magnitude of the slip ratio R. Wherein when
When broken, the fracture is a torsion fracture, as shown in fig. 5. The flow ends.
The quantitative identification method for the torsion and tension faults is simple and convenient to apply, comprehensively considers the moving distance of the faults in the horizontal direction and the vertical direction, and establishes quantitative identification parameters of the torsion and tension faults based on the proportion of the total moving distance of the faults in the horizontal direction. The method has comprehensive consideration factors, realizes quantitative characterization, is reasonable, has wider popularization and application value in exploration of complex construction zones of the torsion basin, and provides favorable support for complex construction trap evaluation and oil and gas exploration.
Finally, it should be noted that: the foregoing description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, but although the present invention has been described in detail with reference to the foregoing embodiment, it will be apparent to those skilled in the art that modifications may be made to the technical solution described in the foregoing embodiment, or equivalents may be substituted for some of the technical features thereof. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Other than the technical features described in the specification, all are known to those skilled in the art.
Claims (1)
1. The quantitative identification method for the land fracture basin torsion Zhang Duanceng is characterized by comprising the following steps of:
step 1, extracting a two-dimensional seismic section based on three-dimensional seismic data, well data and a seismic coherent data volume;
step 2, calculating a section true dip angle alpha based on a two-dimensional seismic section perpendicular to the fault strike;
step 3, determining the horizontal movement distance L of the fault based on a fault cutting relation method, and calculating the vertical movement distance H of the fault and the displacement distance D along the fault tendency;
step 4, calculating the total displacement distance M of the fault;
step 5, calculating the proportion R of the horizontal movement distance L of the fault in the total displacement distance M of the fault, and defining R as the slip ratio;
step 6, quantitatively identifying faults according to the magnitude of the sliding ratio R;
in the step 1, three-dimensional seismic data, speed data and well data of a research area are comprehensively analyzed, fault coherence attribute analysis is carried out on the three-dimensional seismic data, and objective layer structure interpretation is carried out according to coherence data to obtain an objective layer fracture distribution map; selecting a fracture zone to be analyzed, and extracting a two-dimensional seismic section perpendicular to the trend direction of the fracture zone;
in step 2, converting the two-dimensional seismic section into a two-dimensional geological section according to the speed data, and calculating the true inclination angle alpha of the section;
in the step 3, determining the horizontal movement distance L of the fault along the fault trend based on a fault cutting relation method; calculating the vertical movement distance H of the fault by using the drop of two discs of the fault along the fault tendency, and calculating the displacement distance D of the fault along the tendency by using the vertical movement distance H and combining the true inclination angle alpha of the fault, wherein D=H/Sin (alpha);
in step 4, the true inclination angle alpha of the fracture surface is comprehensively considered, and the total displacement distance M of the fracture surface is calculated based on the horizontal movement distance L of the fracture surface and the displacement distance D of the fracture surface along the trend, wherein
In step 5, calculating a ratio R of the horizontal movement distance L of the fault to the total displacement distance M of the fault, defining R as a slip ratio, wherein r=l/M;
in step 6, quantitatively identifying faults according to the magnitude of the slip ratio R; wherein when
When the fault is a torsion fault.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010370976.3A CN113589372B (en) | 2020-04-30 | 2020-04-30 | Quantitative identification method applicable to broken basin torsion Zhang Duanceng |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010370976.3A CN113589372B (en) | 2020-04-30 | 2020-04-30 | Quantitative identification method applicable to broken basin torsion Zhang Duanceng |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113589372A CN113589372A (en) | 2021-11-02 |
| CN113589372B true CN113589372B (en) | 2023-06-30 |
Family
ID=78238033
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010370976.3A Active CN113589372B (en) | 2020-04-30 | 2020-04-30 | Quantitative identification method applicable to broken basin torsion Zhang Duanceng |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113589372B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110618454A (en) * | 2019-10-24 | 2019-12-27 | 西南石油大学 | Identification method for glide fracture development direction in sedimentary basin |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106569268A (en) * | 2015-10-10 | 2017-04-19 | 中国石油化工股份有限公司 | Method for quantitatively identifying type of fault structure |
| FR3043481A1 (en) * | 2015-11-05 | 2017-05-12 | Services Petroliers Schlumberger | |
| CN106199709B (en) * | 2016-06-29 | 2018-02-13 | 中国石油化工股份有限公司 | Turn round the evaluation method of tensile crack activity |
-
2020
- 2020-04-30 CN CN202010370976.3A patent/CN113589372B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110618454A (en) * | 2019-10-24 | 2019-12-27 | 西南石油大学 | Identification method for glide fracture development direction in sedimentary basin |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113589372A (en) | 2021-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109828316B (en) | Calcium-bonded rock type uranium ore prospecting method | |
| CN102200588B (en) | Method for analyzing waveform similarity body curvature of seismic data | |
| CN104747183A (en) | Carbonate reservoir comprehensive classification method | |
| CN110618454B (en) | Identification method for glide fracture development direction in sedimentary basin | |
| WO2023000257A1 (en) | Geological-seismic three-dimensional prediction method for favorable metallogenic site of sandstone-type uranium deposit | |
| CN105158802B (en) | Lacustrine Basins Gravity Flow Sediments well logging quantitative identification method | |
| CN106569268A (en) | Method for quantitatively identifying type of fault structure | |
| Krystopowicz et al. | Tectonic geomorphology and Plio-Quaternary structural evolution of the Tuzgölü fault zone, Turkey: Implications for deformation in the interior of the Central Anatolian Plateau | |
| CN111812707A (en) | Method for detecting boulders in subway line site based on micro-motion and surface wave exploration | |
| CN111679340A (en) | Ore seeking method for black rock series vanadium-molybdenum ore | |
| CN113589372B (en) | Quantitative identification method applicable to broken basin torsion Zhang Duanceng | |
| CN109469477B (en) | Method and device for predicting extension direction of artificial crack | |
| CN106199690A (en) | The Forecasting Methodology in mud shale crack | |
| CN109324171A (en) | A kind of sedimentary facies quantitative identification method based on lithology statistics | |
| Nardone et al. | Shallow velocity structure and site effects at Mt. Vesuvius, Italy, from HVSR and array measurements of ambient vibrations | |
| CN112433253A (en) | Three-dimensional detection method for surrounding rock tunnel with weak broken flowing water | |
| CN102508311B (en) | Tunnel advanced detection data multi-parameter spatial mapping method | |
| Benedicto et al. | Advanced use of borehole acoustic televiewer (ATV) for structural interpretation of unconformity-related uranium deposits | |
| CN114910953B (en) | Method for explaining irregular section of complex fracture zone | |
| Fedorova et al. | Analysis of variance amplitudes of signals for detecting structural permafrost heterogeneities by ground penetrating radar | |
| RU2431869C2 (en) | Seismic prospecting method | |
| Harris et al. | Subseismic Fault Identification Using the Fault Likelihood Attribute: Application to Geosteering in the DJ Basin | |
| CN105676278A (en) | Method and device for processing seismic wave data | |
| Sokolov et al. | GPR investigation of ice-filled cracks in loose deposits | |
| Mei et al. | Research on high‐precision gravity scratch analysis and structure features of Yili Basin, north‐west China |
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 |