CN117872472B - A method for determining the irregular cross-sectional shape of normal faults in sandy mud strata - Google Patents
A method for determining the irregular cross-sectional shape of normal faults in sandy mud strata Download PDFInfo
- Publication number
- CN117872472B CN117872472B CN202311730804.2A CN202311730804A CN117872472B CN 117872472 B CN117872472 B CN 117872472B CN 202311730804 A CN202311730804 A CN 202311730804A CN 117872472 B CN117872472 B CN 117872472B
- Authority
- CN
- China
- Prior art keywords
- seismic
- fault
- target area
- determining
- area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000001788 irregular Effects 0.000 title claims abstract description 22
- 230000005540 biological transmission Effects 0.000 claims abstract description 13
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 8
- 238000012800 visualization Methods 0.000 claims abstract 2
- 238000011156 evaluation Methods 0.000 claims description 37
- 239000010410 layer Substances 0.000 claims description 23
- 238000004321 preservation Methods 0.000 claims description 20
- 230000002159 abnormal effect Effects 0.000 claims description 8
- 238000005553 drilling Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 239000011229 interlayer Substances 0.000 claims description 3
- 238000012512 characterization method Methods 0.000 abstract 1
- 238000011160 research Methods 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 7
- 239000004576 sand Substances 0.000 description 7
- 239000013049 sediment Substances 0.000 description 5
- 238000010276 construction Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000013517 stratification Methods 0.000 description 4
- 238000013519 translation Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000013508 migration Methods 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- 230000007480 spreading Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000005137 deposition process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
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 present invention relates to a method for determining the irregular cross-sectional shape of a normal fault in a sandy mud formation, comprising: obtaining a seismic profile based on seismic data in a target area, evaluating the seismic profile, thereby determining the reliability of the seismic data; determining the location of the seismic zone in the target area through seismic data that meets the reliability requirements, and performing fault interpretation on the seismic data to determine the fault direction in the target area; determining the classification and characteristics of the transmission zone in the target area based on the fault direction, and establishing a three-dimensional visualization model based on the classification and characteristics of the transmission zone to implement the spatial distribution morphology of the fault in the area; gridding the target area, optimizing the seismic profile based on the positional relationship between each grid area and the seismic zone combined with the spatial distribution morphology of the fault, and providing a basis for well site deployment. The present invention solves the problem of fine characterization of the spatial distribution of faults in complex fault zones, and provides a basis for well site deployment.
Description
Technical Field
The invention relates to the technical field of seismic interpretation, in particular to a method for determining the irregular cross-section shape of a positive fault of a sand mud stratum.
Background
Growth faults refer to faults that develop during the deposition process, typically at the edges of the deposition basin. The most typical growth fault develops in the delta basin, when a weak area appears in the silty sediment environment at the front edge of the delta, the weak area is broken due to construction movement, when the river carries silty sand to pass through the weak area, a fault cliff at the front edge of the delta is filled up, then the fault moves again to form a new fault cliff, the delta advances forward, the fault cliff is filled up by the silty sediment, the fault moves continuously, the sediment is filled continuously, the process is repeated, finally, the growth fault is formed, the section shape of the growth fault is irregular, the space spreading form is complex, and the construction interpretation difficulty is high.
The current fault interpretation method is mainly based on the mode that an earthquake section vertical to fault trend and attribute slices such as coherence, ant body and the like are combined in a horizontal section mode, specifically, fault dip angles and fault distances are confirmed by using section earthquake interpretation results, fault plane spreading and combination relations are confirmed by using the attribute slices, and finally, fault space spreading is realized, so that the purpose of describing faults is achieved. The method has good effect on conventional fault interpretation, but is difficult to finely describe faults with irregular section shape and complex space distribution form of sand and mud stratum development.
Chinese patent publication No.: CN108411844a discloses a method and a device for analyzing a mud-rock flow velocity field of a natural channel irregular section, by adopting a gridding idea, the natural channel irregular section is subjected to gridding division, the uneven distribution characteristics of the velocity in the channel section direction and the mud depth direction can be considered in mud-rock flow velocity calculation, and a velocity field model can be obtained by a practical method. Meanwhile, the result is checked by taking the flow calculated by the rain and flood method suggested in the debris flow control standard as a standard, so that the subjectivity and uncertainty of the selection of the roughness coefficient in the calculation of the flow velocity can be overcome.
Therefore, the current method for determining the irregular cross-sectional shape cannot solve the problems that the fine depiction of the fault space distribution of the complex fracture zone cannot well describe the migration evolution process of the fault block high points of different layers, thereby providing a basis for well position deployment.
Disclosure of Invention
Therefore, the invention provides a method for determining the irregular cross-section shape of a positive fault of a sand mud stratum, which is used for solving the problems that in the prior art, the precise description of the fault space spread of a complex fracture zone cannot be well described, and the migration evolution process of high points of fault blocks of different layers cannot be well described, so that basis is provided for well position deployment.
In order to achieve the above object, the present invention provides a method for determining the shape of an irregular cross section of a normal fault of a sand-mud stratum, comprising,
Step S1, obtaining a seismic section according to seismic data of a target area, and evaluating the seismic section so as to determine the reliability of the seismic data;
S2, determining the position of a seismic zone in a target area through seismic data with reliability reaching the requirement, and performing fault interpretation on the seismic data so as to determine the fault trend of the target area;
Step S3, determining the classification and the characteristics of a target area conveyor belt according to the fault trend, and establishing a three-dimensional visual model according to the classification and the characteristics of the conveyor belt so as to realize the fault space distribution form of the area;
And S4, grid division is carried out on the target area, the seismic section is optimized according to the position relation between each grid area and the seismic zone and the combination of fault space distribution form, and a basis is provided for well position deployment.
Further, the evaluation criteria for evaluating the seismic section include the amplitude preservation of the seismic section, the authenticity of the structural form of the seismic section, and the accuracy of fault homing of the seismic section;
and determining the reliability grade of the seismic data according to the amplitude conservation property, the construction morphology and the number of items meeting the single determination condition in the fault homing of the seismic section.
Further, calculating a first difference absolute value according to the reflection time of the in-phase axis where the target layer is located and the corresponding geological stratification conversion time,
Comparing the first difference absolute value with a set structural form authenticity evaluation value, determining whether the authenticity of the structural form of the seismic section meets the requirement, and if so, determining that the seismic section meets the single determination condition.
Further, according to the historical data obtained by collecting the well logging and coring data, lithology and field outcrop data of the research area, determining the stratum structure and stratum lithology characteristics, and adopting the synthetic record layer calibration at the well point to determine the amplitude preservation of the seismic section,
If the consistency coefficient of the synthetic record and the well side seismic channel is greater than or equal to 80%, judging that the amplitude preservation of the seismic section meets the requirement, and determining that the amplitude preservation of the seismic section meets the single determination condition;
And comparing the actual drilling breakpoint position with the seismic data breakpoint position to determine whether the fault homing accuracy of the seismic section meets the requirement, and if so, determining that the fault homing accuracy meets the single determination condition.
Further, in the step S2, the method for performing fault interpretation on the seismic data meeting the requirement of reliability comprises the following steps,
Step S21, determining the arrangement of a broken area on the seismic section and analyzing the size of a vertical broken distance on the seismic section;
Step S22, connecting the break points of the upper and lower seismic reflection layer large break distance faults with the break points of the middle layer small break distance faults in the longitudinal direction to form an irregular section of the same fault;
And S23, realizing the relationship of the arrangement, the falling direction and the plane spread of the fracture area of the interlayer translation fault by identifying the translation fault, and performing the section closure interpretation by utilizing the longitudinal and transverse section closure.
Further, carrying out grid division on the target area to obtain grid areas;
And setting seismic magnitude evaluation values for each grid region, and determining whether the target region meets the stratum stability requirement of well position deployment according to the acquired historical magnitude values of each grid region and the seismic magnitude evaluation values.
Further, for any grid area, performing single-point stability comparison according to the second difference absolute value and the set single-point stability evaluation value, and determining whether the stratum stability in the grid area meets the stratum requirement of well position deployment;
And the second absolute value of the difference is the absolute value of the difference between the historical magnitude value and the seismic magnitude evaluation value of the grid area.
Further, renumbering each grid region conforming to well site deployment and calculating a conforming score;
renumbering each grid area which does not accord with well position deployment and calculating a non-accord score;
and calculating the overall stability evaluation value of the target area according to the obtained coincidence scores and the non-coincidence scores, and comparing the overall stability evaluation value with a set overall stability standard evaluation value to determine whether the target area needs to be replaced for well position deployment.
Further, splicing the grid areas which do not meet the stratum requirements of well position deployment to form a plurality of seismic zones, for any seismic zone, sequentially rounding with 1 time of division length, 2 times of division length and 3 times of division length as radius by taking the centroid of the seismic zone area as a round point, and determining seismic zones with different grades in the seismic zone;
And carrying out three equal division on the distance between the farthest edge point of the seismic zone and the circular center to obtain the division length.
Further, overlapping positions are formed in the overlapping areas of the seismic zone ranges;
And determining the safety state grade or the abnormal state grade of the coincident position according to the number and grade of the seismic bands forming the coincident position.
Compared with the prior art, the invention has the beneficial effects that the invention adopts the irregular section interpretation technology based on mechanical analysis to interpret the section based on the geological stress, the structural evolution and the sediment filling analysis as the guidance, truly describes the fracture structural characteristics of the underground, particularly utilizes the analysis and test data of the drilled well logging, coring and stratum and the open-air outcrop data to establish the combination relation of the stratum grillwork and lithology of the research area, analyzes the structural development characteristics and the formation mechanism of the research area, establishes the fault real interpretation mode, provides the analysis and judgment technical method for the reliability of the seismic data of the research area, synthesizes the record layer calibration and judgment of the conservation of the section through the well point, the method comprises the steps of judging the authenticity of a structural form of an earthquake section through well connection explanation of a well drilling geological stratification on the earthquake section, judging the accuracy of fault homing of the earthquake section according to comparison of a breakpoint encountered by well drilling and a breakpoint position of earthquake data, under the condition that the earthquake data is reliable, guiding by using a geological mode, optimizing the earthquake explanation section, and carrying out fine explanation of a complex fracture zone combining points, lines, planes and the body by combining classification and position confirmation of a transmission zone on the basis of closed calibration of a single well and Lian Jing, so that the combination relation of the complex fracture zone faults is solved, more importantly, the migration evolution process of high points of fracture blocks of different layers is well described, and the explanation result accords with the movement rule and the deposition characteristics of the geological structure, thereby providing a basis for well position deployment.
In particular, by carrying out horizon grid interpretation on a target layer in seismic data, determining the time range of the target layer, obtaining the structural trend of the stratum, slicing along the structural trend surface of the stratum in a three-dimensional seismic data body, and enabling the obtained slices to truly reflect the fault occurrence and trend in the period, the purpose of fault interpretation is achieved, the defect that the conventional slices cannot truly reflect the spread of a fault system along the structural trend is overcome, the influence of human interpretation factors on the slice result along the stratum is avoided, faults are effectively identified, the seismic interpretation is more reasonable, and reliable basis is provided for the combination of a fracture system.
Particularly, grid division is carried out on a target area to obtain grid areas, single-point stability comparison is carried out on the grid areas, whether the grid areas meet stratum requirements of well position deployment or not is determined, grid areas which do not meet the stratum requirements of well position deployment are omitted, the range of the target area is reduced, grid areas are renumbered according to determination results of the grid areas, the coincidence score and the non-coincidence score are calculated according to the number results, so that an overall stability evaluation value of the target area is obtained, whether the target area needs to be replaced or not is determined according to the overall stability evaluation value in combination with an overall stability standard evaluation value, stability of the target area is determined from the single point aspect and the overall aspect, therefore, proper positions of well position deployment are determined, well position deployment inappropriateness caused by inaccurate data obtained by considering the stability of the target area on one hand is avoided, and well position deployment accuracy is improved.
Particularly, the obtained seismic zone areas are classified, the target areas with high seismic intensity levels are omitted, the target areas of well position deployment are further reduced, the accuracy of well position deployment is improved, the safety state level of the coincident positions or the abnormal state level of the areas is determined according to the number and the level of the seismic zones forming the coincident positions in the areas with low seismic intensity levels, the safety of the areas is further determined, whether the areas meet the requirements of well position deployment is judged, well position operation accidents are prevented, and the safety of well position operation is ensured.
Drawings
FIG. 1 is a flow chart of a method for determining irregular cross-sectional shape of a fault in a sand and mud formation according to an embodiment;
FIG. 2 is a flow chart of a method of fault interpretation of the seismic data for which reliability is required, according to an embodiment;
FIG. 3 is a schematic diagram of a mode of co-directional fault growth described in the examples;
FIG. 4 is a schematic diagram of a reverse fault growth mode in an embodiment;
fig. 5 is a schematic view of the opposite fault growth mode in the embodiment.
Detailed Description
In order that the objects and advantages of the invention will become more apparent, the invention will be further described with reference to the following examples; it should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
Preferred embodiments of the present invention are described below with reference to the accompanying drawings. It should be understood by those skilled in the art that these embodiments are merely for explaining the technical principles of the present invention, and are not intended to limit the scope of the present invention.
It should be noted that, in the description of the present invention, terms such as "upper," "lower," "left," "right," "inner," "outer," and the like indicate directions or positional relationships based on the directions or positional relationships shown in the drawings, which are merely for convenience of description, and do not indicate or imply that the apparatus or elements must have a specific orientation, be constructed and operated in a specific orientation, and thus should not be construed as limiting the present invention.
Furthermore, it should be noted that, in the description of the present invention, unless explicitly specified and limited otherwise, the terms "mounted," "connected," and "connected" are to be construed broadly, and may be either fixedly connected, detachably connected, or integrally connected, for example; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communication between two elements. The specific meaning of the above terms in the present invention can be understood by those skilled in the art according to the specific circumstances.
Referring to fig. 1-5, fig. 1 is a flowchart illustrating a method for determining an irregular cross-sectional shape of a normal fault of a sand mud layer according to an embodiment; FIG. 2 is a flow chart of a method of fault interpretation of the seismic data for which reliability is required, according to an embodiment; FIG. 3 is a schematic diagram of a mode of co-directional fault growth described in the examples; FIG. 4 is a schematic diagram of a reverse fault growth mode in an embodiment; fig. 5 is a schematic view of the opposite fault growth mode in the embodiment.
The invention provides a method for determining the irregular cross-section shape of a positive fault of a sand-mud stratum, which comprises the following steps of,
Step S1, obtaining a seismic section according to seismic data of a target area, and evaluating the seismic section so as to determine the reliability of the seismic data;
S2, determining the position of a seismic zone in a target area through seismic data with reliability reaching the requirement, and performing fault interpretation on the seismic data so as to determine the fault trend of the target area;
Step S3, determining the classification and the characteristics of a target area conveyor belt according to the fault trend, and establishing a three-dimensional visual model according to the classification and the characteristics of the conveyor belt so as to realize the fault space distribution form of the area;
And S4, grid division is carried out on the target area, the seismic section is optimized according to the position relation between each grid area and the seismic zone and the combination of fault space distribution form, and a basis is provided for well position deployment.
The invention uses geological stress, structure evolution and sediment filling analysis as guidance, adopts an irregular section interpretation technology based on mechanical analysis to interpret the section, and specifically uses analysis test information of drilled wells, coring and strata and open-air outcrop information to build a combined relation of stratum grilles and lithology of a research area, analyzes development characteristics and formation mechanisms of the research area, builds a fault real interpretation mode, provides an analysis and judgment technical method for the reliability of the seismic information of the research area, and provides an analysis and judgment method for the reliability of the seismic information of the research area.
Specifically, the evaluation criteria for evaluating the seismic section in this embodiment include the amplitude retention of the seismic section, the authenticity of the structural morphology of the seismic section, the accuracy of the fault homing of the seismic section,
When evaluating the amplitude preservation of the seismic section, the stratum structure and stratum lithology characteristics are determined according to the historical data obtained by collecting the well logging and coring data of the well in the research area and the lithology and field outcrop data, the amplitude preservation of the seismic section is determined by adopting the calibration of the synthetic record layer at the well point,
If the consistency coefficient of the synthetic record and the side-well seismic channel is more than or equal to 80%, judging that the amplitude preservation of the seismic section meets the requirement;
If the consistency coefficient of the synthetic record and the well side seismic channel is less than 80%, judging that the amplitude preservation of the seismic section does not meet the requirement;
when evaluating the authenticity of the structural form of the seismic section, calibrating a drilling geological stratification on the seismic section for well connection explanation, calculating a first difference absolute value S1, S1= |t1-t2| according to the reflection time t1 of an in-phase shaft where a target layer is and the corresponding geological stratification conversion time t2,
If S1 is less than or equal to S10, judging that the authenticity of the structural form of the seismic section meets the requirement;
if S1 is more than S10, judging that the authenticity of the structural form of the seismic section does not meet the requirement;
S10 is a set structural form authenticity evaluation value;
When evaluating the fault homing accuracy of the seismic section, the real drilling breakpoint position is compared with the seismic data breakpoint position,
If the real drilling breakpoint position is consistent with the seismic data breakpoint position, judging that the fault homing accuracy of the seismic section meets the requirement;
if the real drilling breakpoint position is inconsistent with the seismic data breakpoint position, judging that the fault homing accuracy of the seismic section does not meet the requirement.
Specifically, in this embodiment, if the amplitude-preserving property of the seismic section meets a requirement, it is determined that the amplitude-preserving property of the seismic section meets a single determination condition;
If the amplitude preservation of the seismic section does not meet the requirement, judging that the amplitude preservation of the seismic section does not meet a single determination condition;
If the authenticity of the structural form of the seismic section meets the requirement, judging that the structural form of the seismic section meets a single determination condition;
if the authenticity of the structural form of the seismic section does not meet the requirement, judging that the structural form of the seismic section does not meet a single determination condition;
if the accuracy of the fault homing of the seismic section meets the requirement, judging that the fault homing of the seismic section meets a single determination condition;
if the accuracy of the fault homing of the seismic section does not meet the requirement, judging that the fault homing of the seismic section does not meet a single determining condition;
If the amplitude preservation performance, the construction form and the fault homing of the seismic section all meet the single determination condition, judging that the reliability level of the seismic data is one level, and providing a basis for well position deployment of the target area;
if the coverage, the structural form and the fault homing of the seismic section meet the single determining condition, judging that the reliability level of the seismic data is two-level, and if part of error information exists in the seismic data, further evaluation is needed;
If the amplitude preservation performance, the structural morphology and the items meeting the single determination condition in the fault homing are smaller than two items, the reliability level of the seismic data is judged to be three-level, and a basis can not be provided for well position deployment of the target area.
In particular, the embodiment performs fault interpretation on the seismic data with the reliability level of one level, including,
Step S21, determining the arrangement of a broken area on the seismic section and analyzing the size of a vertical broken distance on the seismic section;
Step S22, connecting the break points of the upper and lower seismic reflection layer large break distance faults with the break points of the middle layer small break distance faults in the longitudinal direction to form an irregular section of the same fault;
And S23, realizing the relationship of the arrangement, the falling direction and the plane spread of the fracture area of the interlayer translation fault by identifying the translation fault, and performing the section closure interpretation by utilizing the longitudinal and transverse section closure.
The method and the device carry out slicing along the structural trend surface, develop slicing research along the rule of stratum development, develop horizon grid interpretation on a target layer in seismic data, determine the time range of the target layer, obtain the structural trend of the stratum, carry out slicing along the structural trend surface of the stratum in a three-dimensional seismic data body, and the obtained slicing can truly reflect the fault occurrence and trend in the period, thereby achieving the purpose of fault interpretation, overcoming the defect that the traditional slicing cannot truly reflect the spread of a fault system along the structural trend, avoiding the influence of human interpretation factors on the boundary slice result, effectively identifying faults, enabling the seismic interpretation to be more reasonable and providing reliable basis for the combination of fracture systems.
Specifically, the classification of the conveyor belt in this embodiment includes a conveyor belt of the same direction type, a conveyor belt of the opposite direction type,
The homodromous transmission belt is formed by two or more faults with the same trend and approaching each other on a plane; which can be divided into two cases, including a first and a second co-directional fault growth mode,
The first co-directional fault growth mode is that two faults are basically positioned on a straight line, and are naturally connected together or overlapped in the transverse growth lengthening process, and the trend of the tail ends of the faults is not obviously changed;
The second equidirectional fault growth mode is that two faults are not located on the same straight line, the faults are close to each other and overlap in the growth process along the trend, a secondary transverse transmission fault or a main fault trend is generated in the close and overlapping section and is changed, and the faults are connected with each other to form a transmission belt;
the back-facing conveyor belt is characterized in that two adjacent or overlapping faults tend to be opposite and the trend is controlled by basically the same fracture; which can be divided into two cases, including a first reverse-facing fault growth mode and a second reverse-facing fault growth mode,
The first reverse fault growth mode is that two faults are connected with each other to form a similar reverse structure;
The second backward fault growth mode is that two faults are close to each other, a public lower disc is lifted, an upper disc is settled to form a half cutting, a barrier structure is integrally formed, or strain is redistributed among the faults due to fault interaction, and a small fault which is inclined or vertical to a main fault is formed on the public lower disc, so that a small cutting structure is formed;
The opposite type transmission belt is composed of two or two groups of faults which have basically the same trend and are inclined to face each other, and the opposite type transmission belt can be divided into two cases, including a first opposite type fault growth mode and a second opposite type fault growth mode;
the first opposite fault growth mode is that two faults are overlapped, and a transmission belt is formed at a public descent disc, wherein if a cutting structure is formed by common control of two groups of faults, the two faults can be regarded as a combination of the same-direction transmission belt and the opposite-direction transmission belt;
the second opposite fault growth mode is that the tail ends of two faults are close to or partially overlapped and are not connected, anticline between the two faults is obliquely spread, and a deposited water system mainly extends along the root of the fault on the side and extends from the tail ends of the faults to the low-topography part of the fault descent disk.
Specifically, in this embodiment, the target area is meshed, including, a first mesh area A1, a second mesh area A2 … …, and An nth mesh area An, for which An earthquake magnitude evaluation value D0 is set for any mesh area, a historical magnitude value D of each mesh area is obtained, and whether the requirement of stratum stability of well position deployment is met is determined according to the obtained relationship between the historical magnitude value D of each mesh area and the preset earthquake magnitude evaluation value D0, where in the determining process, single-point stability comparison and overall stability comparison are required.
Specifically, the seismic magnitude evaluation value D0 is set for the i-th mesh region Ai, the historical magnitude values Di, i=1, 2,..n,
The single-point stability comparison is carried out on the ith grid region Ai, a second difference absolute value S2i, S2 i= |Di-D0| between the historical magnitude Di of the ith grid region Ai and the seismic magnitude evaluation value D0 is calculated,
If S2i is less than or equal to S2i0, judging that the stratum stability at the grid area Ai meets the stratum requirement of well position deployment;
If S2i is more than S2i0 and Di is more than D0, judging that the stratum stability at the grid area Ai does not meet the stratum requirement of well position deployment;
If S2i is more than S2i0 and Di is less than D0, judging that the stratum stability at the grid area Ai meets the stratum requirement of well position deployment;
where S2i0 is a single-point stability evaluation value at the i-th mesh region Ai.
Specifically, the stability states of all grid areas are integrated, the overall stability of the stratum of the target area is determined, the number of grid areas meeting the stratum requirement of well site deployment is obtained, the meeting score is calculated, the number of grid areas not meeting the stratum requirement of well site deployment is obtained, the non-meeting score is calculated, including,
1) Renumbering the number of conforming grid areas, denoted as first conforming area a11, second conforming area a12 … … mth conforming area A1m,
2) The number of non-conforming grid areas is renumbered, denoted as first non-conforming area a21, second non-conforming area a22 … … p-th non-conforming area A2p,
Judging whether the overall stability of the target area meets the stratum requirement of well position deployment, calculating an overall stability evaluation value K, setting,
K=y-/>
Wherein,Score for the y-th coincidence region,/>Y is a first calculation compensation parameter of the coincidence score of the y-th coincidence region to the overall stability evaluation value K; s22x is the non-conforming score of the x-th non-conforming area, and Dx is the second calculated compensation parameter of the non-conforming score of the x-th non-conforming area to the overall stability evaluation value K;
If K is more than or equal to K0, judging that the overall stability of the target area meets the stratum requirement of well position deployment;
If K is less than K0, judging that the overall stability of the target area does not meet the stratum requirement of well position deployment, and replacing the target area to perform well position deployment;
Wherein K0 is an overall stability standard evaluation value.
According to the embodiment, grid division is carried out on the target area to obtain each grid area, single-point stability comparison is carried out on each grid area, whether each grid area meets stratum requirements of well position deployment or not is determined, grid areas which do not meet the stratum requirements of well position deployment are omitted, the range of the target area is reduced, grid areas are renumbered according to determination results of the grid areas, the coincidence score and the non-coincidence score are calculated according to the number results, so that an overall stability evaluation value of the target area is obtained, whether the target area needs to be replaced or not is determined according to the overall stability evaluation value in combination with an overall stability standard evaluation value, stability of the target area is determined from the single point aspect and the overall aspect, therefore, proper positions of well position deployment are determined, well position deployment inappropriateness caused by inaccurate data obtained by considering the stability of the target area on one hand is avoided, and well position deployment accuracy is improved.
Specifically, in this embodiment, each grid area that does not meet the stratum requirement of well location deployment is spliced to form a plurality of seismic zones,
For any seismic zone, calculating the average intensity of the seismic zone according to the intensity of each seismic source in the seismic zone, determining the vibration range of the seismic zone according to the average intensity, taking the centroid of the seismic zone as a round point, acquiring the distance between the farthest edge point of the seismic zone and the centroid, equally dividing the acquired distance in three directions to obtain dividing lengths, taking the centroid as the round point, sequentially rounding the dividing lengths by 2 times, dividing lengths by 3 times as radius, determining the seismic zones of different grades, including a first-stage seismic zone, a second-stage seismic zone and a third-stage seismic zone,
The first-level earthquake area is an area with the largest earthquake intensity, and the area does not meet stratum requirements of well position deployment;
the secondary earthquake area is an area with larger earthquake intensity, and when judging whether the area meets the stratum requirement of well position deployment, the safety state level of the area or the abnormal state level is required to be further judged;
The three-level earthquake area is an area with smaller earthquake intensity, and when judging whether the area meets the stratum requirement of well position deployment, the safety state level or the abnormal state level of the area needs to be further judged.
In particular, in this embodiment, for any implementation well,
If the well position deployment position is the secondary seismic area in the seismic zone, judging that the well position deployment position of the implementation well is in a first-level state of a safe state;
If the well position deployment position is the superposition position of the secondary seismic area and the tertiary seismic area in the two seismic zones, judging that the well position deployment position of the implementation well is in a secondary state of a safe state;
If the well position deployment position is the superposition position of each secondary seismic area in the two seismic zones, judging that the well position deployment position of the implementation well is in a first-stage state of an abnormal state;
If the well position deployment position is the superposition position formed by each seismic zone and the number of the seismic zones exceeds two, judging that the well position deployment position of the implementation well is in a secondary state of an abnormal state.
According to the embodiment, the obtained seismic zone areas are classified, the target areas with high seismic intensity levels are omitted, the target areas for well position deployment are further reduced, the well position deployment accuracy is improved, the safety state level of the coincident positions is determined according to the number and the level of the seismic zones forming the coincident positions in the areas with low seismic intensity levels, or the safety of the areas is further determined according to the abnormal state level, whether the areas meet the requirement for well position deployment is judged, well position operation accidents are prevented, and well position operation safety is ensured.
Thus far, the technical solution of the present invention has been described in connection with the preferred embodiments shown in the drawings, but it is easily understood by those skilled in the art that the scope of protection of the present invention is not limited to these specific embodiments. Equivalent modifications and substitutions for related technical features may be made by those skilled in the art without departing from the principles of the present invention, and such modifications and substitutions will be within the scope of the present invention.
The foregoing description is only of the preferred embodiments of the invention and is not intended to limit the invention; various modifications and variations of the present invention will be apparent to those skilled in the art. 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.
Claims (3)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311730804.2A CN117872472B (en) | 2023-12-15 | 2023-12-15 | A method for determining the irregular cross-sectional shape of normal faults in sandy mud strata |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311730804.2A CN117872472B (en) | 2023-12-15 | 2023-12-15 | A method for determining the irregular cross-sectional shape of normal faults in sandy mud strata |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117872472A CN117872472A (en) | 2024-04-12 |
| CN117872472B true CN117872472B (en) | 2024-06-18 |
Family
ID=90585480
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311730804.2A Active CN117872472B (en) | 2023-12-15 | 2023-12-15 | A method for determining the irregular cross-sectional shape of normal faults in sandy mud strata |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117872472B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118247280B (en) * | 2024-05-29 | 2024-07-30 | 中国海洋大学 | A method for simulating river cross-section morphology migration based on multi-source remote sensing data |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108411844A (en) * | 2018-03-12 | 2018-08-17 | 中南大学 | A kind of analysis of the debris flow velocity field of the irregular section of natural raceway groove and device |
| CN114910953A (en) * | 2022-04-15 | 2022-08-16 | 中国石油化工股份有限公司 | Interpretation method of irregular section of complex fracture zone |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7486589B2 (en) * | 2006-02-09 | 2009-02-03 | Schlumberger Technology Corporation | Methods and apparatus for predicting the hydrocarbon production of a well location |
| JP5770051B2 (en) * | 2011-09-05 | 2015-08-26 | 公益財団法人鉄道総合技術研究所 | Immediate estimation method of the rupture area of a large earthquake |
| WO2023141350A1 (en) * | 2022-01-24 | 2023-07-27 | Conocophillips Company | Systems and methods for nuclear magnetic resonance (nmr) well logging |
| CN116433002A (en) * | 2023-01-19 | 2023-07-14 | 华北水利水电大学 | A method and system for identification and stability evaluation of dangerous coal seams in multi-layer coal mined-out area |
| CN117192606A (en) * | 2023-09-25 | 2023-12-08 | 中南大学 | Method and system for optimizing support methods based on microseismic signal stability classification |
-
2023
- 2023-12-15 CN CN202311730804.2A patent/CN117872472B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108411844A (en) * | 2018-03-12 | 2018-08-17 | 中南大学 | A kind of analysis of the debris flow velocity field of the irregular section of natural raceway groove and device |
| CN114910953A (en) * | 2022-04-15 | 2022-08-16 | 中国石油化工股份有限公司 | Interpretation method of irregular section of complex fracture zone |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117872472A (en) | 2024-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | Architecture characterization of Ordovician fault-controlled paleokarst carbonate reservoirs in Tuoputai, Tahe oilfield, Tarim Basin, NW China | |
| CN111612899B (en) | Carbonate fracture-cavity type oil reservoir geological modeling method | |
| EP1779145B1 (en) | Method for creating facies probability cubes based upon geologic interpretation | |
| KR101642951B1 (en) | GIS-based real time earthquake prediction method | |
| Eberli et al. | Seismic Imaging of Carbonate Reservoirs and Systems: AAPG Memoir 81 | |
| El Ouahed et al. | Application of artificial intelligence to characterize naturally fractured zones in Hassi Messaoud Oil Field, Algeria | |
| CN114740547B (en) | A method for hierarchical zoning analysis and three-dimensional characterization of thrust-slip fault folds | |
| Kim et al. | A three-dimensional geotechnical spatial modeling method for borehole dataset using optimization of geostatistical approaches | |
| CN102918423A (en) | Method for seismic hydrocarbon system analysis | |
| CN117872472B (en) | A method for determining the irregular cross-sectional shape of normal faults in sandy mud strata | |
| CN113279746A (en) | Method for determining deformation risk area of casing and application | |
| CN106529762A (en) | New mine coal and gas outburst risk assessment method | |
| CN116449423A (en) | Method, device, equipment and storage medium for analyzing single sand body configuration of braided river | |
| CN116680957A (en) | A method for predicting the range of surface movement induced by mining in underground metal mines | |
| Bueno et al. | Structural and stratigraphic feature delineation and facies distribution using seismic attributes and well log analysis applied to a Brazilian carbonate field | |
| CN109469477B (en) | Method and device for predicting extension direction of artificial crack | |
| CN119313835A (en) | A method and system for automatic geological model building and dynamic accuracy correction | |
| CN114690245B (en) | A forward modeling method for geological model of seismic anomalies without drilling | |
| CN118242076A (en) | Fine development method for fracture-cavity carbonate rock oil-gas reservoir | |
| Svanes et al. | Integration of subsurface applications to develop a dynamic stochastic modeling workflow | |
| CN111815769B (en) | Modeling method, computing device and storage medium for thrust covered zone construction | |
| Maisons et al. | Induced microseismicity and procedure for closure of brine production caverns | |
| Singeisen | Initiation and mechanisms of rock slope failures during the 2016 Mw 7.8 Kaikōura earthquake, New Zealand | |
| CN118657387B (en) | A real-time water hazard assessment method and system for mine working face and mine-related areas | |
| CN114859406B (en) | Method for predicting structural hydrocarbon reservoir, electronic equipment and medium |
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 |