CN103941287A - Rapid three-dimensional fault interpretation method based on horizontal navigation - Google Patents
Rapid three-dimensional fault interpretation method based on horizontal navigation Download PDFInfo
- Publication number
- CN103941287A CN103941287A CN201410152774.6A CN201410152774A CN103941287A CN 103941287 A CN103941287 A CN 103941287A CN 201410152774 A CN201410152774 A CN 201410152774A CN 103941287 A CN103941287 A CN 103941287A
- Authority
- CN
- China
- Prior art keywords
- fault
- interpretation
- section
- seismic
- tomography
- 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a rapid three-dimensional fault interpretation method based on horizontal navigation. The method includes the steps that (1), the fault enhancement interpretative processing is carried out on seismic data, and a seismic data volume fully highlighting fault information is obtained; (2), time slices of which the interest interval reflects the fault plane distribution characteristic are extracted from the seismic data volume, and the time slices serve as horizontal navigation slices for fault interpretation; (3), a fault is selected from the horizontal navigation slices; (4), a seismic section is extracted in the direction perpendicular to the fault strike, and the section of the fault is interpreted at the projection point of the section of the fault; (5), meshing processing is carried out on a navigation plane interpretation result of the selected fault and an interpretation result of the fault of the extracted seismic section, and the results are projected to a newly extracted seismic section, so that the fault interpretation of the seismic section is effectively navigated; (6), under the navigation of the fault plane interpretation, seismic sections are extracted in a rolling mode in the direction perpendicular to the fault strike step by step, and the fault interpretation of the seismic sections is finished until all the seismic sections of the fault is finished; (7), the steps (3)-(6) are carried out repeatedly until interpretation of all faults of the whole area is finished.
Description
Technical field
The present invention relates to a kind of fault interpretation method in petroleum exploration domain, particularly about a kind of quick three-dimensional fault interpretation method based on horizontal navigation.
Background technology
Seismic data interpretation comprises that explain fault interpretation, layer position and target geologic body is explained three piths, and what still adopt for the fault interpretation in seismic data interpretation at present is two-dimensional interpretation, and the three dimensions of not realizing is truly explained.Mainly there is following shortcoming in existing fault interpretation technology: 1, not closed-ended question of fault interpretation Existential Space, and error is larger; 2,3-D data volume vacuate two dimension fault interpretation, information utilizes degree low; 3, the fault interpretation based on two-dimentional, plane information application is not enough, and target geologic body understanding is unclear.
Summary of the invention
For the problems referred to above, the object of this invention is to provide a kind of can guarantee and explain closure, in fault interpretation, complete the plane coupling combination of tomography and explain the high quick three-dimensional fault interpretation method based on horizontal navigation of efficiency one by one at seismic section rolling process interrupting layer.
For achieving the above object, the present invention takes following technical scheme: a kind of quick three-dimensional fault interpretation method based on horizontal navigation, comprise the following steps: 1) to the geological data application variance for fault interpretation, curvature, the tomography edge detection method such as relevant, carry out the explanatory processing of tomography reinforcement, so that the breakpoint of tomography and section are more clear on geological data section, obtain fully giving prominence to the seismic data volume of fault information; 2) on seismic data volume, extract the dropping cut slice of objective interval different time, extract the time slice that can reflect fault plane Distribution Characteristics, the horizontal navigation section using this time slice as fault interpretation; 3) in horizontal navigation section, choose a tomography according to the descending principle of fault plane extended length; 4) in the direction perpendicular to selected fault strike, extract seismic section, complete the Interpretation of profile of this tomography at the subpoint place of selected tomography section; 5) on the seismic section of the navigation plane of selected tomography being explained and being extracted, the interpretation results of tomography is carried out gridding processing according to fault interpretation density, fault plane after treatment gridding is projected to new extraction on another seismic section to be explained, on another seismic section, obtain tomographic projection line clearly, with the fault interpretation of the seismic section that effectively navigates; 6) under the navigation directions of explaining at fault plane, according to 3-D seismics fault interpretation density, and in the local refinement of tomography complex region, along the fault interpretation that progressively rolls perpendicular to fault strike direction and extract seismic section and complete this seismic section, until complete the explanation of all seismic sections of this tomography; 7) repeating said steps 3)～6), until complete all fault interpretations of the whole district.
The present invention is owing to taking above technical scheme, and it has the following advantages: 1, because the present invention takes full advantage of horizontal navigation information with fault interpretation one by one, and in fault interpretation one by one, completed the plane coupling combination of tomography; Also utilize in addition horizontal navigation information accurately to locate the fault interpretation of section, therefore, ensured to explain it is closed at seismic section rolling process interrupting layer.2, because the present invention only needs to explain the seismic section perpendicular to fault strike direction, and do not need to explain that seismic section that parallel fault moves towards direction just can ensure the space closure of tomography, therefore greatly save the time of fault interpretation, improved the efficiency of fault interpretation.
Brief description of the drawings
Fig. 1 is that time slice extracts schematic diagram;
Fig. 2 is the horizontal navigation section schematic diagram of certain oil field fault interpretation;
Fig. 3 is that the plane of F1 tomography is explained schematic diagram;
Fig. 4 is perpendicular to F1 fault strike direction and chooses seismic section planimetric position schematic diagram 1., 2., 3.
Fig. 5 is the horizontal navigation Interpretation of profile schematic diagram of F1 tomography, and in figure, a represents navigation section extracting position, the b slice plane fault interpretation subpoint that represents to navigate, and c represents the tomography grid surface of horizontal navigation;
Fig. 6 is the plane projection schematic diagram of F1 Fault profile interpretation results;
Fig. 7 is the achievement diagrammatic cross-section of certain oil field F1 fault interpretation;
Fig. 8 is another achievement diagrammatic cross-section of certain oil field F1 fault interpretation.
Embodiment
Below in conjunction with drawings and Examples, the present invention is described in detail.
The invention provides a kind of quick three-dimensional fault interpretation method based on horizontal navigation, comprise the following steps:
1) to the geological data application variance for fault interpretation, curvature, the tomography edge detection method such as relevant, carry out the explanatory processing of tomography reinforcement, so that the breakpoint of tomography and section are more clear on geological data section, obtain fully giving prominence to the seismic data volume of fault information.
2) on seismic data volume, extract the dropping cut slice of objective interval different time, extract the time slice that can reflect fault plane Distribution Characteristics, the horizontal navigation section using this time slice as fault interpretation.
3) in horizontal navigation section, choose a tomography according to the descending principle of fault plane extended length.
4) in the direction perpendicular to selected fault strike, extract seismic section, complete the Interpretation of profile of this tomography at the subpoint place of selected tomography section.
5) on the seismic section of the navigation plane of selected tomography being explained and being extracted, the interpretation results of tomography is carried out gridding processing according to fault interpretation density, fault plane after treatment gridding is projected to new extraction on another seismic section to be explained, on another seismic section, obtain tomographic projection line clearly, with the fault interpretation of the seismic section that effectively navigates.
6) under the navigation directions of explaining at fault plane, according to 3-D seismics fault interpretation density, and in the local refinement of tomography complex region, along the fault interpretation that progressively rolls perpendicular to fault strike direction and extract seismic section and complete this seismic section, until complete the explanation of all seismic sections of this tomography.
7) repeating step 3)～6), until complete all fault interpretations of the whole district.
Enumerate specific embodiment below, to have a better understanding of the present invention.
Embodiment: the tomography in certain Complicated Fault Block Oilfield district is the main thoroughfare of oil-gas migration, also be one of Dominated Factors of Hydrocarbon Formation Reservoirs, for more clearly showing the faults distribution feature in this Complicated Fault Block Oilfield district, complete fast the fault interpretation in this Complicated Fault Block Oilfield district, adopt the present invention to carry out fault interpretation, comprise the following steps:
1) to the geological data application variance for fault interpretation, curvature, the tomography edge detection method such as relevant, carry out the explanatory processing of tomography reinforcement, so that the breakpoint of tomography and section are more clear on geological data section, obtain fully giving prominence to the seismic data volume of fault information.
2) as shown in Figure 1, on seismic data volume, extract the dropping cut slice of objective interval different time t1, t2 and t3, extract the time slice 1,2 and 3 that can reflect fault plane Distribution Characteristics.As shown in Figure 2, in figure, white portion is the planar distribution form of tomography, the horizontal navigation section using the time slice extracting as fault interpretation.
3) in horizontal navigation section, choose a tomography F1(as shown in Figure 3 according to the descending principle of fault plane extended length).
4) move towards perpendicular to tomography F1 to extract in direction seismic section 1., 2. (as shown in Figure 4), tomography section completes the Interpretation of profile (as shown in Figure 5) of tomography on sympathetic earthquakes section at seismic section subpoint b place 1., 2..
5) navigation plane of F1 tomography is explained and seismic section 1. with seismic section 2. interpretation results carry out gridding processing according to fault interpretation density, 3. the seismic section to be explained that fault plane after treatment gridding is projected in to new extraction is gone up, 3. obtain a tomographic projection line clearly at seismic section on, the effectively fault interpretation of navigation seismic section, as shown in the dotted line c of the seismic section in Fig. 5 in 3..
6) under the navigation directions of explaining at fault plane, according to 3-D seismics fault interpretation density, and in the local refinement of tomography complex region, along the fault interpretation that progressively rolls perpendicular to fault strike direction and extract seismic section and complete this seismic section, until complete the explanation of all seismic sections of F1 tomography, the plane projection of F1 Fault profile interpretation results as shown in Figure 6.
7) repeating step 3)～6), until complete all fault interpretations of the whole district, fault interpretation achievement is as shown in Figure 7, Figure 8.
The various embodiments described above are only for illustrating the present invention; wherein the structure of each parts, connected mode and manufacture craft etc. all can change to some extent; every equivalents of carrying out on the basis of technical solution of the present invention and improvement, all should not get rid of outside protection scope of the present invention.
Claims (1)
1. the quick three-dimensional fault interpretation method based on horizontal navigation, comprises the following steps:
1) to the geological data application variance for fault interpretation, curvature, the tomography edge detection method such as relevant, carry out the explanatory processing of tomography reinforcement, so that the breakpoint of tomography and section are more clear on geological data section, obtain fully giving prominence to the seismic data volume of fault information;
2) on seismic data volume, extract the dropping cut slice of objective interval different time, extract the time slice that can reflect fault plane Distribution Characteristics, the horizontal navigation section using this time slice as fault interpretation;
3) in horizontal navigation section, choose a tomography according to the descending principle of fault plane extended length;
4) in the direction perpendicular to selected fault strike, extract seismic section, complete the Interpretation of profile of this tomography at the subpoint place of selected tomography section;
5) on the seismic section of the navigation plane of selected tomography being explained and being extracted, the interpretation results of tomography is carried out gridding processing according to fault interpretation density, fault plane after treatment gridding is projected to new extraction on another seismic section to be explained, on another seismic section, obtain tomographic projection line clearly, with the fault interpretation of the seismic section that effectively navigates;
6) under the navigation directions of explaining at fault plane, according to 3-D seismics fault interpretation density, and in the local refinement of tomography complex region, along the fault interpretation that progressively rolls perpendicular to fault strike direction and extract seismic section and complete this seismic section, until complete the explanation of all seismic sections of this tomography;
7) repeating said steps 3)～6), until complete all fault interpretations of the whole district.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410152774.6A CN103941287B (en) | 2014-04-16 | 2014-04-16 | A kind of quick three-dimensional fault interpretation method based on horizontal navigation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410152774.6A CN103941287B (en) | 2014-04-16 | 2014-04-16 | A kind of quick three-dimensional fault interpretation method based on horizontal navigation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103941287A true CN103941287A (en) | 2014-07-23 |
| CN103941287B CN103941287B (en) | 2016-10-12 |
Family
ID=51189032
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410152774.6A Active CN103941287B (en) | 2014-04-16 | 2014-04-16 | A kind of quick three-dimensional fault interpretation method based on horizontal navigation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103941287B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104570096A (en) * | 2015-01-22 | 2015-04-29 | 中国海洋石油总公司 | Granite buried hill top surface fault interpretation method |
| CN105911589A (en) * | 2016-06-27 | 2016-08-31 | 中国石油大学(华东) | Quantifying method for confirming fault sliding parameter in three-dimensional space by adopting seismic data |
| CN107870360A (en) * | 2016-09-28 | 2018-04-03 | 中国石油化工股份有限公司 | Earthquake sedimentary evolution animation producing method and device |
| CN108254784A (en) * | 2017-12-11 | 2018-07-06 | 中国石油天然气股份有限公司 | Fault identification method, device and system based on two-dimensional seismic data |
| CN110488353A (en) * | 2019-09-04 | 2019-11-22 | 中海石油(中国)有限公司 | It is a kind of based on the fault interpretation method cuing open flat interaction and being combined with tectonic style guidance |
| CN112505757A (en) * | 2020-11-17 | 2021-03-16 | 中国石油集团工程咨询有限责任公司 | Three-dimensional seismic fault data processing method based on parallel grid control model |
| CN113945973A (en) * | 2020-07-17 | 2022-01-18 | 中国石油化工股份有限公司 | Reservoir characteristic analysis method, storage medium and electronic equipment |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6018498A (en) * | 1998-09-02 | 2000-01-25 | Phillips Petroleum Company | Automated seismic fault detection and picking |
| CN1504763A (en) * | 2002-12-05 | 2004-06-16 | 大庆石油管理局 | Method for extracting coherent slice faultage polygon |
| US6791900B2 (en) * | 2002-06-13 | 2004-09-14 | Exxonmobil Upstream Research Company | Method of calculating a throw volume for quantitative fault analysis |
| CN103245971A (en) * | 2013-04-25 | 2013-08-14 | 中国石油天然气股份有限公司 | Well break point guided seismic minor fault interpretation method and device |
-
2014
- 2014-04-16 CN CN201410152774.6A patent/CN103941287B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6018498A (en) * | 1998-09-02 | 2000-01-25 | Phillips Petroleum Company | Automated seismic fault detection and picking |
| US6791900B2 (en) * | 2002-06-13 | 2004-09-14 | Exxonmobil Upstream Research Company | Method of calculating a throw volume for quantitative fault analysis |
| CN1504763A (en) * | 2002-12-05 | 2004-06-16 | 大庆石油管理局 | Method for extracting coherent slice faultage polygon |
| CN103245971A (en) * | 2013-04-25 | 2013-08-14 | 中国石油天然气股份有限公司 | Well break point guided seismic minor fault interpretation method and device |
Non-Patent Citations (4)
| Title |
|---|
| 伍国勇 等: "高邮凹陷北部斜坡瓦庄东地区断裂特征研究", 《石油天然气学报》, vol. 32, no. 4, 31 August 2010 (2010-08-31), pages 217 - 219 * |
| 岳周会: "孤东地区馆陶组河道砂体油藏勘探技术研究", 《中国矿业》, vol. 19, no. 2, 28 February 2010 (2010-02-28), pages 117 - 120 * |
| 朱庆荣 等: "复杂断层快速三维体解释方法", 《石油物探》, vol. 42, no. 4, 31 December 2003 (2003-12-31), pages 456 - 459 * |
| 李玲 等: "用地震相干数据体进行断层自动解释", 《石油地球物理勘探》, vol. 33, 30 June 1998 (1998-06-30), pages 105 - 111 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104570096A (en) * | 2015-01-22 | 2015-04-29 | 中国海洋石油总公司 | Granite buried hill top surface fault interpretation method |
| CN105911589A (en) * | 2016-06-27 | 2016-08-31 | 中国石油大学(华东) | Quantifying method for confirming fault sliding parameter in three-dimensional space by adopting seismic data |
| CN107870360A (en) * | 2016-09-28 | 2018-04-03 | 中国石油化工股份有限公司 | Earthquake sedimentary evolution animation producing method and device |
| CN108254784A (en) * | 2017-12-11 | 2018-07-06 | 中国石油天然气股份有限公司 | Fault identification method, device and system based on two-dimensional seismic data |
| CN108254784B (en) * | 2017-12-11 | 2020-09-08 | 中国石油天然气股份有限公司 | Fault identification method, device and system based on two-dimensional seismic data |
| CN110488353A (en) * | 2019-09-04 | 2019-11-22 | 中海石油(中国)有限公司 | It is a kind of based on the fault interpretation method cuing open flat interaction and being combined with tectonic style guidance |
| CN113945973A (en) * | 2020-07-17 | 2022-01-18 | 中国石油化工股份有限公司 | Reservoir characteristic analysis method, storage medium and electronic equipment |
| CN113945973B (en) * | 2020-07-17 | 2024-04-09 | 中国石油化工股份有限公司 | Reservoir characteristic analysis method, storage medium and electronic equipment |
| CN112505757A (en) * | 2020-11-17 | 2021-03-16 | 中国石油集团工程咨询有限责任公司 | Three-dimensional seismic fault data processing method based on parallel grid control model |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103941287B (en) | 2016-10-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103941287A (en) | Rapid three-dimensional fault interpretation method based on horizontal navigation | |
| CN104267431B (en) | The craven fault acquisition methods that a kind of seismotectonics is explained | |
| CN103941286B (en) | A kind of quick three-dimensional layer position based on horizontal navigation means of interpretation | |
| CN103941288A (en) | Three-dimensional target geologic body interpretation method based on horizontal navigation | |
| CN108415077A (en) | New edge detection low order fault recognition methods | |
| CN104166163A (en) | Method for automatically extracting fault curved surface based on three-dimensional large-data-volume seismic data cube | |
| CN103792576A (en) | Reservoir heterogeneous detection method and device based on gradient structure tensor | |
| CN105866832A (en) | Method and device for identifying low-order fault of deep tight sandstone gas reservoir | |
| CN104459768A (en) | Three-dimensional space target geologic body tracking method based on visualization | |
| CN104316958A (en) | Coherent processing method for identifying different scales of formation fractures | |
| HAO et al. | A study on the extension of fault zones in Yellow Sea and its adjacent areas based on gravity data | |
| O'Briain et al. | Improved subsalt imaging and salt interpretation by RTM scenario testing and image partitioning | |
| CN103777241B (en) | Three dimensional seismic data rapid edge-detection method based on time domain generalized Hilbert transformation | |
| CN105607129A (en) | Seismic identification method for complex fault block unconformity surfaces | |
| Taixian et al. | Techniques for high-efficient development of offshore fluvial oilfields | |
| CN104777513A (en) | Non-continuity boundary detection method for seismic data gradient information | |
| CN104199099B (en) | A kind of method utilizing multiwindow serial section efficiently to carry out seismic interpretation | |
| CN102914790A (en) | Observing system method for primary collection of two-dimensional observing system and three-dimensional observing system | |
| Zhao et al. | An automatical infill shot method for uniform imaging of target layer | |
| CN104714248A (en) | Formation overlap-denudation line recognition method based on iteration earthquake DNA detection | |
| Gao et al. | Curvature attribute based on dip scan with eccentric window | |
| CN112213777B (en) | Geosteering phase interpretation method and device | |
| Coimbra et al. | Migration velocity analysis with diffraction events using residual moveout: Application to SIGSBEE 2B data | |
| Yakymchuk et al. | Experience of Mobile Geophysical Methods Using for Hydrocarbon Potential Assessment of Some Structures in Barents Sea | |
| CN105629302A (en) | Three-dimensional seismic directional volumetric curvature analysis method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee after: CNOOC research institute limited liability company Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee before: CNOOC Research Institute Patentee before: China National Offshore Oil Corporation |
|
| CP01 | Change in the name or title of a patent holder |