CN112859159A - Calculation method of goose-row type positive fault sliding belt sliding displacement - Google Patents
Calculation method of goose-row type positive fault sliding belt sliding displacement Download PDFInfo
- Publication number
- CN112859159A CN112859159A CN202110044550.3A CN202110044550A CN112859159A CN 112859159 A CN112859159 A CN 112859159A CN 202110044550 A CN202110044550 A CN 202110044550A CN 112859159 A CN112859159 A CN 112859159A
- Authority
- CN
- China
- Prior art keywords
- sliding
- belt
- fault
- slip
- trend
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 51
- 238000004364 calculation method Methods 0.000 title claims abstract description 15
- 239000010410 layer Substances 0.000 claims abstract description 60
- 238000000034 method Methods 0.000 claims abstract description 30
- 230000001427 coherent effect Effects 0.000 claims abstract description 29
- 230000000694 effects Effects 0.000 claims abstract description 15
- 239000012791 sliding layer Substances 0.000 claims abstract description 14
- 238000005553 drilling Methods 0.000 claims abstract description 6
- 238000010586 diagram Methods 0.000 claims description 31
- 238000011161 development Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 5
- 241000272808 Anser Species 0.000 claims description 3
- 238000007670 refining Methods 0.000 claims description 3
- 241000272814 Anser sp. Species 0.000 claims description 2
- 238000011160 research Methods 0.000 abstract description 5
- 238000012512 characterization method Methods 0.000 abstract description 4
- 238000012795 verification Methods 0.000 abstract description 3
- 230000007547 defect Effects 0.000 abstract description 2
- 239000004576 sand Substances 0.000 description 31
- 230000015572 biosynthetic process Effects 0.000 description 8
- 238000005755 formation reaction Methods 0.000 description 8
- 239000013049 sediment Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 101100400452 Caenorhabditis elegans map-2 gene Proteins 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 241000195474 Sargassum Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 239000013589 supplement Substances 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/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
- 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
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)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The goose-row-type normal-fault sliding belt sliding displacement calculation method is characterized by combining well drilling data and well logging data, carrying out horizon calibration on seismic data, carrying out horizon interpretation to obtain a primary calibration seismic profile, carrying out mutual verification according to the primary calibration seismic profile and an attribute map of an along-layer coherent body to obtain a fine sliding belt plane distribution map and a fine calibration seismic profile of a sedimentary stratum, then combining the fine sliding belt plane distribution map and the fine calibration seismic profile, determining the sliding belt trend, the sliding layer trend and the horizontal offset Li of the sliding layer along the trend in an active period, and finally obtaining the positive-fault sliding belt sliding displacement sigma (Li/sin alpha). Under the condition that the research area lacks obvious two-disc markers or lacks a tail-end broom-shaped fault, quantitative characterization of the glide displacement amount can be carried out on the glide belt, the defects of methods such as obvious fault-breaking markers and the tail-end broom-shaped fault are overcome, and quantitative research on the activity of underground cases of the goose-type glide belt is effectively promoted.
Description
Technical Field
The invention belongs to the technical field of tectonic geological exploration, and particularly relates to a calculation method of a sliding displacement of an goose-rank type normal fault sliding belt.
Background
With the major breakthrough of oil and gas in the sliding zone of the Tarim basin along the north, geologists gradually recognized the importance of the sliding zone inside the basin. Usually, the sliding displacement amount of a sliding fracture zone (sliding zone for short) in the basin is not more than thousands of meters, and the sliding zone is an underground structure pattern with a complex structure. No matter in the early stage or the final stage of the sliding belt movement, the wild goose-row type sliding layer appears in the sliding belt.
The goose-row type positive fault sliding belt refers to a sliding shearing fracture belt consisting of goose-row type positive faults, and particularly relates to an interval range limited by envelope lines of endpoints of all sliding fault layers in a strip-shaped spreading sliding fault system formed under the action of sliding shearing, wherein all the sliding fault layers in a sliding fault group are staggered or inclined in parallel, identical in shape and property, close in length and short in length, and the end parts of adjacent sections of the sliding fault layers are arranged in sequence at the time of pinch-off.
The slip displacement is the slip displacement generated by the slip layer when the movement is generated along the trend of the slip layer. In the basin, the walk-slip displacement of the walk-slip fault layer is an important index for describing the activity intensity and the influence range of the walk-slip fault layer.
The slip fault in this context refers to a fracture plane that is either staggered or skew to the formation.
The following two methods are currently used for the published quantitative characterization of the amount of slip displacement:
the first method is to calculate the displacement of the slip fault by using the slip fault cutting markers (such as river channel, fold, forebody, igneous rock, two-disc markers in outcrop area, etc.). For example: in the underground case of the sedimentary basin, time slices (or depth slices) of three-dimensional seismic data are utilized to restore the river channel outlines on two sides of the fault to original corresponding positions, and further the plane sliding displacement of the sliding fault is calculated. This method has the following limitations: the method needs obvious markers on the two disks of the fault to calculate the sliding displacement, but is influenced by factors such as basin sedimentation environment, sedimentary facies distribution, fault occurrence and the like, and the two disks of the fault have very high possibility of not having the corresponding markers, namely the method does not have the condition for calculating the sliding displacement of the goose-type sliding belt.
The second method is to calculate the sliding distance of the sliding fracture by using the relation between the fracture strike sliding distance and the inclined sliding distance. For example: the broom-shaped fault is usually generated at the tail end of the glide fracture, the trend of the broom-shaped fault is generally nearly perpendicular to the trend of the main fracture, and the glide distance of the glide fracture is quantitatively represented by calculating the horizontal displacement of a plurality of broom-shaped faults at the tail end of the glide fracture. However, the method is suitable for the glide fracture with the tail-end broom-shaped branch fault on a large scale, but for the glide fracture with the small scale in the basin and without the development of the tail-end broom-shaped glide fault, the condition for calculating the glide displacement of the goose-type glide belt by using the method is not provided.
Therefore, the problem that the walk-slip displacement of the walk-slip fault in the basin is small, the seismic data quality precision is low, the relative marker is short and the like is solved, and the technical difficulty which is not overcome in the field is the activity research of the underground case of the goose-row walk-slip zone, particularly the quantitative characterization of the walk-slip displacement.
Disclosure of Invention
The invention provides a calculation method of the walking slip displacement of an goose-row type normal fault walking slip belt, which utilizes well drilling data, well logging data and three-dimensional seismic data to quantitatively calculate the walking slip displacement of the goose-row type normal fault walking slip belt on the basis of seismic interpretation of the goose-row type normal fault walking slip belt.
The technical scheme of the invention is as follows:
a calculation method of the walking slide displacement of an Yanglian type normal fault walking slide belt comprises the following steps:
s1, combining well drilling data and well logging data, performing horizon calibration on each stratum in a sliding zone in the seismic data, performing horizon interpretation on each stratum of two sliding layers in the sliding zone in a seismic profile of the seismic data to obtain a primary calibration seismic profile, and further determining primary profile distribution characteristics of each sliding layer in the sliding zone; the primary section distribution characteristics of the walk-slide layer in the walk-slide zone comprise the trend of the walk-slide layer and the deformation period of main walk-slide activity in the walk-slide zone;
s2, making an attribute map of an in-layer coherent body of a sedimentary stratum in a sliding belt, and based on the attribute map of the in-layer coherent body and in combination with the primary calibration seismic profile map, defining the primary plane distribution characteristics of the sliding belt in the sedimentary stratum and drawing a primary sliding belt plane distribution map of the sedimentary stratum; further refining the development characteristics of the upper and lower disc stratums of each sliding fault in the sliding zone in the primary calibration seismic section map to obtain a fine calibration seismic section map; then, the primary sliding zone plane distribution diagram is adjusted according to the fine calibration seismic profile map in a feedback mode to obtain a fine sliding zone plane distribution diagram, wherein the plane distribution characteristics of the sliding zone on the sedimentary stratum comprise the trend of each sliding zone on the sedimentary stratum and the relative position relation of each sliding zone;
s3, according to the trend of each walk-slip fault in the plane distribution characteristics of the walk-slip zone in the sedimentary stratum, obtaining the horizontal offset L of the walk-slip fault along the trend in the active period in the fine calibration seismic section mapi(ii) a i represents a slip fault index, the trend is a direction perpendicular to the trend of the slip fault;
s4, determining the trend of the sliding belt according to the plane characteristics of the sliding belt in the sedimentary stratum in the fine sliding belt plane distribution map; thenThe total slip displacement of the slip belt during the active period is sigma (L)i/sinαi) (ii) a Wherein alpha is an included angle between the trend of the sliding fault layer and the trend of the sliding belt.
Preferably, in S2, creating an attribute map of the coherent body along the layer of some sedimentary earth formations in the slip zone specifically includes: according to the thickness of some sedimentary strata in the actual active period, time windows with proper widths are respectively selected from the top and the bottom of the top surface of the sedimentary strata for coherent body attribute extraction, and then an attribute map of the coherent body along the stratum of the top surface of the sedimentary strata is obtained.
Preferably, in S2, the combined primary calibration seismic profile includes at least a plurality of primary vertical calibration seismic profiles, and the primary vertical calibration seismic profile is a profile of a strike-slip fault in a primary profile distribution characteristic of the primary calibration seismic profile, which is perpendicular to the strike-slip fault.
Preferably, in S2, the specifying the planar distribution characteristics of the slip band in the sedimentary earth formation specifically includes: re-determining the plane distribution of the slip fault displayed on the attribute map of the interbed coherent body by using the primary vertical calibration seismic profile so as to obtain a primary plane distribution characteristic of the slip zone on the top surface of the sedimentary stratum; the primary plane distribution characteristics of the slip zones on the top surface of the sedimentary earth layer comprise the trend of each slip layer on the top surface of the sedimentary earth layer and the relative position relationship of each slip layer.
Preferably, in S2, further refining development characteristics of the upper and lower tray levels of each walk-slip layer in the walk-slip band in the primary calibration seismic section to obtain a fine calibration seismic section, which specifically includes: and performing feedback adjustment and fine drawing on the development characteristics of the stratums of the upper and lower discs of the slip fault of at least the primary vertical calibration seismic section in the primary calibration section through the slip zone plane distribution diagram to obtain a fine calibration seismic section.
Preferably, in S3, the horizontal offset Li of the slip layer along the dip during the active period is the horizontal offset of the top surface of the sedimentary earth formation between the upper and lower plates of the slip layer during the active period.
Preferably, in S4, the determining the running and sliding belt trend specifically includes: and on the fine sliding belt plane distribution diagram, applying a trial and error method, and when the trend of each sliding layer in the sliding belt is consistent with the included angle of a certain direction, the direction is the trend of the sliding belt.
Preferably, in S4, the total slip displacement amount of the slip belt in the active period is Σ (Li/sin α), and specifically includes: in the fine sliding belt plane distribution diagram, determining the sliding component Si ═ L of the sliding belt on the sliding layer according to the included angle alpha between the sliding layer trend and the sliding belt trendi/sinαiThen the total sliding displacement of the sliding belt in the active period is Σ (Si).
Preferably, in S2, the certain sedimentary earth formations include sedimentary earth formations in which deformation periods of main sliding activities in the sliding belt exist, and in the attribute map of the coherent body along the layer, the part with a large coherence coefficient indicates a fracture.
Preferably, in the step S2, when a time window with a proper width is selected from top and bottom of the top surface of some sedimentary formations for coherent body attribute extraction, the width of the time window with the proper width is 1/20-1/5 times the thickness of the sedimentary formations.
Compared with the prior art, the invention has the advantages that: the goose-row type positive fault sliding belt sliding displacement calculation method is characterized by combining well drilling data and well logging data, carrying out horizon calibration on seismic data, carrying out horizon interpretation to obtain a primary calibration seismic profile, mutually verifying the primary calibration seismic profile and a along-layer coherent body attribute map to obtain a fine sliding belt plane distribution map and a fine calibration seismic profile of a sedimentary stratum, then combining the fine sliding belt plane distribution map and the fine calibration seismic profile to determine the sliding belt trend, the sliding layer trend and the horizontal offset Li of the sliding layer along the trend in the active period, and finally obtaining positive fault sliding belt sliding displacement sigma (Li/sin alpha). Under the condition that the research area lacks obvious two-disc markers or lacks a tail-end broom-shaped fault, the quantitative characterization of the sliding displacement of the sliding belt can be carried out, and further data support is provided for the moving characteristics and the influence range of the sliding belt. The method makes up the defects of obvious methods of fault-breaking markers and tail-end broomstick faults and the like, supplements the calculation of the running slip displacement quantity of the goose-type running-sliding belt, and effectively promotes the quantitative research of the activity of underground cases of the goose-type running-sliding belt.
Drawings
FIG. 1 is a flow chart of a method for calculating the amount of slide displacement of an orthostictus normal fault slide belt according to the invention;
FIG. 2 is a primary calibration seismic section diagram of stratum calibration and stratum interpretation of the goose-column type gliding belt of the Wang Home in the three-dimensional seismic data coverage area of the east-Yingzhuang depression according to the calculation method of the goose-column type normal-fault gliding belt gliding displacement amount of the invention;
FIG. 3 is an attribute diagram of an along-layer coherence of the top surface of a Shake river street group sand 2 segment in a goose-row type gliding zone of a Wang Home guard in a covering area of a three-dimensional earthquake data with a Yang depression east cave according to the calculation method of the goose-row type normal-fault gliding zone gliding displacement amount of the invention, which is a geodetic coordinate system;
FIG. 4 is a walking and sliding belt plane distribution diagram of the top surface of the sand in the goose-row walking and sliding belt group of sand 2 segment of the sand in the goose-row walking and sliding belt of the king house post in the three-dimensional seismic data coverage area with the depression of yang down in east China according to the calculation method of the walking and sliding displacement of the goose-row type positive-fault walking and sliding belt of the invention, which is a geodetic coordinate system;
FIG. 5 is a finely calibrated seismic section diagram of stratum calibration and stratum interpretation of the goose-row type gliding belt of the Wang Home in the three-dimensional seismic data coverage area of the east-Yingzhuang depression according to the calculation method of the glide displacement amount of the goose-row type normal-fault gliding belt of the invention;
FIG. 6 is a schematic diagram of the slip displacement of the goose-row type normal-fault slip belt in the calculation method of the slip displacement of the goose-row type normal-fault slip belt.
Detailed Description
To facilitate an understanding of the invention, the invention is described in more detail below with reference to figures 1-6 and the specific examples.
Fig. 1 is a flowchart of a method for calculating a sliding displacement amount of an anser positive fault sliding belt according to the present invention, and the following embodiment is an analysis and calculation process of a king home-watch positive fault sliding belt in a three-dimensional depression-east earthquake data coverage area according to the method for calculating the sliding displacement amount of the anser positive fault sliding belt according to the present invention. The method for calculating the walking and sliding displacement of the Wangjiagang positive fault walking and sliding belt in the east-sagged Dongying three-dimensional seismic data coverage area of the Jiyang sag comprises the following steps of:
s1, combining the well drilling data and the well logging data, calibrating the horizon of each stratum in the sliding zone in the seismic data, and interpreting the horizons of each two stratums of each sliding zone in the seismic profile of the seismic data, especially interpreting each stratum of the upper and lower plates of the sliding zone to obtain a primary calibration seismic profile as shown in figure 2, wherein in figure 2, T respectively represents the standard reflection interfaces of the tops (or bottoms) of different stratums, namely, calibrating the horizons of each stratum, wherein, the horizon calibration and the horizon interpretation of each horizon are as follows:
T0a bottom surface of the Minghua ballast assembly,
T1a bottom surface of a ceramic group of a library,
T2' bottom of the Dongying group,
T2the bottom surface of the sand 1 section of the sand river street group (the top surface of the sand 2 section of the sand river street group),
T3the bottom surface of the sand 2 section of the sand river street group,
T4the sand 3 is arranged on the bottom surface of the upper section of the sand,
T6a lower bottom surface of the sand 3 of the sand river street group,
T7a sand section 4 bottom surface of the sand river street group,
tr is the bottom surface of the hole shop group,
according to the primary calibration seismic profile map 2, determining primary profile distribution characteristics of all walk-slide faults in a walk-slide zone, wherein the primary profile distribution characteristics of the walk-slide faults comprise the trend of the walk-slide faults; in fig. 2, the horizontal broken line is different strata and their deformation, and the thick solid line close to the vertical direction is a rough marked slip fault. Determining the deformation period of main sliding activities in the goose rank belt as the sediment periods of sand 3 and sand 2 of the sand group of the sand river street according to the cutting relation between the roughly marked sliding fault layer and the stratum;
s2, making an attribute map of the along-layer coherent body, wherein the attribute map of the along-layer coherent body can be made for each stratum in the walking sliding belt, and the attribute map of the along-layer coherent body is analyzed; and (3) making an attribute map of the along-layer coherent body of the stratum in the main sliding activity deformation period in the Yanglian zone, and analyzing the attribute map of the along-layer coherent body.
In the attribute diagram of the along-layer coherent body, a part with a large coherence coefficient is a part with a fracture, and a solid line part shown in fig. 3 is a part where a slip fault of the attribute analysis of the along-layer coherent body is located.
Fig. 3 is an attribute diagram of the along-layer coherent body of the top surface of the sand 2 segment of the sand group in the sliding zone, which is obtained by selecting time windows with a width of 20ms from the top to the bottom of the top surface of the sand 2 segment of the sand group for coherent body attribute extraction according to the actual activity deformation period and the thickness of the deposited stratum of the sand 2 segment of the sand group of 200 ms.
Combining a profile which vertically cuts the trend of each slip fault in the primary calibration seismic profile, namely the primary vertical calibration seismic profile, re-defining the distribution and the trend of the slip faults in the slip fault distribution range displayed on the along-layer coherent body attribute map, defining the primary plane distribution characteristics of the slip zones in the sedimentary stratum and drawing a primary slip zone plane distribution map of the sedimentary stratum; in the re-determination process, the primary vertical calibration seismic profile and the attribute map of the coherent body along the stratum need to be used, and mutual verification is repeated and carried out for many times, so that the primary plane distribution characteristic of the sliding belt on the top surface of the sedimentary stratum is obtained; the planar distribution characteristics of the slip zones on the top surface of the sedimentary earth layer comprise the trend of each slip layer on the top surface of the sedimentary earth layer and the relative position relationship of each slip layer. And drawing a primary sliding belt plane distribution diagram of the sedimentary stratum according to the primary plane distribution characteristics of the sliding belt on the top surface of the sedimentary stratum, and when mutual verification is carried out for multiple times, performing feedback adjustment on the upper and lower disc level development characteristics of each sliding layer in the sliding belt in the primary calibration seismic section diagram to obtain a fine depiction according to the primary sliding belt plane distribution diagram, so as to obtain the fine calibration seismic section diagram shown in the figure 5. And then adjusting the primary sliding zone plane distribution diagram according to the fine calibration seismic profile to obtain a fine sliding zone plane distribution diagram as shown in fig. 4, wherein fig. 4 is a fine sliding zone plane distribution diagram of the top surface of the sand 2 segment of the intra-sliding-zone sargassum group repeatedly verified and drawn according to fig. 2, fig. 3 and fig. 5. In the fine sliding strip plane distribution diagram, each realization in the fine sliding strip plane distribution diagram is a sliding fault which is mutually verified and drawn, and the direction of a short line perpendicular to the trend of the fine sliding strip plane distribution diagram is the inclination of the fine sliding strip plane distribution diagram, wherein the short line is above the sliding fault and indicates that the sliding fault is a north inclined fault, and the short line is below the sliding fault and indicates that the sliding fault is a south inclined fault.
S3, according to the trend of each walk-slip fault, obtaining the horizontal offset of the top surface of the sand 2-segment sedimentary earth layer of the sand street group in the active period between the upper plate and the lower plate of the walk-slip fault, namely the horizontal offset L along the inclination in the fine calibration seismic section diagram 5i(ii) a Where i denotes a slip fault index, as shown in fig. 4, there are 32 slip faults, of which 18 are north dip faults and 14 are south dip faults, so i is 1,2,3 … …, 32.
The tendency is a direction perpendicular to the strike of the slip fault;
s4, by means of the relative position relationship between the trend of the slip zone in the fine slip zone plane distribution diagram and the slip zone, when the included angle between the trend of each slip zone in the slip zone and a certain direction tends to be consistent, the direction is the trend of the slip zone by using a trial and error method, and is shown as a dotted line in fig. 4.
According to the trend of the sliding belt and the trend of the sliding fault layer in the graph 4, the included angle between the trend of the sliding belt and the trend of the sliding fault layer is determined to be alpha. The total slip displacement of the slip zone on the top surface of the sediment layer of the sand 2 section of the sand river street group in the activity period is sigma (S)i)=∑(Li/sinαi)。
Wherein, the formula S of the sliding component of the sliding belt at the sliding fault layer ii=Li/sinαiThe principle is shown in fig. 6, in the goose-row type sliding fault model of the sliding belt, the direction shown by the dotted line is the trend of the sliding belt and is along the sliding beltThe left side is a footwall (as indicated by an arrow 1 in fig. 6), the right side is an upper wall (as indicated by an arrow 2 in fig. 6), the shaded part is a slip layer, the thick solid lines of the upper boundary of the shaded part are a slip layer footwall cutting line and an upper wall cutting line respectively, the small arrow of the shaded part is a horizontal offset L of the obtained slip layer i along the inclination in the active period in the fine calibration seismic profilei(ii) a Therefore, the formula of the sliding component of the sliding belt at the sliding fault layer i is Si=LiA/sin alpha; the total sliding displacement of the sliding belt in the activity period on the top surface of the sediment layer of the sand 2 section of the sand river street group is sigma (S)i)=∑(Li/sinαi)。
Table 1 shows the amount of slip on the top surface of the sedimentary earth layer at sand 2 in the group of the sand roads, which is statistically obtained according to fig. 3.
It should be noted that the above-described embodiments may enable those skilled in the art to more fully understand the present invention, but do not limit the present invention in any way. Therefore, although the present invention has been described in detail with reference to the drawings and examples, it will be understood by those skilled in the art that various changes and modifications can be made therein without departing from the spirit and scope of the invention.
Claims (10)
1. A calculation method of the walking slide displacement of an goose-rank type normal fault walking slide belt is characterized by comprising the following steps:
s1, combining well drilling data and well logging data, performing horizon calibration on each stratum in a sliding zone in the seismic data, performing horizon interpretation on each stratum of two sliding layers in the sliding zone in a seismic profile of the seismic data to obtain a primary calibration seismic profile, and further determining primary profile distribution characteristics of each sliding layer in the sliding zone; the primary section distribution characteristics of the walk-slide layer in the walk-slide zone comprise the trend of the walk-slide layer and the deformation period of main walk-slide activity in the walk-slide zone;
s2, making an attribute map of an in-layer coherent body of a sedimentary stratum in a sliding belt, and based on the attribute map of the in-layer coherent body and in combination with the primary calibration seismic profile map, defining the primary plane distribution characteristics of the sliding belt in the sedimentary stratum and drawing a primary sliding belt plane distribution map of the sedimentary stratum; further refining the development characteristics of the upper and lower disc stratums of each sliding fault in the sliding zone in the primary calibration seismic section map to obtain a fine calibration seismic section map; then, the primary sliding zone plane distribution diagram is adjusted according to the fine calibration seismic profile map in a feedback mode to obtain a fine sliding zone plane distribution diagram, wherein the plane distribution characteristics of the sliding zone on the sedimentary stratum comprise the trend of each sliding zone on the sedimentary stratum and the relative position relation of each sliding zone;
s3, according to the trend of each walk-slip fault in the plane distribution characteristics of the deposit stratum, obtaining the horizontal offset Li of the walk-slip fault along the trend in the active period in the fine calibration seismic profile; i represents a slip fault index, the trend is a direction perpendicular to the trend of the slip fault;
s4, determining the trend of the sliding belt according to the plane characteristics of the sliding belt in the sedimentary stratum in the fine sliding belt plane distribution map; the total slip displacement amount of the slip belt in the active period is sigma (Li/sin alpha i); wherein alpha is an included angle between the trend of the sliding fault layer and the trend of the sliding belt.
2. The method for calculating the walking-slip displacement of the anser-type normal fault walking-slip belt according to claim 1, wherein in the step S2, the making of the along-layer coherent body attribute map of some sedimentary strata in the walking-slip belt specifically comprises: according to the thickness of some sedimentary strata in the actual active period, time windows with proper widths are respectively selected from the top and the bottom of the top surface of the sedimentary strata for coherent body attribute extraction, and then an attribute map of the coherent body along the stratum of the top surface of the sedimentary strata is obtained.
3. The method of claim 2, wherein the integrated primary calibration seismic profile of S2 at least comprises a plurality of primary vertical calibration seismic profiles, and the primary vertical calibration seismic profiles are profiles of the strike-slip fault in the primary profile distribution feature perpendicular to the strike-slip fault in the primary calibration seismic profiles.
4. The method for calculating the walking-slip displacement of the anser-type normal fault walking-slip belt according to claim 3, wherein in the step S2, the step S of defining the planar distribution characteristics of the walking-slip belt in the sedimentary stratum specifically comprises the following steps: re-determining the plane distribution of the slip fault displayed on the attribute map of the interbed coherent body by using the primary vertical calibration seismic profile so as to obtain the plane distribution characteristic of the slip zone on the top surface of the sedimentary stratum; the planar distribution characteristics of the slip zones on the top surface of the sedimentary earth layer comprise the trend of each slip layer on the top surface of the sedimentary earth layer and the relative position relationship of each slip layer.
5. The method for calculating the walking-sliding displacement amount of the anser-type normal fault walking-sliding belt according to claim 4, wherein in the step S2, the development characteristics of the upper and lower plate levels of each walking-sliding belt in the primary calibration seismic profile are further refined to obtain a fine calibration seismic profile, and the method specifically comprises the following steps: and performing feedback adjustment and fine drawing on the development characteristics of the stratums of the upper and lower discs of the slip fault of at least the primary vertical calibration seismic section in the primary calibration section through the slip zone plane distribution diagram to obtain a fine calibration seismic section.
6. The method for calculating the amount of the aquaplaning offset of the orthonormal fault aquaplaning zone according to claim 4 or 5, wherein in the step S3, the horizontal offset Li of the aquaplaning zone along the inclination during the active period is the horizontal offset of the top surface of the sedimentary earth strata during the active period between the upper plate and the lower plate of the aquaplaning zone.
7. The method for calculating the walking displacement amount of the anser type normal fault walking belt according to claim 4 or 5, wherein in the step S4, the step of determining the walking belt trend specifically comprises the steps of: and on the fine sliding belt plane distribution diagram, applying a trial and error method, and when the trend of each sliding layer in the sliding belt is consistent with the included angle of a certain direction, the direction is the trend of the sliding belt.
8. The method for calculating the aquatype normal slice gliding band gliding displacement amount of the goose according to claim 7, wherein in S4, the total gliding displacement amount of the gliding band during the activity period is Σ (Li/sin α), which specifically comprises: and in the fine sliding belt plane distribution diagram, determining a sliding component Si of the sliding belt in the sliding belt according to an included angle alpha between the trend of each sliding belt and the trend of the sliding belt, wherein the sliding component Si is Li/sin alpha i, and the total sliding displacement of the sliding belt in the activity period is sigma (Si).
9. The method for calculating the aquaplaning offset of an orthonormal fault aquaplaning zone according to claim 1, wherein in S2, the sedimentary strata includes sedimentary strata in which a deformation period of main aquaplaning activity is located in the aquaplaning zone, and a fracture is indicated in the attribute diagram of the coherent body along the layer where the portion with a large coherence coefficient is.
10. The method for calculating the walking displacement of the goose-type normal fault walking belt according to claim 1, wherein in S2, when a time window with a proper width is selected from top and bottom of the top surface of some sedimentary strata respectively for coherent body attribute extraction, the width of the time window with the proper width is 1/20-1/5 times the thickness of the sedimentary strata.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110044550.3A CN112859159B (en) | 2021-01-13 | 2021-01-13 | Calculation method of goose-row type positive fault sliding belt sliding displacement |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110044550.3A CN112859159B (en) | 2021-01-13 | 2021-01-13 | Calculation method of goose-row type positive fault sliding belt sliding displacement |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112859159A true CN112859159A (en) | 2021-05-28 |
CN112859159B CN112859159B (en) | 2022-03-11 |
Family
ID=76003517
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110044550.3A Expired - Fee Related CN112859159B (en) | 2021-01-13 | 2021-01-13 | Calculation method of goose-row type positive fault sliding belt sliding displacement |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112859159B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114117827A (en) * | 2022-01-24 | 2022-03-01 | 四川省公路规划勘察设计研究院有限公司 | Random simulation and parameter sensitivity analysis method for Yanglian fault model |
CN114545496A (en) * | 2022-02-25 | 2022-05-27 | 中国地质大学(北京) | Method for judging sliding direction of underground sliding fracture |
CN115809397A (en) * | 2023-02-06 | 2023-03-17 | 西南交通大学 | Method and device for calculating dislocation stratum deformation of movable fracture zone and readable storage medium |
CN114647004B (en) * | 2022-02-25 | 2024-03-15 | 中国地质大学(北京) | Method for confirming sliding directions of underground sliding fracture at different periods |
CN117805902A (en) * | 2024-02-28 | 2024-04-02 | 山东石油化工学院 | Method and system for quantitatively representing relation between basin stretching and sliding action strength |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001090783A1 (en) * | 2000-05-25 | 2001-11-29 | Petrecon Australia Pty Ltd | Method for detecting direction and relative magnitude of maximum horizontal stress in earth's crust |
CN106556869A (en) * | 2017-01-19 | 2017-04-05 | 中国石油大学(华东) | Method that tomography walk sliding displacement is portrayed under a kind of extention environment quantitatively |
KR101859650B1 (en) * | 2017-08-21 | 2018-05-18 | 충북대학교 산학협력단 | A procedure of fault interpretation from 3D seismic data using displacement patterns of normal faults |
CN108106582A (en) * | 2017-12-14 | 2018-06-01 | 中国石油化工股份有限公司 | Sliding displacement computational methods are walked based on Strike-slip faulted under multi-stage tectonic movement background |
CN108107476A (en) * | 2016-11-25 | 2018-06-01 | 中国石油化工股份有限公司 | Broom shape strike-slip fault walks the computational methods of sliding displacement |
US20180321404A1 (en) * | 2015-11-05 | 2018-11-08 | Schlumberger Technology Corporation | Generation of fault displacement vector and/or fault damage zone in subsurface formation using stratigraphic function |
CN109856678A (en) * | 2019-03-07 | 2019-06-07 | 中国石油化工股份有限公司 | A kind of calculation method of the pressure property turned round strike-slip fault Strike slip rate |
-
2021
- 2021-01-13 CN CN202110044550.3A patent/CN112859159B/en not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001090783A1 (en) * | 2000-05-25 | 2001-11-29 | Petrecon Australia Pty Ltd | Method for detecting direction and relative magnitude of maximum horizontal stress in earth's crust |
US20180321404A1 (en) * | 2015-11-05 | 2018-11-08 | Schlumberger Technology Corporation | Generation of fault displacement vector and/or fault damage zone in subsurface formation using stratigraphic function |
CN108107476A (en) * | 2016-11-25 | 2018-06-01 | 中国石油化工股份有限公司 | Broom shape strike-slip fault walks the computational methods of sliding displacement |
CN106556869A (en) * | 2017-01-19 | 2017-04-05 | 中国石油大学(华东) | Method that tomography walk sliding displacement is portrayed under a kind of extention environment quantitatively |
KR101859650B1 (en) * | 2017-08-21 | 2018-05-18 | 충북대학교 산학협력단 | A procedure of fault interpretation from 3D seismic data using displacement patterns of normal faults |
CN108106582A (en) * | 2017-12-14 | 2018-06-01 | 中国石油化工股份有限公司 | Sliding displacement computational methods are walked based on Strike-slip faulted under multi-stage tectonic movement background |
CN109856678A (en) * | 2019-03-07 | 2019-06-07 | 中国石油化工股份有限公司 | A kind of calculation method of the pressure property turned round strike-slip fault Strike slip rate |
Non-Patent Citations (3)
Title |
---|
周卿等: "走滑断裂和裂缝发育带的地震地质综合识别", 《物探化探计算技术》 * |
房璐等: "辽东湾坳陷东部地区走滑断裂发育特征及走滑量的估算", 《地质力学学报》 * |
李伟等: "辽东湾坳陷新生代主干断裂走滑量的估算及其地质意义", 《大地构造与成矿学》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114117827A (en) * | 2022-01-24 | 2022-03-01 | 四川省公路规划勘察设计研究院有限公司 | Random simulation and parameter sensitivity analysis method for Yanglian fault model |
CN114117827B (en) * | 2022-01-24 | 2022-04-15 | 四川省公路规划勘察设计研究院有限公司 | Random simulation and parameter sensitivity analysis method for Yanglian fault model |
CN114545496A (en) * | 2022-02-25 | 2022-05-27 | 中国地质大学(北京) | Method for judging sliding direction of underground sliding fracture |
CN114545496B (en) * | 2022-02-25 | 2024-01-19 | 中国地质大学(北京) | Method for judging sliding direction of underground sliding fracture |
CN114647004B (en) * | 2022-02-25 | 2024-03-15 | 中国地质大学(北京) | Method for confirming sliding directions of underground sliding fracture at different periods |
CN115809397A (en) * | 2023-02-06 | 2023-03-17 | 西南交通大学 | Method and device for calculating dislocation stratum deformation of movable fracture zone and readable storage medium |
CN115809397B (en) * | 2023-02-06 | 2023-04-18 | 西南交通大学 | Method and equipment for calculating dislocation stratum deformation of movable fracture zone and readable storage medium |
CN117805902A (en) * | 2024-02-28 | 2024-04-02 | 山东石油化工学院 | Method and system for quantitatively representing relation between basin stretching and sliding action strength |
CN117805902B (en) * | 2024-02-28 | 2024-05-24 | 山东石油化工学院 | Method and system for quantitatively representing relation between basin stretching and sliding action strength |
Also Published As
Publication number | Publication date |
---|---|
CN112859159B (en) | 2022-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112859159B (en) | Calculation method of goose-row type positive fault sliding belt sliding displacement | |
Tian et al. | Multi-layered Ordovician paleokarst reservoir detection and spatial delineation: A case study in the Tahe Oilfield, Tarim Basin, Western China | |
CN104047597B (en) | Oily mud shale stratum log standardized method | |
CN110824557B (en) | Quantitative karst ancient landform recovery method based on well-seismic combination | |
Sullivan et al. | An integrated approach to characterization and modeling of deep-water reservoirs, Diana field, western Gulf of Mexico | |
CN114114459B (en) | Deep-ultra-deep carbonate rock thin reservoir prediction method under phase control constraint | |
CN111624677B (en) | Compact sandstone gas reservoir single-phase river channel sand body boundary depicting method | |
CN112505754B (en) | Method for collaborative partitioning sedimentary microfacies by well-seismic based on high-precision sequence grid model | |
CN116856845A (en) | Tight sandstone gas reservoir horizontal well geosteering method | |
CN105089659A (en) | Conglomerate oil reservoir seepage unit recognizing method | |
CN112180463A (en) | Prediction method for sandstone spread in transition layer section | |
CN115685327A (en) | Method for predicting distribution range of deep-water turbid accumulated sand body | |
CN104345337B (en) | A kind of time control reservoir parameter modeling method for seismic inversion | |
Martinez et al. | Utility of ground-penetrating radar in near-surface, high-resolution imaging of Lansing-Kansas City (Pennsylvanian) limestone reservoir analogs | |
Nielsen et al. | Pilot phase of the Aguada Federal Block, black-oil window | |
Duvbiama et al. | 3D Static Modelling of an offshore field in the Niger-delta | |
Kondratyev et al. | Submarine fan reservoir architecture and heterogeneity influence on hard-to-recover reserves. Achimov Fm | |
CN113031065B (en) | Oilfield sand geological modeling method | |
Rocco | 2. Comparative Geological and Geophysical Study of the Po Basin (Italy) | |
Podnebesnykh et al. | New approach to the evaluation of the structure of initial reserves in Ozhginskoe gas-oil field | |
ADEBOYE | SEISMIC RESERVOIR CHARACTERIZATION AND 3D MODELLING OF OY FIELD, OFFSHORE NIGER DELTA | |
Shrivastava et al. | Delineation of structural features by gravity and magnetic surveys in Than-Somasar tract, Saurashtra peninsula, Gujarat (India)-Implications on uranium mineralization | |
Jabir | High-resolution 3D geocellular facies and property modeling of bentiu-1 reservoir interval, Jack Oilfield, Fula Sub-basin, Sudan | |
Langenheim et al. | Thickness and geometry of Cenozoic deposits in California Wash area, Nevada, based on gravity and seismic-reflection data | |
CN114687725A (en) | Identification method for cardiac beach scale |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220311 |