CN110488353B - Fault interpretation method based on combination of profile interaction and construction style guidance - Google Patents
Fault interpretation method based on combination of profile interaction and construction style guidance Download PDFInfo
- Publication number
- CN110488353B CN110488353B CN201910831798.7A CN201910831798A CN110488353B CN 110488353 B CN110488353 B CN 110488353B CN 201910831798 A CN201910831798 A CN 201910831798A CN 110488353 B CN110488353 B CN 110488353B
- Authority
- CN
- China
- Prior art keywords
- fault
- section
- interpretation
- construction
- pattern
- 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
- 238000010276 construction Methods 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000003993 interaction Effects 0.000 title claims abstract description 25
- 238000005516 engineering process Methods 0.000 claims abstract description 26
- 238000012545 processing Methods 0.000 claims abstract description 9
- 238000003908 quality control method Methods 0.000 claims abstract description 9
- 238000012800 visualization Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 9
- 238000011161 development Methods 0.000 claims description 7
- 230000002159 abnormal effect Effects 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000005520 cutting process Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 3
- 230000001629 suppression Effects 0.000 claims description 3
- 238000002224 dissection Methods 0.000 claims 2
- 230000000694 effects Effects 0.000 description 4
- 244000258825 Scaevola taccada Species 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000007794 visualization technique Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 1
- 230000009189 diving Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000011179 visual inspection 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
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)
- Analysing Materials By The Use Of Radiation (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a fault interpretation method based on combination of profile interaction and construction style guidance, which comprises the following steps of: 1) performing explanatory processing on the seismic data for fault interpretation; 2) making an attribute slice for explaining the vicinity of a target layer by using a variance body and ant body technology; 3) analyzing a construction stress field and a section typical construction pattern of the work area, and establishing a construction pattern library; 4) carrying out fault interpretation on a section perpendicular to the fault based on section interaction and section gridding technology; 5) carrying out fault interpretation quality control by using a three-dimensional visualization technology; 6) after the interpretation of the single fault is finished, the established construction pattern library is used as prior information, and the combination relation between the faults on the section is restrained and guided by a classical construction pattern template; 7) repeating the steps 4) to 6) until the explanation of all fault in the whole work area is completed; 8) and establishing a fault frame model of the work area to finish fault interpretation work.
Description
Technical Field
The invention relates to a seismic fault interpretation method for adjusting a diving stage in the middle and later stages of oil and gas field development, in particular to a three-dimensional fault interpretation method for a complex fault block oil and gas field.
Background
The fault of the complex fault block oil and gas field develops for multiple times, and the fracture system is complex. After the oil field enters the middle and later development stages, along with the requirements of residual oil distribution prediction and well position optimization, the precision of seismic structure interpretation is also more strictly required. In recent years, with the progress of three-dimensional seismic data acquisition and processing technology, a full three-dimensional interpretation technology taking full-automatic interpretation of a three-dimensional space as a core is rapidly developed. The technology sets spatial seed points according to geological knowledge, and realizes the explanation of a horizon and a fault and the carving of the top and the bottom of a sand body through automatic tracking. However, in practical applications, since the applicable conditions are too ideal, the technical effect is not good and the practicability is poor when the data quality is poor or the geological conditions are complex. Under the background, the seismic structure interpretation returns to the track of the three-dimensional data two-dimensional interpretation again, the interpretation is carried out section by section, then the closure of the horizon and the fault is sought in the space, the main survey line and the connecting line after the three-dimensional work area is thinned are taken as a frame through line, surface and body ideas, and the survey line interpretation is gradually encrypted to realize the structure interpretation. For a large three-dimensional data volume, the problems of layer position, fault unclosed, unreasonable fault plane combination and the like can exist when the layer position, the fault and the fault plane are combined according to a two-dimensional interpretation method, the layer position and the fault plane are required to be repeatedly modified, and the working efficiency is seriously influenced. Therefore, how to effectively mine the intrinsic information contained in the three-dimensional data and improve the utilization rate and the working efficiency of the three-dimensional data information is a problem which needs to be solved urgently in the current seismic structure interpretation work.
Disclosure of Invention
In view of the above problems, the present invention aims to provide a fault interpretation method based on a combination of profile interaction and structural style guidance, which can provide high enough interpretation precision when fault interpretation is performed in a complex fault block oil field, ensure fault space closure, and improve work efficiency.
In order to achieve the purpose, the invention adopts the following technical scheme that a fault interpretation method based on combination of flattening interaction and construction style guidance is characterized by comprising the following steps of:
1) performing interpretative processing on the seismic data for fault interpretation to obtain a seismic data body capable of fully highlighting fault information;
2) utilizing a variance body and ant body technology, optimizing calculation parameters, and making attribute slices near an interpretation target layer to be used as a fault interpretation plane navigation chart;
3) analyzing a construction stress field and a profile typical construction pattern of a work area, establishing a construction pattern library, and using the construction pattern library as fault interpretation prior information and constraint;
4) carrying out fault interpretation on a section perpendicular to the fault based on section interaction and section gridding technology;
5) carrying out fault interpretation quality control by using a three-dimensional visualization technology to complete fault space closure;
6) after the interpretation of the single fault is finished, the construction pattern library established in the step 3) is used as prior information, and a fault interpretation pattern in a real area is implemented by using a classic construction pattern template in the construction pattern library to constrain and guide the combination relationship between faults on the section;
7) repeating the steps 4) to 6) until the explanation of all fault in the whole work area is completed;
8) and establishing a fault frame model of the work area to finish fault interpretation work.
Further, in step 1) above, the explanatory processes include resolution enhancement, noise suppression, and structure-guided structure filtering.
Further, in the step 2), the calculation parameters include a time window length, an operator radius, an operator shape and a seismic data volume; and (3) optimizing the calculation parameters by adopting a limit trial parameter method, firstly segmenting an optional parameter range by using a preset total sample number, and obtaining an optimal variance body by experimenting different values of the calculation parameters so as to determine the optimal value of the calculation parameters.
Further, in the step 3), a specific process of establishing the construction style library is as follows:
firstly, performing structural stress field analysis on a work area, researching a structural stress background, a main structural action, a structural position and a structural pattern of a basin where the work area is located, and determining basic structural elements and structural keywords contained in the work area; and searching according to the construction key words of the work area, determining the construction style of the work area, performing analogy with the construction characteristics of other oil fields in adjacent blocks, determining the typical construction style of the development of the area, and establishing a fault interpretation construction style library.
Further, in the step 4), fault interpretation is performed on a section perpendicular to the fault based on the profile interaction and the section gridding technology, and the specific process is as follows:
firstly, navigating to a region capable of clearly depicting a fault layer based on the planar navigation map formed in the step 2), and interpreting a fault on 2-3 adjacent attribute slices at preset time intervals;
secondly, the fault is explained on the section, the fault explanation result on the middle plane is projected on the seismic section in the first step, fault projection points are formed on the seismic section, and the fault explanation of the fault on the section is completed by connecting the projection points;
thirdly, performing section gridding prediction on the fault interpretation result on the section in the second step by using a section gridding technology so as to refine fault interpretation;
and further extracting the section along the direction vertical to the fault trend, generating a fault grid according to the fault interpretation result on the second-step section, and performing interpolation interpretation on the section further extracted along the direction vertical to the fault trend so as to refine the fault interpretation.
Further, in the first step, the predetermined time interval is 50 to 80 milliseconds.
Further, in the step 5), the three-dimensional visualization technology is used for fault interpretation quality control, and the specific content includes:
generating a section based on the fault interpretation result in the step 4), and in a three-dimensional view, detecting whether an abnormal interpretation line exists on the section and whether the section is smooth by naked eyes along the trend of a fault combination formed by the fault interpretation result in the plane in the step 4) and the interpretation result in the section in the second step, so as to judge whether the fault interpretation is reasonable;
if the section is not smooth, returning to the step 4), correcting the fault interpretation result on the section, re-interpreting the non-smooth fault line, and then generating the section;
and thirdly, repeating the steps of the first step and the second step until the fault surface is smooth, and finishing the fault space closing.
Further, in the step 6), the construction style library established in the step 3) is used as prior information to construct a combination relationship between classical construction style template constraints in the style library and faults on the guidance profile, and the specific content includes:
firstly, fixing the longitudinal and transverse proportion of section display, comparing the existing main section layer combination pattern in a window with the classic pattern in a construction pattern library window by window, and preferably selecting the classic pattern with the highest matching degree as an expected pattern;
secondly, correcting the cutting and combination relation of main faults on the section based on the expected pattern;
and adjusting the actually interpreted fault pattern according to the overlapping relation and the pattern form of the expected fault pattern, so that the difference between the actually interpreted fault pattern and the expected fault pattern is minimized, and the implementation of the section fault layer combination is completed.
By adopting the technical scheme, the invention has the following advantages: the method comprises the steps of performing interpretative processing on seismic data for fault interpretation, manufacturing attribute slices near an interpretation target layer by using a variance body and ant body technology, analyzing a structural stress field and a section typical structural pattern of a work area, establishing a structural pattern library, interpreting prior information and constraint for faults by using the structural pattern library, navigating, fault gridding and three-dimensional visualization as quality control by using a plane attribute slice, and realizing fault interpretation by means of profile interaction, so that high interpretation precision can be provided when fault interpretation is carried out on a complex fault block oil field, fault space closure is ensured, and working efficiency is improved.
Drawings
FIG. 1 is a schematic view of the flow structure of the present invention;
FIG. 2a is an explanatory pre-processed seismic data set at fault interpretation for a particular field, and FIG. 2b is an explanatory post-processed seismic data set;
FIG. 3 is a plan navigation attribute extracted from the seismic data after the explanatory processing at the time of interpretation of the oilfield fault;
FIG. 4 is a library of formation patterns created by formation pattern analysis when the oilfield faults are interpreted, FIG. 4a is a large spade positive fault, FIG. 4b is an X-type conjugate positive fault, FIG. 4c is a half-fancy slip fault (Y-shaped fault), FIG. 4d is a broken step positive fault, FIG. 4e is a negative fancy slip fault, and FIG. 4f is a spade half-flower fault;
FIG. 5 is a fault interpreted by the flattening interaction at the time of the oilfield fault interpretation, FIG. 5a is a plan fault interpretation, and FIG. 5b is a cross-sectional fault interpretation;
FIG. 6 is an explanation using fault gridding techniques in explaining the fault;
FIG. 7 is a quality control of fault interpretation using three-dimensional visualization techniques, FIG. 7a is an original fault interpretation, FIG. 7b is a corrected fault interpretation, FIG. 7c is an original fault plane, and FIG. 7d is a corrected section;
FIG. 8 is a combination of faults on a section under the constraint of prior information of a construction pattern library, FIG. 8A is a combination of a shovel type positive fault and a negative pattern fault, FIG. 8B is a combination graph of a fault order and a semi-pattern fault, FIG. 8C is a combination of a negative pattern fault and an X-type conjugate fault layer, and FIG. 8D is a combination of a shovel type semi-pattern fault and a fault order;
FIG. 9 is a fault framework model established after the interpretation of all faults in the work area is completed.
Detailed Description
The preferred embodiments of the present invention will be described in detail below with reference to the accompanying drawings so that the objects, features and advantages of the invention can be more clearly understood. It should be understood that the embodiments shown in the drawings are not intended to limit the scope of the present invention, but are merely intended to illustrate the spirit of the technical solution of the present invention.
As shown in FIG. 1, the invention provides a fault interpretation method based on combination of profile interaction and structural style guidance, which comprises the following steps:
1) performing interpretative processing on the seismic data for fault interpretation to obtain a seismic data body capable of fully highlighting fault information;
the explanatory processing comprises improving the resolution, suppressing noise, filtering a structure guide structure and the like so as to improve the data quality, enhance the transverse continuity of the seismic reflection event, highlight fault display and obtain a seismic data volume capable of fully highlighting fault information.
2) Making attribute slices near an interpretation target layer by using a variance body and ant body technology and optimizing calculation parameters, and taking the attribute slices as a fault interpretation plane navigation chart;
the generation of the plane navigation graph relates to a plurality of important parameters, and the important parameters mainly comprise time window length, operator radius, operator shape, seismic data volume and the like. The optimal value of the calculation parameter is optimized by adopting a limit trial parameter method, firstly, the optional parameter range is divided by the preset total sample number, and different values of the calculation parameter are tested to obtain the optimal variance body, so that the optimal value of the calculation parameter is determined.
The following examples are given for illustrative purposes: presetting the total sample number to be 4, preferably selecting the optimal value of the time window length, wherein the time window length selectable by software is 1-160 milliseconds, respectively testing the planar navigation map effects under the time window lengths of 40 milliseconds, 80 milliseconds, 120 milliseconds and 160 milliseconds at the interval of (160-1)/4-40 milliseconds, and selecting the time window length value corresponding to the optimal planar navigation map effect as the optimal value; then, the optimal values of other parameters are optimized in the manner.
3) Analyzing a construction stress field and a profile typical construction pattern of a work area, establishing a construction pattern library, and using the construction pattern library as fault interpretation prior information and constraint;
firstly, performing structural stress field analysis on a work area, researching a structural stress background, a main structural action, a structural position and a structural pattern of a basin where the work area is located, and determining basic structural elements and structural keywords contained in the work area; and searching according to the construction key words of the work area, determining the construction style of the work area, performing analogy with the construction characteristics of other oil fields in adjacent blocks, determining the typical construction style of the development of the area, and establishing a fault interpretation construction style library.
4) Carrying out fault interpretation on a section perpendicular to the fault trend direction based on section-level interaction (section-level interaction) and a section gridding technology; the method specifically comprises the following steps:
firstly, navigating to a region capable of clearly depicting a fault area based on the plane navigation map formed in the step 2), and interpreting a fault on 2-3 adjacent attribute slices at a preset time interval to form a fault interpretation result on a plane; wherein the predetermined time interval may be 50ms to 80 ms;
secondly, the fault is explained on the section, the fault explanation result on the middle plane is projected on the seismic section in the first step, fault projection points are formed on the seismic section, and the fault explanation of the fault on the section is completed by connecting the projection points;
thirdly, performing section gridding prediction on the fault interpretation result on the section in the second step by using a section gridding technology so as to refine fault interpretation;
and further extracting the section along the direction vertical to the fault trend, generating a fault grid according to the fault interpretation result on the second-step section, and performing interpolation interpretation on the section further extracted along the direction vertical to the fault trend so as to refine the fault interpretation.
5) Carrying out fault interpretation quality control by using a three-dimensional visualization technology to complete fault space closure;
generating a section based on the fault interpretation result in the step 4), and in a three-dimensional view, visually checking whether an abnormal interpretation line exists on the section and whether the section is smooth or not along the combined trend of the sections formed by the fault interpretation result in the plane in the step 4) and the interpretation result in the section in the second step, thereby judging whether the fault interpretation is reasonable or not.
If the section is not smooth, returning to the step 4), correcting the fault interpretation result on the section, re-interpreting the non-smooth fault line, and then generating the section.
And thirdly, repeating the steps of the first step and the second step until the fault surface is smooth, and finishing the fault space closing.
6) After the interpretation of the single fault is finished, the construction pattern library established in the step 3) is used as prior information, and a fault interpretation pattern in a real area is implemented by using a classic construction pattern template in the construction pattern library to constrain and guide the combination relationship between faults on the section; the method specifically comprises the following steps:
firstly, fixing the longitudinal and transverse proportion of section display, comparing the existing main section layer combination pattern in a window with the classic pattern in a construction pattern library window by window, and preferably selecting the classic pattern with the highest matching degree as an expected pattern;
secondly, correcting the cutting and combination relation of the main fault on the section based on the expected pattern,
and adjusting the actually interpreted fault pattern according to the overlapping relation and the pattern form of the expected fault pattern, so that the difference between the actually interpreted fault pattern and the expected fault pattern is minimized, and the implementation of the section fault layer combination is completed.
7) Repeating the steps 4) to 6) until the explanation of all fault in the whole work area is completed;
8) and establishing a fault frame model of the work area to finish fault interpretation work.
The invention is illustrated below with reference to specific examples:
the fault of the complex fault block oil and gas field develops for multiple times, and the fracture system is complex. The structural implementation degree directly determines the oil field development and adjustment effect, and is the key for residual oil distribution prediction and well design adjustment implementation. Taking a complex fault block oil field fault explanation entering the middle and later stages of development as an example, the fault explanation is developed by adopting the method of the invention, which comprises the following steps:
s1, performing interpretative processing such as resolution improvement, noise suppression, structural guide structure filtering and the like on the seismic data of the oil field for fault interpretation so as to improve the seismic data quality, enhance the transverse continuity of seismic reflection event, highlight fault display and obtain a seismic data body capable of fully highlighting fault information. The comparison of the seismic data before and after the explanatory process is shown in fig. 2, where fig. 2a is the seismic data before the explanatory process and fig. 2b is the seismic data after the explanatory process.
S2, combining the variance body and ant body technology, optimizing calculation parameters including time window length, calculation radius, operator shape, seismic data body and other influence factors, and making attribute slices near the interpretation target layer as a fault interpretation plane navigation chart, as shown in FIG. 3.
S3, researching the fracture background of the oil field, analyzing the section construction pattern, establishing a construction pattern library, and using the construction pattern library as fault interpretation prior information and constraint;
the oil field is located in a Bohai Bay basin which is a typical newly-born generation main basin, and a middle-newly-born generation superposed trap breaking basin has an early-stage stretching and superposing late-stage sliding structure stress field. In addition, the oil field structure is located on the boundary of a yellow river mouth recess and a Bohai south low bulge, and has a structure pattern that two mountains sandwich one deep valley. The Bohai Bay basin early stretching and later stage sliding construction stress background determines that the construction style contained in the oil field is formed by combining four positive fault layers (non-rotating plane type, shovel type and sloping plateau type positive fault layers) and two sliding fault layers (positive pattern type and negative pattern type sliding fault layers) basic construction elements. On the basis, the oil field structure keywords of 'yellow river mouth sunken', 'Bohai south low bulge' and 'two hill-in-one deep valley pattern' are researched and researched, and the oil field is locked to have 6 typical structure patterns (as shown in fig. 4), wherein fig. 4a is a large shovel-type positive fault, fig. 4b is an X-type conjugate positive fault, fig. 4c is a half-fancy slip fault (Y-shaped fault), fig. 4d is a broken step-type positive fault, fig. 4e is a negative fancy slip fault, and fig. 4f is a shovel-type half-flower fault.
And S4, developing true three-dimensional fault interpretation based on the interaction between the profile and the plane (profile-plane interaction).
Firstly, based on the planar navigation map obtained in step S2, a fault area capable of being clearly depicted is navigated, and a fault F1 is interpreted on 2 to 3 adjacent attribute slices at a time interval of 50ms to 80ms (as shown in fig. 5 a).
Secondly, explaining F1 on the section, projecting the fault F1 explained in the first step on the seismic section to form fault projection points on the seismic section, and connecting the projection points to finish the fault explanation of the fault F1 on the section (as shown in FIG. 5 b);
thirdly, performing section gridding prediction on the fault interpretation result on the section in the second step by using a section gridding technology so as to refine fault interpretation;
further extracting the section along the direction perpendicular to the fault F1, generating a fault grid according to the fault interpretation result on the second step section, and carrying out interpolation interpretation on the section further extracted along the direction perpendicular to the fault F1, thereby refining the fault interpretation (as shown in FIG. 6).
And S5, performing fault interpretation quality control by using a three-dimensional visualization technology.
A cross section of F1 is generated based on the fault interpretation result in step S4, and it is judged that there is an unreasonable interpretation of the fault F1 by finding with a visual inspection that there is an abnormal interpretation line and the cross section is not smooth along the direction of the combination of the fault layers formed by the fault interpretation result in the plane of the first step and the interpretation result on the cross section of the second step in step 4) in the three-dimensional view. Returning to the step 4), the F1 interpretation result on the corresponding section is corrected until the fault plane of the F1 is smooth. Fig. 7a shows the original fault interpretation, fig. 7b shows the corrected fault interpretation under the three-dimensional visualization technique, fig. 7c shows the original fault plane, and fig. 7d shows the corrected fault plane.
S6, after completing the interpretation of F1, using the priori information of the construction pattern library established in step S3 to constrain and guide the cutting and combination relationship between each fault on the section by using the specific 6 typical construction patterns in the construction pattern library, so as to implement the fault interpretation pattern in the working area (as shown in fig. 8), which roughly includes a shovel type positive fault + negative pattern fault combination shown in fig. 8A, a fault order + half pattern fault combination shown in fig. 8B, a negative pattern fault + X type conjugate fault combination shown in fig. 8C, and a shovel type half pattern fault + fault order combination shown in fig. 8D.
And S7, repeating the steps S4 to S6 until the explanation of all faults in the oil field full area is completed.
And S8, building the high-precision fault framework model of the oil field by using the interpreted fault, and completing fault interpretation work as shown in FIG. 9.
The present invention has been described with reference to the above embodiments, and the structure, arrangement, and connection of the respective members may be changed. On the basis of the technical scheme of the invention, the improvement or equivalent transformation of the individual components according to the principle of the invention is not excluded from the protection scope of the invention.
Claims (8)
1. A fault interpretation method based on combination of profile interaction and structural style guidance is characterized by comprising the following steps:
1) performing interpretative processing on the seismic data for fault interpretation to obtain a seismic data body capable of fully highlighting fault information;
2) utilizing a variance body and ant body technology, optimizing calculation parameters, and making attribute slices near an interpretation target layer to be used as a fault interpretation plane navigation chart;
3) analyzing a construction stress field and a profile typical construction pattern of a work area, establishing a construction pattern library, and using the construction pattern library as fault interpretation prior information and constraint;
4) carrying out fault interpretation on a section perpendicular to the fault based on section interaction and section gridding technology;
5) carrying out fault interpretation quality control by using a three-dimensional visualization technology to complete fault space closure;
6) after the interpretation of the single fault is finished, the construction pattern library established in the step 3) is used as prior information, and a fault interpretation pattern in a real area is implemented by using a classic construction pattern template in the construction pattern library to constrain and guide the combination relationship between faults on the section;
7) repeating the steps 4) to 6) until the explanation of all fault in the whole work area is completed;
8) and establishing a fault frame model of the work area to finish fault interpretation work.
2. The fault interpretation method based on the combination of the profile interaction and the construction style guidance as claimed in claim 1, wherein: in step 1) above, the explanatory processes include resolution enhancement, noise suppression, and structure-guided structure filtering.
3. The fault interpretation method based on the combination of the profile interaction and the construction style guidance as claimed in claim 1, wherein: in the step 2), the calculation parameters comprise time window length, operator radius, operator shape and seismic data volume; and (3) optimizing the calculation parameters by adopting a limit trial parameter method, firstly segmenting an optional parameter range by using a preset total sample number, and obtaining an optimal variance body by experimenting different values of the calculation parameters so as to determine the optimal value of the calculation parameters.
4. The fault interpretation method based on the combination of the profile interaction and the construction style guidance as claimed in claim 1, wherein in the step 3), the specific process of establishing the construction style library is as follows:
firstly, performing structural stress field analysis on a work area, researching a structural stress background, a main structural action, a structural position and a structural pattern of a basin where the work area is located, and determining basic structural elements and structural keywords contained in the work area; and searching according to the construction key words of the work area, determining the construction style of the work area, performing analogy with the construction characteristics of other oil fields in adjacent blocks, determining the typical construction style of the work area development, and establishing a fault interpretation construction style library.
5. The method for fault interpretation based on combination of profile interaction and structural style guidance as claimed in claim 1, wherein in the step 4), fault interpretation is performed on a section perpendicular to the fault based on profile interaction and section gridding technology, and the specific process is as follows:
firstly, navigating to a region capable of clearly depicting a fault layer based on the planar navigation map formed in the step 2), and interpreting a fault on 2-3 adjacent attribute slices at preset time intervals;
secondly, the fault is explained on the section, the fault explanation result on the middle plane is projected on the seismic section in the first step, fault projection points are formed on the seismic section, and the fault explanation of the fault on the section is completed by connecting the projection points;
thirdly, performing section gridding prediction on the fault interpretation result on the section in the second step by using a section gridding technology so as to refine fault interpretation;
and further extracting the section along the direction vertical to the fault trend, generating a fault grid according to the fault interpretation result on the second-step section, and performing interpolation interpretation on the section further extracted along the direction vertical to the fault trend so as to refine the fault interpretation.
6. The fault interpretation method based on the combination of dissection and leveling interaction and construction style guidance as claimed in claim 5, wherein: in the first step, the predetermined time interval is 50-80 milliseconds.
7. The fault interpretation method based on the combination of the dissection interaction and the construction style guidance as claimed in claim 5, wherein in the step 5), the fault interpretation quality control is performed by using a three-dimensional visualization technology, and the specific content includes:
generating a section based on the fault interpretation result in the step 4), and in a three-dimensional view, detecting whether an abnormal interpretation line exists on the section and whether the section is smooth by naked eyes along the trend of a fault combination formed by the fault interpretation result in the plane in the step 4) and the interpretation result in the section in the second step, so as to judge whether the fault interpretation is reasonable;
if the section is not smooth, returning to the step 4), correcting the fault interpretation result on the section, re-interpreting the non-smooth fault line, and then generating the section;
and thirdly, repeating the steps of the first step and the second step until the fault surface is smooth, and finishing the fault space closing.
8. The fault interpretation method based on the combination of profile interaction and construction pattern guidance as claimed in claim 1, wherein in the step 6), the construction pattern library established in the step 3) is used as prior information to construct a combination relationship between classical construction pattern template constraints in the pattern library and faults on the guidance profile, and the concrete contents include:
firstly, fixing the longitudinal and transverse proportion of section display, comparing the existing main section layer combination pattern in a window with the classic pattern in a construction pattern library window by window, and preferably selecting the classic pattern with the highest matching degree as an expected pattern;
secondly, correcting the cutting and combination relation of main faults on the section based on the expected pattern;
and adjusting the actually interpreted fault pattern according to the overlapping relation and the pattern form of the expected fault pattern, so that the difference between the actually interpreted fault pattern and the expected fault pattern is minimized, and the implementation of the section fault layer combination is completed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910831798.7A CN110488353B (en) | 2019-09-04 | 2019-09-04 | Fault interpretation method based on combination of profile interaction and construction style guidance |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910831798.7A CN110488353B (en) | 2019-09-04 | 2019-09-04 | Fault interpretation method based on combination of profile interaction and construction style guidance |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110488353A CN110488353A (en) | 2019-11-22 |
| CN110488353B true CN110488353B (en) | 2020-11-03 |
Family
ID=68556324
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910831798.7A Active CN110488353B (en) | 2019-09-04 | 2019-09-04 | Fault interpretation method based on combination of profile interaction and construction style guidance |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110488353B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111474583A (en) * | 2020-06-03 | 2020-07-31 | 中国石油化工股份有限公司 | Fault interpretation method and structural trap identification method for fault block oil reservoir |
| CN111856572B (en) * | 2020-07-06 | 2021-07-20 | 中国石油大学(北京) | Method and device for determining width of fault fracture belt |
| CN111948709B (en) * | 2020-07-30 | 2021-08-06 | 中国石油大学(北京) | Fracture zone structure identification method, device and system |
| CN112505757B (en) * | 2020-11-17 | 2021-08-10 | 中国石油集团工程咨询有限责任公司 | Three-dimensional seismic fault data processing method based on parallel grid control model |
| CN113933895A (en) * | 2021-10-29 | 2022-01-14 | 中国石油大学(华东) | Three-dimensional transparent display method for seismic interpretation horizon |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7069149B2 (en) * | 2001-12-14 | 2006-06-27 | Chevron U.S.A. Inc. | Process for interpreting faults from a fault-enhanced 3-dimensional seismic attribute volume |
| CN103941287B (en) * | 2014-04-16 | 2016-10-12 | 中国海洋石油总公司 | A kind of quick three-dimensional fault interpretation method based on horizontal navigation |
| CN104199099B (en) * | 2014-09-19 | 2016-08-24 | 东北石油大学 | A kind of method utilizing multiwindow serial section efficiently to carry out seismic interpretation |
| CN105334534A (en) * | 2015-10-21 | 2016-02-17 | 中国石油大学(华东) | Low order fault interpretation method based on construction mode guidance |
| CN107765301A (en) * | 2017-10-13 | 2018-03-06 | 中国煤炭地质总局地球物理勘探研究院 | The method for quickly identifying and device of coal seam craven fault |
-
2019
- 2019-09-04 CN CN201910831798.7A patent/CN110488353B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110488353A (en) | 2019-11-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110488353B (en) | Fault interpretation method based on combination of profile interaction and construction style guidance | |
| CN105469443B (en) | Based on geotraverse（PRB）The method that process Data fuse generates three-dimensional geological figure | |
| CN102918423B (en) | Method for earthquake hydrocarbon system anlysis | |
| CN107515423B (en) | Method for constructing three-dimensional geological model with high and steep structure | |
| CN109633745B (en) | Three-dimensional construction diagram drawing method and device | |
| NO328482B1 (en) | Method, product and system for automated, non-artificial, three-dimensional fault surface-based horizon modeling of the subsurface | |
| US10890688B2 (en) | Method for generating secondary data in geostatistics using observed data | |
| CN103941286A (en) | Rapid three-dimensional horizon interpretation method based on horizontal navigation | |
| AU2009264967A1 (en) | Well log correlation | |
| CN107728227B (en) | Method for rapidly distinguishing sand bodies of underground river in well pattern mature area | |
| CN112394404A (en) | Progressive reservoir fine characterization method | |
| CN108119133B (en) | Method and device for determining reservoir deposition process | |
| CN110428497B (en) | Braided river training image generation method | |
| CN105487116B (en) | FEM layer model method for building up | |
| CN110703334B (en) | Method for quickly explaining complex fault block oil field horizon | |
| CN102353990B (en) | Method for manufacturing seismic slice with overthrust fault | |
| CN101587597B (en) | Construction method of complex-structure block geological model based on geological rule constraint | |
| CN113821906A (en) | Single sand body period division method and multi-period sand body fine carving method | |
| CN111474583A (en) | Fault interpretation method and structural trap identification method for fault block oil reservoir | |
| CN110441821B (en) | Fault fast interpretation method based on variable interpretation net density | |
| CN114740547A (en) | Back-flushing sliding fracture fold grading and partitioning analysis and three-dimensional characterization method | |
| CN111815769B (en) | Modeling method, computing device and storage medium for thrust covered zone construction | |
| Kasmaee et al. | Risk reduction in Sechahun iron ore deposit by geological boundary modification using multiple indicator Kriging | |
| CN112782760B (en) | Method for dissecting braided river reservoir structure by using discontinuous boundaries of seismic reservoir | |
| Moore et al. | Evaluation of Nd isotope data for the Grenville Province of the Laurentian shield using a geographic information system |
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 |