CN117111143A - Method for realizing reliability of structure trap of mountain front punching belt - Google Patents
Method for realizing reliability of structure trap of mountain front punching belt Download PDFInfo
- Publication number
- CN117111143A CN117111143A CN202210533424.9A CN202210533424A CN117111143A CN 117111143 A CN117111143 A CN 117111143A CN 202210533424 A CN202210533424 A CN 202210533424A CN 117111143 A CN117111143 A CN 117111143A
- Authority
- CN
- China
- Prior art keywords
- layer
- speed
- geological
- trap
- reliability
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000004080 punching Methods 0.000 title claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 16
- 238000010586 diagram Methods 0.000 claims abstract description 13
- 238000006243 chemical reaction Methods 0.000 claims abstract description 11
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 10
- 238000005553 drilling Methods 0.000 claims description 31
- 238000009933 burial Methods 0.000 claims description 30
- 230000008859 change Effects 0.000 claims description 20
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000001228 spectrum Methods 0.000 claims description 11
- 238000009826 distribution Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000008021 deposition Effects 0.000 claims description 5
- 238000010276 construction Methods 0.000 abstract description 30
- 239000010410 layer Substances 0.000 description 169
- 238000013508 migration Methods 0.000 description 12
- 230000005012 migration Effects 0.000 description 12
- 238000005755 formation reaction Methods 0.000 description 6
- 238000013507 mapping Methods 0.000 description 5
- 235000015125 Sterculia urens Nutrition 0.000 description 4
- 240000001058 Sterculia urens Species 0.000 description 4
- 238000004587 chromatography analysis Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 241001412225 Firmiana simplex Species 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000002679 ablation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000005253 cladding Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 229910052900 illite Inorganic materials 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- VGIBGUSAECPPNB-UHFFFAOYSA-L nonaaluminum;magnesium;tripotassium;1,3-dioxido-2,4,5-trioxa-1,3-disilabicyclo[1.1.1]pentane;iron(2+);oxygen(2-);fluoride;hydroxide Chemical compound [OH-].[O-2].[O-2].[O-2].[O-2].[O-2].[F-].[Mg+2].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[K+].[K+].[K+].[Fe+2].O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2 VGIBGUSAECPPNB-UHFFFAOYSA-L 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- VYMDGNCVAMGZFE-UHFFFAOYSA-N phenylbutazonum Chemical compound O=C1C(CCCC)C(=O)N(C=2C=CC=CC=2)N1C1=CC=CC=C1 VYMDGNCVAMGZFE-UHFFFAOYSA-N 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
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/282—Application of seismic models, synthetic seismograms
-
- 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/643—Horizon tracking
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 technical method for realizing the reliability of a mountain front broken belt structure trap, which comprises the following steps: step 1, establishing an earthquake interpretation work area; step 2, utilizing geological outcrop, regional drainage and synthesis record calibration to determine stratum attribution; step 3, determining the structural form, defining stratum filling of blocks with different structures, and establishing a geological structure model of the region; step 4, establishing an initial velocity field V (x, y, z) of the whole layer system of the work area according to the geological structure model; step 5, correcting the velocity field V (x, y, z) to obtain a corrected velocity field V1 (x, y, z); step 6, interpreting the destination layer to be T 0 And carrying out deep conversion on the data to obtain a depth domain structure diagram. The technical method for realizing the reliability of the construction trap of the front-mountain broken belt effectively solves the problem of difficult implementation of the reliability of the construction trap of the front-mountain broken belt oil-gas exploration, has a propulsion significance for the front-mountain broken belt oil-gas exploration, has a large popularization and application prospect, and has remarkable economic and social benefits.
Description
Technical Field
The invention relates to the technical field of seismic exploration, in particular to a technical method for realizing reliability of a structure trap of a front-mountain broken belt.
Background
The mountain front belt is broken by multiple periods of broken belt development and is back-flushed, a plurality of old strata are exposed to the ground surface, the structural deformation is serious, the earthquake speed is influenced by the double complex characteristics of the ground surface and the underground, the speed spectrum energy mass is not concentrated enough, the speed trend of the complex part is difficult to master, the reliability of the structural trap is difficult to realize, and therefore, a geological structure model and a variable speed map are required to be combined, the reliability of the structural trap is realized, and a reliable basis is provided for drilling well positions.
The geological structure model characterizes the oil and gas reservoir elements and the reservoir key period by analyzing the structure superposition evolution characteristics of the formation unconformity surface, the fracture system, the structure-formation sequence, the sediment filling and the different structure layers of the research area.
Variable mapping is one of the key factors in implementing the trap configuration. The current research method of variable mapping is more, mainly comprising the steps of well point speed difference fitting, model inversion based, average speed field calculation from superimposed speed spectrum by a model chromatography method, and time depth conversion directly by a prestack depth migration speed field method, so as to obtain a final structural map. The method for utilizing the speed field of the single well speed difference and the model-based inversion method cannot accurately show the change rule of the underground speed in a work area with uneven well distribution or sparse well pattern density. The three-dimensional average speed calculated from the superimposed speed field can truly show the change rule of the underground speed under the condition of no well. The prestack depth migration can accurately realize trap, prestack depth data migration in a complex region in front of a mountain is better than poststack time migration and prestack time migration, and construction can be accurately realized. The pre-stack depth migration has a long processing period and is expensive. Aiming at a low exploration degree area with only two-dimensional measuring lines in a front-mountain punching zone, no better variable graph forming method exists at present, so that the reliability degree of construction trap is difficult to implement, and the progress of oil and gas exploration is restricted.
In application number: in CN201611019456.8, a method for modeling depth domain velocity for processing seismic data in the front of mountain is related. The method comprises the following steps: 1) Local single-point control is carried out on the prestack time migration root mean square speed by utilizing speed constraint inversion, then the prestack time migration root mean square speed is converted into a time domain layer speed by a DIX formula, then deep conversion is carried out, a depth domain layer speed body is obtained, smoothing processing is carried out, and a prestack depth migration depth domain initial layer speed model is obtained; 2) Inverting the optimized speed model by using grid chromatography; 3) Performing prestack depth migration on all data of a work area by using an optimized speed model to obtain a depth domain superposition data body, then scaling the depth domain superposition data body to a time domain, and performing construction interpretation in the time domain to obtain a construction model; 4) And performing chromatographic inversion optimization on the structural model by adopting well constraint to obtain the method. The method improves the precision of the depth domain velocity model of the depth migration imaging processing of the high steep structure of the mountain front belt, ensures that the precision of the homing of the high steep structure of the mountain front belt is higher, and improves the degree of well earthquake coincidence.
In application number: in the chinese patent application of cn2017101101002. X, a method for making a structure diagram under a belt by thrust and thrust before mountain is related to, comprising the following steps: 1) Calculating the average speed of stratum with normal deposition area in front of mountain; 2) The influence of the pushing body is not considered, the reflection time of the target layer is converted into depth, and the depth of the reflecting layer is obtained; 3) Calculating a thrust covering body correction amount; 4) Correcting the depth of the reflecting layer obtained in the step 2) by adopting the correction amount of the thrust cladding body obtained in the step 3), so as to obtain the corrected depth of the reflecting layer; 5) And manufacturing a structural diagram according to the corrected depth of the reflecting layer. The method for manufacturing the structure map under the mountain front thrust covered belt improves the mapping precision of the structure map under the mountain front thrust covered belt, solves the problem that the existing complex structure under the mountain front thrust covered belt is difficult to finely map, is suitable for mapping the mountain front thrust covered belt with low signal-to-noise ratio and weak reflection of three-dimensional seismic data, and is easy to popularize and use.
In application number: in the chinese patent application of cn201811171969.X, a method and a device for constructing a velocity field are related. The method comprises the following steps: calculating speed values of junctions of a plurality of well points and a plurality of geological interfaces; performing interpolation operation on the speed value of each geological interface juncture to obtain a first speed value of the geological interface; acquiring a second speed value of each geological interface based on a reference speed field obtained by prestack depth migration; constructing an error velocity field based on the first velocity value and the second velocity value; and correcting the reference speed field by using the error speed field to obtain a target speed field for time-depth conversion. By utilizing the method, the depth domain related data with accurate well points and reasonable geological rules can be obtained in the complex fractured basin.
The prior art is greatly different from the invention, and the technical problem which is needed to be solved by the user is not solved, so that the invention provides a novel technical method for realizing the reliability of the structure trap of the front-mountain broken belt.
Disclosure of Invention
The invention aims to provide a technical method for realizing the reliability of the structure trap of the front-mountain broken belt, which solves the problem of difficult implementation of the reliability of the structure trap of the front-mountain broken belt with a thrust meaning for the oil-gas exploration of the front-mountain broken belt.
The aim of the invention can be achieved by the following technical measures: the technical method for realizing the reliability of the structure trap of the front-mountain punching belt comprises the following steps:
step 1, establishing an earthquake interpretation work area;
step 2, utilizing geological outcrop, regional drainage and synthesis record calibration to determine stratum attribution;
step 3, determining the structural form, defining stratum filling of blocks with different structures, and establishing a geological structure model of the region;
step 4, establishing an initial velocity field V (x, y, z) of the whole layer system of the work area according to the geological structure model;
step 5, correcting the velocity field V (x, y, z) to obtain a corrected velocity field V1 (x, y, z);
step 6, interpreting the destination layer to be T 0 And carrying out deep conversion on the data to obtain a depth domain structure diagram.
The aim of the invention can be achieved by the following technical measures:
in step 1, a seismic interpretation work area is established in a LandMark interpretation system, and seismic data, namely two-dimensional seismic lines, are input to obtain a mountain front zone seismic data section view and a seismic network plan view.
In the step 2, the horizon attribution of the exposed stratum is clearly shown through geological outcrop calibration; and (3) through well drilling synthesis record calibration, the zone guiding layer clearly belongs to the layer position of the deep stratum.
In step 3, fine and reasonable fault and horizon interpretation is performed on the seismic section by combining with the structural style to obtain each target horizon interpretation T 0 And the data and fracture system is used for determining the structural form, defining stratum filling of different structural blocks and establishing a geological structure model of the region.
In step 4, error analysis is carried out on the original velocity spectrum, so that reasonable availability is ensured; obtaining an average layer velocity of each layer from the earth's surface to the target layer based on the logging data and the geological survey data; an initial velocity field V (x, y, z) of the full layer system of the work area is established according to the geological structure model.
In step 5, the velocity field V (x, y, z) is corrected using the target layer-by-layer velocity of the drilled and virtual control points to obtain a corrected velocity field V1 (x, y, z).
In step 5, correcting a velocity field by using the velocity of the drilling layer in the adjacent region for the region with only two-dimensional seismic lines; the longitudinal change of the drilling stratum speed and the transverse change of the reservoir are used for constraining and correcting the speed field V (x, y, z), so that the longitudinal change and the transverse change of the stratum speed are ensured to accord with geological rules, and the corrected speed field V1 (x, y, z) is obtained; the drilling layer speed is calculated by: v=ht/2
Wherein H is the burial depth of the drilling target layer, and t is the double-journey travel time of the target layer.
In step 5, for the region with only two-dimensional seismic lines, selecting a virtual control point on a key position of the geological structure change according to the geological structure characteristics of the research region, calculating the layer speed of the point, and adding the upper limit and the lower limit of the stratum statistical speed to restrict and correct the speed field V1 (x, y, z).
In step 5, the selection of the key positions of the virtual control points is arranged according to the destination layer burial depths of the blocks with different structures, and if the burial depths of the two blocks are greatly different, the control points can be respectively arranged at the two blocks.
In step 5, the virtual point speed obtains the layer speed based on the target layer and the thickness of each set of stratum calculated when the two-way travel of each set of geological reflection layer measured at the virtual control point of the target layer is performed, the thickness of each set of stratum from the earth surface to the target layer is accumulated, the burial depth of the target layer is obtained, and the average speed of the target layer at the virtual well is obtained after the burial depth is divided by the one-way travel.
In step 5, the calculation formula of the thickness of the virtual control point layer of each geological structure layer is as follows:
ΔH i ＝V i *Δt i the method comprises the steps of (1),
in formula 1, ΔH i The thickness of the ith layer is the virtual control point; v (V) i Layer speed for the corresponding i-th layer; Δt (delta t) i Travel time differences for each set of geological reflecting layers;
H i ＝ΔH i +ΔH i-1 +ΔH i-2 + … +ΔH type 2
In formula 2, H i Is the burial depth of the ith layer at the virtual control point, delta H i The thickness is corresponding to the i layer;
in the formula 3, the components are mixed,the average layer speed corresponding to the ith layer at the virtual control point is obtained; t is t i Is the double journey travel time of the ith layer at the virtual control point.
In step 5, the average speed of the destination layer is obtained based on the average speed of the destination layer of the virtual control point, and the layer speed of the destination layer corresponding to the drilling of the neighboring cell can be used; when the reservoir is not changed greatly and the destination layer is buried quite deeply, the reservoir speed is not changed greatly at the position with simple structure.
In step 5, overall analysis is carried out on the layer speed, and whether the speed trend is reasonable or not is analyzed; whether the distribution of the layer speed on the plane is consistent with the relation of the deposition environment, the structure and the burial depth or not; after the layer speed is ensured to accord with the geological rule in the longitudinal and transverse directions, the corrected speed field V1 (x, y, z) is considered to be a reasonable available speed field.
In step 6, T will be explained 0 Performing time-depth conversion on the data and the corrected average layer velocity field V1 (x, y, z) to obtain depth domain horizon data H (x, y, z); and drawing a depth domain structure diagram together with the depth domain horizon data H (x, y, z) and the corresponding target layer fracture system.
The technical method for realizing the reliability of the structure trap of the mountain front broken belt further comprises the steps of after the step 6, checking whether the position, the amplitude and the space homing of the structure trap high points of the structure diagram are consistent with the structure form of the geological structure model, if the structure trap high points, the amplitude and the form are inconsistent with the geological structure model, returning to the step 5, and revising the speed field by adjusting the position and the speed of the virtual control point until the structure trap amplitude and the form of the target layer are close to reality; a better depth domain structure map is obtained.
The technical method for realizing the reliability of the mountain front broken belt structure trap is simple and convenient, and can effectively realize the reliability of the structure trap by utilizing the geological structure model to restrict the longitudinal and transverse changes of the stratum speed and utilizing the drilling logging speed and increasing the change rule of the virtual control point speed to restrict the plane speed. The technical method for realizing the reliability of the mountain front broken belt structure trap is characterized by establishing a geological structure model of a work area and utilizing the geological structure model to restrict the longitudinal and transverse changes of a speed field. The planar change of the layer speed is constrained by the method of drilling wells in the work area (adjacent area) or adding virtual control points. The method is a technical method for realizing construction trap reliability implementation by continuously correcting a refined speed field. The method effectively solves the problem of difficult implementation of reliability of the trap of the mountain front broken belt oil and gas exploration structure, has a propulsion meaning for mountain front belt oil and gas exploration, has a large popularization and application prospect, and has remarkable economic and social benefits.
Drawings
FIG. 1 is a flow chart of an embodiment of a method for implementing reliability of a mountain front break strap configuration trap in accordance with the present invention;
FIG. 2 is a cross-sectional view of a seismic interpretation of a geological outcrop of the line Z10N_ASL_NE008 in an Aldrich region according to an embodiment of the invention;
FIG. 3 is a cross-sectional view of an ASL008_NS2002_1 seismic interpretation of an Aldrich region in accordance with an embodiment of the invention;
FIG. 4 is a cross-sectional view of the geological structure of the Ashri region Z10_ASL_NE007 in accordance with an embodiment of the present invention;
FIG. 5 is a schematic view of the bottom surface of an eight-bay set of the dwarf system in the Ashli region in accordance with an embodiment of the present invention;
FIG. 6 is a bottom view of a set of two-stack karaya grooves in the Aldrich region according to one 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.
The invention discloses a technical method for realizing reliability of a mountain front punching belt structure trap. The planar change of the layer speed is constrained by the method of drilling wells in the work area (adjacent area) or adding virtual control points. The method is a technical method for realizing construction trap reliability implementation by continuously correcting a refined speed field. The method effectively solves the problem of difficult implementation of reliability of the trap of the mountain front broken belt oil and gas exploration structure, has a propulsion meaning for mountain front belt oil and gas exploration, has a large popularization and application prospect, and has remarkable economic and social benefits.
The following are several embodiments of the invention
Example 1
In a specific embodiment 1 of the present invention, the method for implementing the reliability of the front-mountain break belt structure trap includes the following steps (fig. 1):
1. and (3) establishing a seismic interpretation work area in the LandMark interpretation system, and inputting seismic data (two-dimensional seismic lines) to obtain a front-mountain break zone seismic data section view and a seismic network plan view.
2. The stratum attribution is defined by using geological outcrop, regional drainage and synthesis record calibration, and the step is the basis of the invention. The horizon attribution of the exposed stratum is clearly shown through geological outcrop calibration; the well drilling synthesis record is calibrated, and the zone guiding layer definitely belongs to the layer position of the deep stratum.
3. Performing fine and reasonable fault and horizon interpretation on the seismic section by combining with the structural style to obtain each target horizon interpretation T 0 And the data and fracture system is used for determining the structural form, defining stratum filling of different structural blocks and establishing a geological structure model of the region.
4. And carrying out error analysis on the original velocity spectrum to ensure reasonable availability. Based on the logging data and the geological survey data, an average layer velocity is obtained for each layer from the surface to the target layer. An initial velocity field V (x, y, z) of the full layer system of the work area is established according to the geological structure model.
5. The velocity field V (x, y, z) is corrected using the destination layer velocity of the drilled and virtual control points to obtain a corrected velocity field V1 (x, y, z).
6. In the DoubleFox geologic mapping software, T will be explained 0 The data and the corrected average layer velocity field V1 (x, y, z) are subjected to time-depth conversion to obtain depth domain horizon data H (x, y, z). And drawing a depth domain structure diagram together with the depth domain horizon data H (x, y, z) and the corresponding target layer fracture system.
7. Checking whether the position, the amplitude and the space homing of the construction trap high points of the construction diagram are consistent with the construction form of the geological structure model, if the construction trap high points, the amplitude and the form are inconsistent with the geological structure model, returning to the step 5, and revising the speed field by adjusting the position and the speed of the virtual control points until the construction trap amplitude and the form of the target layer approach reality. A better depth domain structure map is obtained.
Example 2
In an embodiment 2 of the present invention, the velocity field V (x, y, z) is modified by using the target layer-by-layer velocity of the drilled and virtual control points to obtain a modified velocity field V1.
And correcting the speed field by using the adjacent well drilling layer speed aiming at the area with only two-dimensional seismic survey lines. The longitudinal change of the drilling stratum speed and the transverse change of the reservoir are used for constraining and correcting the speed field V (x, y, z), so that the longitudinal change and the transverse change of the stratum speed are ensured to accord with geological rules, and the corrected speed field V1 (x, y, z) is obtained.
The drilling layer speed is calculated by: v=ht/2
Wherein H is the burial depth of the drilling target layer, and t is the double-journey travel time of the target layer.
For the area with only two-dimensional seismic lines, selecting a virtual control point on a key position of geological structure change according to geological structure characteristics of a research area, calculating the layer speed of the point, and adding upper and lower limit constraint correction speed fields V1 (x, y, z) of stratum statistical speed.
The selection of the key positions of the virtual control points is arranged according to the destination layer burial depths of blocks with different structures, and if the burial depths of two blocks are greatly different, the control points can be respectively arranged at the two blocks.
And the virtual point speed is obtained based on the layer speed of the target layer and the thickness of each set of stratum is calculated when the target layer travels in double-pass of each set of geological reflection layers measured at the virtual control point of the target layer, the thickness of each set of stratum from the earth surface to the target layer is accumulated, the burial depth of the target layer is obtained, and the average speed of the target layer at the virtual well is obtained when the burial depth is divided by the single-pass travel.
(1) The calculation formula of the thickness of the virtual control point layer of each geological structure layer is as follows:
ΔH i ＝V i *Δt i the method comprises the steps of (1),
in formula 1, ΔH i I is 1,2,3,4, etc. for the thickness of the i-th layer of the virtual control point; v (V) i Layer speed for the corresponding i-th layer; Δt (delta t) i Travel time differences for each set of geological reflection layers.
(2)H i ＝ΔH i +ΔH i-1 +ΔH i-2 + … +ΔH type 2
In formula 2, H i Is the burial depth of the ith layer at the virtual control point, delta H i The thickness corresponding to the i-th layer.
(3)
In the formula 3, the components are mixed,the average layer speed corresponding to the ith layer at the virtual control point is obtained; t is t i Is the double journey travel time of the ith layer at the virtual control point.
The average speed of the target layer is obtained based on the average speed of the target layer of the virtual control point, and the layer speed of the target layer corresponding to the drilling of the adjacent area can be used. When the reservoir is not changed greatly and the destination layer is buried quite deeply, the reservoir speed is not changed greatly at the position with simple structure.
And then carrying out overall analysis on the layer speed, and analyzing whether the speed trend is reasonable. Whether the distribution of layer velocity on a plane is consistent with the relationship of deposition environment, construction and burial depth. After the layer speed is ensured to accord with the geological rule in the longitudinal and transverse directions, the corrected speed field V1 (x, y, z) is considered to be a reasonable available speed field.
Example 3
In a specific embodiment 3 of the present invention, taking the example of the front impact zone of the pseudo-south mountain as the front fold impact zone of the pseudo-south mountain, the front fold impact zone of the pseudo-south mountain located in the illite black biped root mountain and the bogida mountain, the main geological structure is the back slope of the first Tun river, the back slope of the first Tun river is divided into a large-sized reverse mask-pushing covering structure by fracture, the north is steep and the south is steep, and the thickness distribution of the two-wing stratum is uneven.
Before the invention is implemented, the structural stress state and the structural evolution background of the Alshili region are studied, and the structural style of the region is determined so as to guide the establishment of a geological structure model.
Step 101, in a Landmark interpretation system, a seismic interpretation work area of an Aldrich region is established, and a two-dimensional seismic survey line is input to obtain a section view of seismic data of a front-mountain break zone and a plane view of a seismic survey network of the Aldrich region.
And 102, determining stratum attribution by using geological outcrop, regional drainage and synthesis record calibration, wherein the step is the basis of the invention. The horizon attribution of the exposed stratum is clearly shown through geological outcrop calibration; the well drilling synthesis record is calibrated, and the zone guiding layer definitely belongs to the layer position of the deep stratum.
In the step 102, determining the stratum attribution by using the geological outcrop area layer guiding and the synthetic record calibration may include: and (5) the outcrop on the earth surface in the quasi-south area is manufactured into a graph cut section by observing the field geological section. And (3) splicing, comparing and analyzing the outcrop map section adjacent to the working area survey line with the surface outcrop stratum on the working area two-dimensional seismic survey line, and realizing the horizon attribution of the surface outcrop stratum on the working area seismic survey line (figure 2).
In the step 102, the method may further include: and (3) drilling synthesis record calibration is carried out in a land mark workstation Synthool module, the time-depth relation and the earthquake-geological reflection characteristic corresponding relation of each target layer are defined, and the attribution and reflection characteristics of the lower layer are defined through the drainage layer of the continuous well section area (figure 3).
In the step 102, the geological-seismic correspondence between the earth surface and the stratum of the target layer in the ashray region can be obtained through the geological outcrop and the well connection calibration. Meanwhile, the positions and fault interpretation can be provided by identifying and marking the seismic reflection points and fault break points of each set of stratum on the seismic section, such as ablation, superelevation, top elevation, cutting and the like (figure 4).
Step 103, carrying out fine and reasonable fault and horizon interpretation on the seismic section by combining with the structural style to obtain horizon interpretation T of each stratum from the earth surface to the target layer 0 The data and fracture system determines the structural form, determines the stratum filling of different structural blocks and establishes a geological structure model of the region (figure 3).
In step 103, fault and horizon interpretation is performed in an interpretation system to build a geological structure model of the work area, including: in a seiworks module of a Landmark interpretation system, interpreting the layers of each stratum from the fracture and the ground surface to the target layer according to the investigated regional construction style to obtain a fracture system and each set of target layer interpretation T 0 Data. And establishing a reasonable geological structure model according to the fracture system and horizon data.
In step 103, fault and horizon interpretation is performed in an interpretation system to build a geological structure model of the work area, including: and a Fault module on a mapping base diagram of the Landmark interpretation system edits Fault Polygons, performs Fault plane combination, and checks Fault combination relations of all target layers from old to new in a flat-profile combined way, so that the implementation reliability of a fracture system is ensured.
And 104, checking and comparing analysis on the original velocity spectrum, and obtaining the layer velocity of each set of stratum from the earth surface of the work area to the target layer based on logging data and geological investigation. An initial velocity field V (x, y, z) of the whole layer system of the work area is established according to the geological structure model. The geologic structure model controls the longitudinal and transverse variation of the layer velocity.
In step 104, an initial velocity field of the work area is established from the geologic structure model, comprising: the layer velocities of the various formations in the ashlar region are obtained based on intra-zone (neighbor) drilling log data, seismic data, and geological survey data. No well is drilled in the ashi region, and only the drilling layer speed of the adjacent region can be used. Estimating the dwarf system speed interval of the target layer at 2800-3600 m/s according to the ashi and peripheral geological statistics and experience; the speed interval of the triad is 3600-4200 m/s; the speed interval of the two-fold system phoenix tree ditch group is 3300-4200 m/s; the speed range of the carboloy is 4800-5500m/s.
In step 104, a work area initial velocity field is established according to the geological structure model, and the method further comprises: and (3) checking, comparing and analyzing the original speed spectrum of the Aldrich region, wherein the Aldrich region only has 10 two-dimensional seismic lines, the speed spectrum data is poor in quality, and the rationality of the original speed spectrum value is ensured by checking the speed spectrum error of the closing point of the two-dimensional seismic lines.
In step 104, a work area initial velocity field is established according to the geological structure model, and the method further comprises: and (3) establishing a stratum lattice according to the contact relation of the stratum from old to new from bottom to top by using the well explained geological structure model, filling the stratum lattice with the stratum velocity of the corresponding stratum, and establishing a full-stratum velocity structure model of the Ashli region.
In step 104, a work area initial velocity field is established according to the geological structure model, and the method further comprises: and (3) carrying out model chromatographic velocity analysis by using loop wave software, and establishing an initial velocity field V (x, y, z) of the work area according to the model chromatographic velocity analysis.
The filling speed in the grid is reversely pushed from top to bottom layer by layer according to the principle of a model chromatography by utilizing the original superposition speed spectrum data. The principle is that when the layer speed of the n-1 layer and the n-1 th speed interface are known, the layer speed of the n-1 th layer is obtained through iteration and the position of the n-th reflecting interface is determined by utilizing the principle that the angle of incidence is equal to the angle of reflection through curved ray tracking, and meanwhile, the position of the reflecting point deviating from the incident point is obtained, and the like, and layer-by-layer calculation is performed.
The model chromatography is used for obtaining the layer speed, so that the influence of the stratum with a large overlying dip angle on the structure form and the structure high point can be eliminated, and the real structure form can be obtained. However, due to the complex structure of the earth surface and the underground of the mountain front belt, the single layer speed is difficult to meet the requirement of realizing the reliability of the construction trap.
Step 105, the velocity field is corrected by the target layer velocity of the drilled and virtual control points to obtain a corrected velocity field V1 (x, y, z).
In step 105, the velocity field is modified using the target layer-by-layer velocity of the drilled and virtual control points, including: the method is characterized in that the drilling and logging layer speeds in the Ashli region are researched and analyzed, the same layer is different in stratum speed according to the differences of stratum lithology, construction positions and stratum burial depths, and therefore a section value exists in the stratum speed. The drilled reference well can be selected according to lithology combinations of different structural blocks, the buried depth position and the distance.
In step 105, the velocity field is corrected using the target layer-by-layer velocity of the well and virtual control points, further comprising: according to the structural characteristics of the Aldrich region, the adjacent region exploratory well small 3-well drilling stratum is equivalent to the Aldrich region through region layer-guiding contrast analysis, the structural background is similar, the distance is relatively short, and reference can be made to the contrast analysis. The eight-channel bay group in the Aldrich area develops a plaited river delta, the reservoir lithofacies is characterized by the combination of gravel-containing fine sandstone, fine sandstone and gray mudstone, the distribution is stable, and the burial depth of the eight-channel bay group in the area is less influenced by the impact of the break structure. The initial velocity field V1 (x, y, z) can thus be modified by using the average bed velocity 3560m/s of the 3-well octal group.
In step 105, the velocity field is corrected using the target layer-by-layer velocity of the well and virtual control points, further comprising: and combining the characteristics of deep fracture slipping structures, and arranging virtual control points at the key positions of the structures for controlling different burial depths of the same stratum according to the characteristics of different burial depths of the fracture stratum. According to the construction form and the survey line distribution of the Aldrich area, 5-10 virtual control points are arranged according to the burial depth of a target layer, the filling characteristics of block strata with different constructions and the ground surface development condition, and the morphological characteristics of construction trap are controlled (table 1).
TABLE 1 virtual control Point coordinates and destination layer speed filling Table in Alshi region
In step 105, the velocity field is corrected using the target layer-by-layer velocity of the well and virtual control points, further comprising: the virtual control point average layer speed acquisition is based on theIs measured at the virtual control point of the target layer 0 And calculating to obtain the thickness of each stratum, accumulating the thickness of each stratum from the earth surface to the target layer to obtain the burial depth of the target layer, and dividing the burial depth by the average speed of the target layer at the virtual control point when traveling in one pass.
(1) The calculation formula of the thickness of the virtual control point layer of each geological structure layer is as follows:
ΔH i ＝V i *Δt i the method comprises the steps of (1),
in formula 1, ΔH i I is 1,2,3,4, etc. for the thickness of the i-th layer of the virtual control point; v (V) i Layer speed for the corresponding i-th layer; Δt (delta t) i Travel time differences for each set of geological reflection layers.
(2)，H i ＝ΔH i +ΔH i-1 +ΔH i-2 + … +ΔH type 2
In formula 2, H i Is the burial depth of the ith layer at the virtual control point, delta H i I is 1,2,3,4, etc. for the thickness corresponding to the i-th layer.
(3)，
In the formula 3, the components are mixed,the average speed corresponding to the ith layer at the virtual control point; t is t i Is the double journey travel time of the ith layer at the virtual control point.
In step 105, the virtual control point average velocity acquisition may also be based on the average formation velocity corresponding to the similar formation of the desired formation relative to the borehole. The speed interval of the two-stack karaya ditch group by using the small 3-well drilling layer speed is 3300-4200 m/s, and the average layer speed is 3800m/s.
In step 105, the velocity field is corrected using the target layer-by-layer velocity of the well and virtual control points, further comprising: the initial velocity field is modified using the average layer velocity at the virtual control point to obtain a modified velocity field V1 (x, y, z).
In step 105, overall analysis is performed on the layer speed, whether the trend of the speed is reasonable is analyzed, and the relationship between the layer speed plane distribution and the deposition environment, structure and burial depth is analyzed on the layer speed plane graph. After the layer speed is ensured to accord with the geological rule in the longitudinal and transverse directions, the corrected speed field V1 (x, y, z) is considered to be a reasonable available speed field.
Step 106, using the interpreted destination layer level T 0 T of data and fault combined sketch 0 The diagram will explain the good destination layer T 0 The data is subjected to time-depth conversion by using the corrected velocity field V1 (x, y, z), and a depth domain structure diagram is obtained.
In step 106, further includes: in the Doublefox geologic mapping software, the interpretation horizon T of the target layer is utilized 0 T of data and fault plane combined drawing target layer 0 And (5) constructing a drawing.
In step 106, it includes: for the destination layer of the Ashtray set bottom, the two-stack karaya groove set T in the Ashli region by using the corrected average velocity field V1 (x, y, z) 0 The data is subjected to time-depth conversion to obtain depth domain horizon data H (x, y, z).
In step 106, further includes: in the DoubleFox geological mapping software, depth domain structure maps are drawn by applying the depth domain data H (x, y, z) of the target layer and the corresponding fault plane combination of the target layer, and the depth domain structure maps of the Bay group and the two-fold karaya ditch group of the Ashli region are obtained.
And 107, constructing trap reliability implementation. Checking whether the position, the amplitude and the form of the high point of the construction map trap are consistent with the construction form and the amplitude of the geological structure model, obtaining a better depth domain construction map if the position, the amplitude and the form of the high point of the construction map trap are consistent with the construction form and the amplitude of the geological structure model, and returning to correct the velocity field V1 (x, y, z) again if the position, the amplitude and the form of the high point of the construction map trap are inconsistent with the construction form and the amplitude of the geological structure model until the construction map of the better depth domain is obtained.
The technical method for realizing the reliability of the structure trap of the front-mountain broken belt achieves good geological effects in the application of an Alshi region, and well solves the problem that the reliability of the structure trap of the front-mountain low exploration degree region is difficult to realize.
By using the technical method, the construction trap of each target layer in the Aldrich area is effectively realized, and the risk exploratory well 1 is successfully deployed in the area.
Drilling confirmed the reliability of the construction trap in the ashira region, the design depth of the target layer Badao bay group was 250m, the real drilling depth was 253m, the depth error was 3m, and the relative error was only 1.2% (fig. 5). The phoenix tree ditch group was designed to have a well depth of 1450m, a true well depth of 1587m, and a relative error of 9% (fig. 6). The error value is larger than that of the shallow layer due to the development of deep layer break.
Overall, the application of the technical method for realizing the reliability of the mountain front broken belt structure trap achieves good effects.
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 (15)
1. The technical method for realizing the reliability of the structure trap of the front-mountain punching belt is characterized by comprising the following steps:
step 1, establishing an earthquake interpretation work area;
step 2, utilizing geological outcrop, regional drainage and synthesis record calibration to determine stratum attribution;
step 3, determining the structural form, defining stratum filling of blocks with different structures, and establishing a geological structure model of the region;
step 4, establishing an initial velocity field V (x, y, z) of the whole layer system of the work area according to the geological structure model;
step 5, correcting the velocity field V (x, y, z) to obtain a corrected velocity field V1 (x, y, z);
step 6, interpreting the destination layer to be T 0 And carrying out deep conversion on the data to obtain a depth domain structure diagram.
2. The method for realizing the reliability of the mountain front broken belt structure trap according to claim 1, wherein in step 1, a seismic interpretation work area is established in a Landmark interpretation system, and seismic data, namely two-dimensional seismic lines, are input to obtain a mountain front belt seismic data section view and a seismic survey net plan view.
3. The method for implementing the reliability of the mountain front broken belt structure trap according to claim 1, wherein in step 2, the horizon attribution of the exposed stratum is clearly shown through geological outcrop calibration; and (3) through well drilling synthesis record calibration, the zone guiding layer clearly belongs to the layer position of the deep stratum.
4. The method for realizing the reliability of a mountain front break belt structure trap as claimed in claim 1, wherein in step 3, fine and reasonable fault and horizon interpretation is performed by combining the structure patterns on the seismic section to obtain each target layer interpretation T 0 And the data and fracture system is used for determining the structural form, defining stratum filling of different structural blocks and establishing a geological structure model of the region.
5. The method for implementing the reliability of the structure trap of the pre-mountain break belt according to claim 1, wherein in step 4, an error analysis is performed on an original velocity spectrum to ensure reasonable availability; obtaining an average layer velocity of each layer from the earth's surface to the target layer based on the logging data and the geological survey data; an initial velocity field V (x, y, z) of the full layer system of the work area is established according to the geological structure model.
6. The method of claim 1, wherein in step 5, the velocity field V (x, y, z) is corrected by the target layer velocity of the drilled and virtual control points to obtain a corrected velocity field V1 (x, y, z).
7. The method of implementing reliability of a mountain front broken string structure trap as claimed in claim 6, wherein in step 5, the velocity field is corrected by using the adjacent zone drilling layer velocity for the two-dimensional seismic line area only; the longitudinal change of the drilling stratum speed and the transverse change of the reservoir are used for constraining and correcting the speed field V (x, y, z), so that the longitudinal change and the transverse change of the stratum speed are ensured to accord with geological rules, and the corrected speed field V1 (x, y, z) is obtained; the drilling layer speed is calculated by: v=ht/2
Wherein H is the burial depth of the drilling target layer, and t is the double-journey travel time of the target layer.
8. The method of claim 6, wherein in step 5, for the two-dimensional seismic survey line area only, virtual control points are selected at key locations of the geologic structure change according to the geologic structure characteristics of the survey area, and the corrected velocity field V1 (x, y, z) is constrained by calculating the layer velocity of the points and adding the upper and lower limits of the formation statistical velocity.
9. The method according to claim 8, wherein in step 5, the selection of the key positions of the virtual control points is based on the destination layer burial depth arrangement of the blocks of different structures, and if the burial depths of the two blocks are greatly different, the control points can be respectively arranged at the two blocks.
10. The method for realizing the reliability of the mountain front broken belt structure trap according to claim 8, wherein in step 5, the virtual point speed is obtained, the thickness of each set of stratum is calculated based on the layer speed of the target layer and the double travel time of each set of geological reflection layers measured at the virtual control point of the target layer, the thicknesses of each set of stratum from the earth surface to the target layer are accumulated, the burial depth of the target layer is obtained, and the average speed of the target layer at the virtual well is obtained by dividing the burial depth by the single travel time.
11. The method for implementing reliability of a mountain front break belt structure trap as claimed in claim 10, wherein in step 5, the calculation formula of the virtual control point layer thickness of each geological structure layer is:
ΔH i ＝V i *Δt i the method comprises the steps of (1),
in formula 1, ΔH i The thickness of the ith layer is the virtual control point; v (V) i Layer speed for the corresponding i-th layer; Δt (delta t) i Travel time differences for each set of geological reflecting layers;
H i ＝ΔH i +ΔH i-1 +ΔH i-2 + … +ΔH type 2
In formula 2, H i Is the burial depth of the ith layer at the virtual control point, delta H i The thickness is corresponding to the i layer;
in the formula 3, the components are mixed,the average layer speed corresponding to the ith layer at the virtual control point is obtained; t is t i Is the double journey travel time of the ith layer at the virtual control point.
12. The method for implementing reliability of a mountain front broken belt structure trap according to claim 10, wherein in step 5, the average speed of the destination layer is obtained based on the average speed of the destination layer of the virtual control point, and the layer speed of the destination layer corresponding to the drilling well in the neighboring area can be used; when the reservoir is not changed greatly and the destination layer is buried quite deeply, the reservoir speed is not changed greatly at the position with simple structure.
13. The method for implementing reliability of a mountain front broken belt structure trap as claimed in claim 6, wherein in step 5, overall analysis is performed on the layer speed to analyze whether the speed trend is reasonable; whether the distribution of the layer speed on the plane is consistent with the relation of the deposition environment, the structure and the burial depth or not; after the layer speed is ensured to accord with the geological rule in the longitudinal and transverse directions, the corrected speed field V1 (x, y, z) is considered to be a reasonable available speed field.
14. The method for implementing reliability of a mountain break belt structure trap as claimed in claim 1, wherein in step 6, T is explained 0 Performing time-depth conversion on the data and the corrected average layer velocity field V1 (x, y, z) to obtain depth domain horizon data H (x, y, z); and drawing a depth domain structure diagram together with the depth domain horizon data H (x, y, z) and the corresponding target layer fracture system.
15. The method for realizing the reliability of the structure trap of the mountain front punching belt according to claim 1, wherein the method further comprises the steps of checking whether the position, the amplitude and the space homing of the structure trap high points of the structure chart are consistent with the structure form of a geological structure model after the step 6, and if the structure trap high points, the amplitude and the form are inconsistent with the geological structure model, returning to the step 5 to revise the speed field by adjusting the position and the speed of the virtual control point until the structure trap amplitude and the form of a target layer are close to reality; a better depth domain structure map is obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210533424.9A CN117111143A (en) | 2022-05-17 | 2022-05-17 | Method for realizing reliability of structure trap of mountain front punching belt |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210533424.9A CN117111143A (en) | 2022-05-17 | 2022-05-17 | Method for realizing reliability of structure trap of mountain front punching belt |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117111143A true CN117111143A (en) | 2023-11-24 |
Family
ID=88811571
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210533424.9A Pending CN117111143A (en) | 2022-05-17 | 2022-05-17 | Method for realizing reliability of structure trap of mountain front punching belt |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117111143A (en) |
-
2022
- 2022-05-17 CN CN202210533424.9A patent/CN117111143A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107526101B (en) | A kind of acquisition and processing method obtaining earthquake reflected wave | |
| CN104698497B (en) | A kind of fine calibration of thin sand-mud interbed layer position | |
| CN103472484B (en) | Horizontal well path optimization method based on RS three-dimensional sensitive earthquake attributive analysis | |
| CN101872024B (en) | Method for carrying out well design by using time-lapse seismic | |
| CN106597533A (en) | Depth domain velocity modeling method for piedmont zone seismic data processing | |
| CN105044770A (en) | Compact glutenite gas reservoir quantificational prediction method | |
| CN105629325B (en) | Foreland basin alluvial fan meticulous depiction and Forecasting Methodology | |
| CN106875471A (en) | Coal measures contains or water barrier Visualization Modeling method | |
| CN105093307B (en) | The oblique true formation thickness acquiring method of Lower Paleozoic strata | |
| CN109738945A (en) | A method of structural map is directly generated using pre-stack depth migration achievement | |
| CN104049283A (en) | Method for identification and prediction of fore-land basin extremely-thick conglomerate body | |
| CN105911592B (en) | A kind of real three dimensional seismic data long wavelength's static correcting method for boring constraint | |
| CN111179415A (en) | Three-dimensional geological model construction method for calcium-bonded rock type uranium ore | |
| Zhao et al. | Relationship between fractures, stress, strike-slip fault and reservoir productivity, China Shunbei oil field, Tarim Basin | |
| CN104375187B (en) | Method for designing cannon well depths based on microlog and non-seismic data | |
| CN117111143A (en) | Method for realizing reliability of structure trap of mountain front punching belt | |
| Deng et al. | Hydrocarbon accumulation conditions and key exploration and development technologies for PL 19–3 oilfield | |
| Jansen | Seismic investigation of wrench faulting and fracturing at Rulison Field, Colorado | |
| Martin et al. | Memoir 71, Chapter 7: Advanced Reservoir Characterization for Improved Oil Recovery in a New Mexico Delaware Basin Project | |
| Martinez et al. | Utility of ground-penetrating radar in near-surface, high-resolution imaging of Lansing-Kansas City (Pennsylvanian) limestone reservoir analogs | |
| Fainstein et al. | Structure and stratigraphy of AVS field, Java sea | |
| Shakir et al. | Prospect generation studies of cretaceous sands, Khipro Area by integrating seismic and well data | |
| Nwankwo et al. | Seismic Refraction Investigation for Groundwater Potential in Parts of Rivers State, Nigeria. | |
| Curkan | Reservoir characterization of channel-belt strata, McMurray Formation, northeastern Alberta | |
| CN109581499A (en) | A method of structural map is generated using anisotropy pre-stack depth migration |
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 |