CN113495294B - Quantitative characterization and evaluation method and device for sliding fracture - Google Patents

Quantitative characterization and evaluation method and device for sliding fracture Download PDF

Info

Publication number
CN113495294B
CN113495294B CN202010265385.XA CN202010265385A CN113495294B CN 113495294 B CN113495294 B CN 113495294B CN 202010265385 A CN202010265385 A CN 202010265385A CN 113495294 B CN113495294 B CN 113495294B
Authority
CN
China
Prior art keywords
fracture
sliding
attribute
dimensional
horizon
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
Application number
CN202010265385.XA
Other languages
Chinese (zh)
Other versions
CN113495294A (en
Inventor
李世银
杨凤英
沈春光
关宝珠
于红枫
赵龙飞
王彭
刘瑞东
何军
熊昶
李会元
丁肇媛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202010265385.XA priority Critical patent/CN113495294B/en
Publication of CN113495294A publication Critical patent/CN113495294A/en
Application granted granted Critical
Publication of CN113495294B publication Critical patent/CN113495294B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/301Analysis for determining seismic cross-sections or geostructures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/301Analysis for determining seismic cross-sections or geostructures
    • G01V1/302Analysis for determining seismic cross-sections or geostructures in 3D data cubes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6169Data from specific type of measurement using well-logging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/64Geostructures, e.g. in 3D data cubes
    • G01V2210/641Continuity of geobodies
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/64Geostructures, e.g. in 3D data cubes
    • G01V2210/643Horizon tracking

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention provides a quantitative characterization and evaluation method and device for sliding fracture, which acquire state data according to drilled well and three-dimensional post-stack fidelity seismic data, wherein the state data comprises the following steps: constructing three-dimensional seismic data after guide filtering, horizon interpretation results, three-dimensional body attributes, layer root mean square attributes, fracture interpretation results, fracture plane combinations and residual trend surface data; converting the state data to obtain quantitative characterization characteristic parameters of the slip fracture; quantitatively characterizing the sliding fracture according to the sliding fracture quantitative characterization characteristic parameters to obtain sliding fracture evaluation characteristic parameters, wherein the sliding fracture evaluation characteristic parameters comprise: fracture type, fracture level, fracture strain type segmentation; according to the characteristic parameters of the sliding fracture evaluation, the result of the sliding fracture evaluation is obtained, and the accuracy of the sliding fracture evaluation is improved.

Description

Quantitative characterization and evaluation method and device for sliding fracture
Technical Field
The invention relates to the computer technology, in particular to a quantitative characterization and evaluation method and device for slip fracture.
Background
The slip fracture refers to the fracture generated by the relative horizontal movement of two fault discs under the action of couple of force under the action of torsional stress or shear stress field of the crust; it is also called a transverse fault, and can cause horizontal sliding of two sides of the fault relative to each other. In the process of oil gas formation and accumulation, a sliding fracture geological structure generally affects the transportation of oil gas and the formation of a reservoir, and in particular, in a carbonate stratum, karst action often develops along a sliding fracture zone to form a karst body reservoir body, so that an oil gas accumulation space is formed. In the aspect of oil and gas reservoir formation, deep and large sliding fracture is a main medium for communicating deep-barren hydrocarbon source rocks, is a dominant channel for oil and gas migration, the fracture activity is favorable for the hydrocarbon source rocks and an oil and gas accumulation area to form a space configuration relation, and the fracture multi-stage activity has an important influence on the oil and gas migration path, the reservoir formation period and the space distribution of the oil and gas reservoir, so that a high-quality reservoir and an oil and gas enrichment area need to be determined through sliding fracture in oil and gas reservoir exploration.
In the prior art, the skid fracture is qualitatively identified only through various methods or only a certain characteristic parameter of the skid fracture is quantitatively characterized, and then the skid fracture is evaluated through the obtained result.
However, in the prior art, the accuracy of evaluating the slip break is low.
Disclosure of Invention
The embodiment of the invention provides a quantitative characterization and evaluation method and device for a sliding fracture, which improve the accuracy of evaluating the sliding fracture.
In a first aspect of the embodiment of the present invention, a method for quantitatively characterizing and evaluating a skid fracture is provided, including:
acquiring state data according to the drilled well and the three-dimensional post-stack fidelity seismic data, wherein the state data comprises: constructing three-dimensional seismic data after guide filtering, horizon interpretation results, three-dimensional body attributes, layer root mean square attributes, fracture interpretation results, fracture plane combinations and residual trend surface data;
performing conversion processing on the state data to obtain quantitative characterization characteristic parameters of the slip fracture, wherein the quantitative characterization characteristic parameters of the slip fracture comprise: longitudinal break horizon, plane extension length, average fracture zone width, average break distance, fracture density, and remaining trend surface type;
quantitatively characterizing the sliding fracture according to the sliding fracture quantitative characterization characteristic parameters to obtain sliding fracture evaluation characteristic parameters, wherein the sliding fracture evaluation characteristic parameters comprise: fracture type, fracture level, fracture strain type segmentation;
And acquiring an evaluation result of the sliding fracture according to the sliding fracture evaluation characteristic parameters.
In a second aspect of the embodiments of the present invention, there is provided a device for quantitatively characterizing and evaluating a skid fracture, including:
the state data module is used for acquiring state data according to the drilled well and the three-dimensional post-stack fidelity seismic data, wherein the state data comprises: constructing three-dimensional seismic data after guide filtering, horizon interpretation results, three-dimensional body attributes, layer root mean square attributes, fracture interpretation results, fracture plane combinations and residual trend surface data;
the characteristic parameter module is used for carrying out conversion processing on the state data to obtain a quantitative characterization characteristic parameter of the sliding fracture, wherein the quantitative characterization characteristic parameter of the sliding fracture comprises: longitudinal break horizon, plane extension length, average fracture zone width, average break distance, fracture density, and remaining trend surface type;
the quantitative characterization module is used for quantitatively characterizing the sliding fracture according to the sliding fracture quantitative characterization characteristic parameter to obtain a sliding fracture evaluation characteristic parameter, wherein the sliding fracture evaluation characteristic parameter comprises: fracture type, fracture level, fracture strain type segmentation;
And the evaluation module is used for acquiring the evaluation result of the sliding fracture according to the characteristic parameter of the sliding fracture evaluation.
In a third aspect of the embodiments of the present invention, there is provided a skid break quantitative characterization and evaluation apparatus, including: a memory, a processor and a computer program stored in the memory, the processor running the computer program to perform the method of the first aspect and the various possible designs of the first aspect.
In a fourth aspect of embodiments of the present invention, there is provided a readable storage medium having stored therein a computer program for implementing the method of the first aspect and the various possible designs of the first aspect when the computer program is executed by a processor.
The invention provides a quantitative characterization and evaluation method and a quantitative characterization and evaluation device for sliding fracture, which acquire state data according to drilled well and three-dimensional post-stack fidelity seismic data, wherein the state data comprises the following steps: constructing three-dimensional seismic data after guide filtering, horizon interpretation results, three-dimensional body attributes, layer root mean square attributes, fracture interpretation results, fracture plane combinations and residual trend surface data; performing conversion processing on the state data to obtain quantitative characterization characteristic parameters of the slip fracture, wherein the quantitative characterization characteristic parameters of the slip fracture comprise: longitudinal break horizon, plane extension length, average fracture zone width, average break distance, fracture density, and remaining trend surface type; quantitatively characterizing the sliding fracture according to the sliding fracture quantitative characterization characteristic parameters to obtain sliding fracture evaluation characteristic parameters, wherein the sliding fracture evaluation characteristic parameters comprise: fracture type, fracture level, fracture strain type segmentation; and acquiring an evaluation result of the sliding fracture according to the sliding fracture evaluation characteristic parameters. The method comprises the steps of converting state data in oil and gas reservoir exploration and development, and obtaining parameters for quantitatively representing the slip fracture characteristics; the quantitative characterization of the sliding fracture is realized by utilizing the quantitative characterization characteristic parameters of the sliding fracture, the characteristic parameters of the sliding evaluation are further obtained on the basis, the sliding fracture is finally evaluated by utilizing the characteristic parameters of the sliding fracture, and the accuracy of evaluating the sliding fracture is improved. The method effectively guides the implementation of trap and the optimization of high-efficiency well positions, improves the drilling success rate and the contribution rate of single well productivity, guides the optimization of well type and drilling track, improves the drilling meeting rate of a reservoir and avoids the drilling risk. The method has good application effect and promotion value to the same industry in the oil and gas reservoir exploration and development of fracture-cavity type carbonate reservoirs.
Drawings
Fig. 1 is a schematic view of an application scenario provided in an embodiment of the present invention;
FIG. 2 is a schematic flow chart of a quantitative characterization and evaluation method for slip fracture provided by an embodiment of the invention;
FIG. 3 is a schematic structural diagram of a quantitative characterization and evaluation device for skid fracture provided by an embodiment of the invention;
fig. 4 is a schematic hardware structure diagram of a skid fracture quantization characterization and evaluation device according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are only some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The terms "first," "second," "third," "fourth" and the like in the description and in the claims and in the above drawings, if any, are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used may be interchanged where appropriate such that the embodiments of the invention described herein may be implemented in sequences other than those illustrated or otherwise described herein.
It should be understood that, in various embodiments of the present invention, the sequence number of each process does not mean that the execution sequence of each process should be determined by its functions and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present invention.
It should be understood that in the present invention, "comprising" and "having" and any variations thereof are intended to cover non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements that are expressly listed or inherent to such process, method, article, or apparatus.
It should be understood that in the present invention, "plurality" means two or more. "and/or" is merely an association relationship describing an association object, and means that three relationships may exist, for example, and/or B may mean: a exists alone, A and B exist together, and B exists alone. The character "/" generally indicates that the context-dependent object is an "or" relationship. "comprising A, B and C", "comprising A, B, C" means that all three of A, B, C comprise, "comprising A, B or C" means that one of the three comprises A, B, C, and "comprising A, B and/or C" means that any 1 or any 2 or 3 of the three comprises A, B, C.
It should be understood that in the present invention, "B corresponding to a", "a corresponding to B", or "B corresponding to a" means that B is associated with a, from which B can be determined. Determining B from a does not mean determining B from a alone, but may also determine B from a and/or other information. The matching of A and B is that the similarity of A and B is larger than or equal to a preset threshold value.
As used herein, "if" may be interpreted as "at … …" or "at … …" or "in response to a determination" or "in response to detection" depending on the context.
The technical scheme of the invention is described in detail below by specific examples. The following embodiments may be combined with each other, and some embodiments may not be repeated for the same or similar concepts or processes.
The slip fracture refers to the fracture generated by the relative horizontal movement of two fault discs under the action of couple of force under the action of torsional stress or shear stress field of the crust; it is also called a transverse fault, and can cause horizontal sliding of two sides of the fault relative to each other. In the process of oil gas formation and accumulation, a sliding fracture geological structure generally affects the transportation of oil gas and the formation of a reservoir, and in particular, in a carbonate stratum, karst action often develops along a sliding fracture zone to form a karst body reservoir body, so that an oil gas accumulation space is formed. In the aspect of oil and gas reservoir formation, deep and large sliding fracture is a main medium for communicating deep-barren hydrocarbon source rocks, is a dominant channel for oil and gas migration, the fracture activity is favorable for the hydrocarbon source rocks and an oil and gas accumulation area to form a space configuration relation, and the fracture multi-stage activity has an important influence on the oil and gas migration path, the reservoir formation period and the space distribution of the oil and gas reservoir, so that a high-quality reservoir and an oil and gas enrichment area need to be determined through sliding fracture in oil and gas reservoir exploration. In the prior art, the sliding fracture is qualitatively identified only through various methods or only a certain characteristic parameter of the sliding fracture is quantitatively characterized, then the sliding fracture is evaluated through the obtained result, and the sliding breaking system is not quantitatively characterized and evaluated. The current production practice proves that the characteristics of oil gas display, gas-oil ratio, water energy, single well yield, descending rule and the like of the drilling of different sliding fracture zones are obviously different, and the characteristics of the actual drilling of different sections of the same sliding fracture zone are also greatly different, so that the qualitative or single-parameter quantitative description of the sliding fracture in the prior art cannot meet the actual production requirement, and the trap evaluation, well position optimization and drilling track optimization cannot be accurately guided. Therefore, the accuracy of the prior art in evaluating the slip break is low.
Referring to fig. 1, an application scenario is schematically shown in the embodiment of the present invention. The measuring instrument 11 is used for measuring certain characteristic parameter data of the sliding fracture, then the server 12 is used for processing the acquired parameter data to acquire evaluation data of the sliding fracture, and finally the data is used for guiding trap evaluation, well position optimization and drilling track optimization. The current production practice proves that the characteristics of oil gas display, gas-oil ratio, water energy, single well yield, descending rule and the like of the drilling of different sliding fracture zones are obviously different, and the characteristics of the actual drilling of different sections of the same sliding fracture zone are also greatly different, so that the qualitative or single-parameter quantitative description of the sliding fracture in the prior art cannot meet the actual production requirement, and the trap evaluation, well position optimization and drilling track optimization cannot be accurately guided. Therefore, the accuracy of the prior art in evaluating the slip break is low.
Referring to fig. 2, a flowchart of a method for quantitatively characterizing and evaluating a skid break according to an embodiment of the present invention is shown, where an execution body of the method shown in fig. 2 may be a software and/or hardware device. The execution bodies of the present application may include, but are not limited to, at least one of: user equipment, network equipment, etc. The user equipment may include, but is not limited to, computers, smart phones, personal digital assistants (Personal Digital Assistant, abbreviated as PDA), and the above-mentioned electronic devices. The network device may include, but is not limited to, a single network server, a server group of multiple network servers, or a cloud of a large number of computers or network servers based on cloud computing, where cloud computing is one of distributed computing, and a super virtual computer consisting of a group of loosely coupled computers. The method comprises the steps of S101 to S104, and specifically comprises the following steps:
S101, acquiring state data according to drilled well and three-dimensional post-stack fidelity seismic data, wherein the state data comprises: constructing three-dimensional seismic data after guide filtering, horizon interpretation results, three-dimensional body attributes, layer root mean square attributes, fracture interpretation results, fracture plane combinations and residual trend surface data.
Specifically, three-dimensional seismic data, horizon interpretation results, three-dimensional body attributes, layer root mean square attributes, fracture interpretation results, fracture plane combinations and residual trend surface data after construction guide filtering are obtained according to the drilled and three-dimensional post-stack fidelity seismic data, so as to obtain quantitative characterization characteristic parameters of the sliding fracture. It can be understood that the scheme can perform conversion processing on the state data, and can obtain the quantitative characterization characteristic parameters of the slip fracture.
The following processes of acquiring three-dimensional seismic data, horizon interpretation results, three-dimensional body attributes, layer-following root mean square attributes, fracture interpretation results, fracture plane combinations and residual trend surface data after construction guide filtering are respectively described in detail, and specifically are as follows:
constructing guided filtered three-dimensional seismic data:
and acquiring the three-dimensional seismic data after the structure-oriented filtering according to the drilled well and the three-dimensional post-stack fidelity seismic data, wherein the three-dimensional seismic data after the structure-oriented filtering comprises the structure-oriented filtering processing of carrying out fracture enhancement on the three-dimensional post-stack fidelity seismic data, and acquiring the three-dimensional seismic data after the structure-oriented filtering.
Specifically, the three-dimensional post-stack fidelity seismic data, namely the three-dimensional seismic data, can be subjected to fracture-enhanced structural filtering treatment, and further improves the signal to noise ratio of the seismic data on the premise of not changing structural morphology, so that the continuity and the discontinuous characteristics of the same phase axis of the seismic data are more obvious, and a foundation is laid for later horizon interpretation, fracture interpretation and three-dimensional attribute body extraction.
Horizon interpretation results:
obtaining the horizon interpretation results according to the drilled and three-dimensional post-stack fidelity seismic data, including: performing well vibration calibration processing on the drilled well to obtain the horizon interpretation scheme; and then, according to the horizon interpretation scheme, conducting horizon fine interpretation on the three-dimensional seismic data body after construction guide filtering, and obtaining a horizon interpretation result.
Specifically, starting from the well drilling, establishing a corresponding relation between the same phase axis on the seismic section and the underground geological interface through fine well earthquake calibration, and finally determining a horizon interpretation scheme; and then according to the horizon interpretation scheme, conducting horizon fine interpretation on the three-dimensional seismic data body after construction guide filtering, and finishing fine interpretation of the top of a target horizon (oil and gas exploration development main force layer) and the top horizon of a chilblain system, wherein the interpretation precision is 1 interval of main survey lines and 1 interval of connecting lines, and the horizon interpretation result is obtained.
Three-dimensional volume attributes:
and acquiring the three-dimensional body attribute according to the drilled well and the three-dimensional post-stack fidelity seismic data.
Specifically, attribute extraction processing is performed on the three-dimensional seismic data after construction guide filtering, and the three-dimensional body attribute is obtained, wherein the three-dimensional body attribute comprises a third-generation coherent body attribute, a maximum positive curvature body attribute, a structure tensor body attribute and an ant detector attribute.
Layer-wise root mean square attribute:
and acquiring the layer-along root mean square attribute according to the drilled well and the three-dimensional post-stack fidelity seismic data.
Specifically, according to the three-dimensional body attribute and the layer interpretation result, the top edge layer root mean square attribute of the target layer is obtained, wherein the top edge layer root mean square attribute comprises a third-generation coherent edge layer root mean square attribute, a maximum positive curvature edge layer root mean square attribute, an ant detection edge layer root mean square attribute and a structure tensor edge layer root mean square attribute.
Fracture interpretation results:
and acquiring the fracture interpretation result according to the drilled well and the three-dimensional post-stack fidelity seismic data.
Specifically, according to the three-dimensional seismic data after the construction guide filtering processing and the isochronous slices of the third-generation coherence attribute, performing fine interpretation of the sliding fracture, wherein the interpretation precision is 2 lines apart from a main line and 2 lines apart from a connecting line, and obtaining the fracture interpretation result.
Fracture plane combination:
and acquiring the fracture plane combination according to the drilled well and the three-dimensional post-stack fidelity seismic data.
Specifically, according to the fracture interpretation result, the third-generation coherent layer-edge root mean square attribute of the top of the target layer and the maximum positive curvature layer-edge root mean square attribute of the top of the target layer, the target layer fracture plane combination is obtained.
Remaining trend surface data:
acquiring the residual trend surface data according to the drilled and three-dimensional post-stack fidelity seismic data, including: smoothing the interpretation result of the top level of the target layer to obtain a level after smoothing; and acquiring the residual trend surface data according to the difference value between the horizon interpretation result and the smoothed horizon.
Specifically, smoothing the top-level bit interpretation data of the destination layer, wherein smoothing parameters are 31 channels of main line intervals and 31 channels of tie line intervals, smoothing the top-level bit interpretation data of the destination layer continuously for 3 times according to the smoothing parameters to obtain level data after the top-level smoothing processing of the destination layer, and subtracting the level data after the smoothing processing from the top-level bit interpretation data of the destination layer to obtain the remaining trend surface data of the top-level of the destination layer.
S102, converting the state data to obtain quantitative characterization characteristic parameters of the sliding fracture.
Specifically, the slip fracture quantitative characterization characteristic parameters include longitudinal fracture horizon, plane extension length, average fracture zone width, average fracture distance, fracture density, and remaining trend surface type. It can be understood that after the state data is obtained, the state data can be converted to obtain the quantitative characterization characteristic parameters of the slip fracture.
The specific process of obtaining the longitudinal fracture horizon, the plane extension length, the average fracture zone width, the average fracture distance, the fracture density and the residual trend surface type is as follows:
longitudinal disconnect horizon:
and converting the state data to obtain the longitudinal disconnection horizon.
Specifically, according to the three-dimensional seismic data after construction guiding filtering, the third-generation coherence attribute and the horizon interpretation result, a longitudinal breaking horizon of the walk-slip fracture is determined, namely a place on the seismic section where the in-phase axis is obviously broken, twisted or the coherence value on the third-generation coherence attribute section is obviously reduced relative to the surrounding is a fracture development position, and then the longitudinal breaking horizon of the walk-slip fracture is determined according to the horizon interpretation result.
Plane extension length:
and converting the state data to obtain the plane extension length.
The method comprises the steps of determining the extension length of a top sliding fracture plane of a target layer according to a fracture plane combination result, and measuring a length value from a sliding fracture starting development position to a sliding fracture disappearance position along a main sliding fracture trend to obtain the extension length of the sliding fracture plane.
Average breaker width:
and converting the state data to obtain the average broken belt width.
Specifically, according to the root mean square attribute of the objective layer top structure tensor along the layer, referring to the three-dimensional seismic data after construction guide filtering, determining the sliding fracture breaking belt width of the objective layer top. The objective layer top structure tensor plane attribute value is obviously increased, and the position where the phase axis of the seismic section is disordered and broken is a fracture breaking area. And (3) perpendicular to the trend of the sliding fracture, starting from the development starting position of the sliding fracture, counting the widths of fracture zones at intervals of 2 tracks until the sliding fracture disappears, obtaining the width data of the fracture zones of the sliding fracture, and averaging the obtained width data of the fracture zones of the sliding fracture to obtain the average fracture zone width of the sliding fracture on the top of the target layer.
Average breaking distance:
and converting the state data to obtain the average break distance.
Specifically, according to the fracture interpretation result and the top layer position interpretation result of the destination layer, the method is perpendicular to the running and sliding fracture trend, and every 2 times of statistics of the top breaking distance of the destination layer are carried out from the running and sliding fracture starting development position until the running and sliding fracture disappears, running and sliding fracture breaking distance data are obtained, and then the running and sliding fracture breaking distance data are averaged to obtain the average breaking distance of the running and sliding fracture on the top of the destination layer.
Crack density:
and converting the state data to obtain the crack density.
Specifically, according to the root mean square attribute of the target layer top ant detection along the layer, counting the number of cracks in the range of 1 km (including 1 km) at the two sides of the target layer top sliding fracture, dividing the number of cracks by the plane extension length to obtain the number of cracks in the unit length range, and obtaining the crack density.
Remaining trend surface types:
and converting the state data to obtain the type of the residual trend surface.
Specifically, the type of the residual trend surface is obtained according to the data of the residual trend surface of the top of the target layer, wherein a positive value of the residual trend surface of the top of the target layer represents that the stratum bulge corresponds to the extrusion type, a negative value of the residual trend surface of the top of the target layer represents that the stratum bulge corresponds to the tension type, and a zero value of the residual trend surface of the top of the target layer represents that the stratum bulge does not exist or the stratum bulge corresponds to the concave translation type.
S103, quantitatively characterizing the sliding fracture according to the characteristic parameters of the sliding fracture, and obtaining the characteristic parameters of the sliding fracture evaluation.
Specifically, the slip fracture evaluation characteristic parameters include fracture type, fracture level, fracture strain type segmentation.
The specific process of obtaining the fracture type, fracture grade and fracture strain type segments is as follows:
fracture type:
and quantitatively characterizing the sliding fracture according to the characteristic parameters of the sliding fracture quantification characterization, and obtaining the fracture type.
Specifically, the fracture type is determined according to whether the sliding fracture longitudinal fracture horizon is fractured to a chilly line (hydrocarbon source rock development zone): the sliding fracture which is not broken to the chilla system (hydrocarbon source rock development area) in the longitudinal direction is non-oil source sliding fracture; the skid break in the longitudinal direction down to the chills (hydrocarbon-based rock development zone) is an oil source skid break.
Fracture grade:
and quantitatively characterizing the sliding fracture according to the sliding fracture quantitative characterization characteristic parameters, obtaining a sliding fracture quantitative characterization result, and obtaining the fracture level according to the sliding fracture quantitative characterization result and preset conditions.
Specifically, according to the quantitative characterization result of the sliding fracture and the preset condition, the sliding fracture is determined to be a primary sliding fracture, a secondary sliding fracture and a tertiary sliding fracture.
The classification process is as follows:
according to the quantitative characterization result of the skid break and the preset condition, determining the skid break as a primary skid break comprises the following steps:
the quantitative characterization result of the skid break meets the following two or more preset conditions to determine the skid break as a primary skid break,
the plane extends over a length of greater than 10 kilometers;
the average breaker belt width is greater than 200 meters;
the average break distance is greater than 50 meters;
the crack density is greater than 50.
According to the quantitative characterization result of the skid break and the preset condition, determining the skid break as a secondary skid break comprises the following steps:
the quantitative characterization result of the skid break meets the following two or more preset conditions to determine the skid break as a secondary skid break,
the plane extension length is greater than or equal to 5 kilometers and less than or equal to 10 kilometers;
the average breaker strip width is greater than or equal to 100 meters and less than or equal to 200 meters;
the average break distance is greater than or equal to 30 meters and less than or equal to 50 meters;
the crack density is greater than or equal to 30 and less than or equal to 50.
Dividing the sliding fracture into three-stage sliding fracture according to the sliding fracture quantitative characterization result and the preset condition, wherein the sliding fracture comprises the following steps:
And determining the skid break which does not meet the primary skid break preset condition and the secondary skid break preset condition as a tertiary skid break.
Fracture strain type segmentation:
and quantitatively characterizing the sliding fracture according to the sliding fracture quantitative characterization characteristic parameters, and obtaining the fracture strain type segment.
Specifically, the fracture type and the fracture level represent the differences of the fracture type and the level of the fracture of different sliding fracture, and the strain type in the same sliding fracture is also different, so that the same sliding fracture can be segmented according to the strain type difference, and particularly the fracture strain type of the same sliding fracture is segmented according to the type of the residual trend surface of the top of the target layer at the development position of the sliding fracture: the extrusion type residual trend surface corresponds to the extrusion section; the tension type residual trend surface corresponds to the tension section; the translation type residual trend surface corresponds to the translation segment.
S104, according to the sliding fracture evaluation characteristic parameters, acquiring an evaluation result of the sliding fracture.
Specifically, in practical application, different fracture types have different roles in source control, storage control and storage control, the oil source sliding fracture has the roles of source control, storage control and storage control, and the non-oil source sliding fracture only has the roles of storage control and storage control, so that the importance of the oil source sliding fracture in the fracture type is obviously higher than that of the non-oil source sliding fracture in the fracture type in oil and gas exploration and development practice.
In addition, the fracture grades are different in representation of the fracture activity intensity, the higher the fracture grade is, the stronger the fracture activity is, and the formation and the oil gas filling of a reservoir are facilitated, so that the primary sliding fracture is optimal, the secondary sliding fracture is inferior, and the tertiary sliding fracture is worst in the oil gas exploration and development practice.
The difference of different strain types in the same sliding fracture is also used for reservoir transformation and oil gas filling, the stress of the whole upper tensile section is released, and the tensile cracks around the fracture develop, so that the method is beneficial to reservoir transformation and oil gas filling; the stress of the extrusion section is concentrated, the cracks around the fracture are few and mainly compressive cracks, so that the reservoir transformation and the oil gas filling are not facilitated; the translation section is a stretch section and an extrusion section transition section, and the functions of reservoir transformation and oil gas filling are between the two sections, so that the stretch section is optimal, the translation section is inferior and the extrusion section is worst in the same sliding fracture in oil gas exploration and development practice.
Comprehensively referencing the slip fracture evaluation characteristic parameters: and (5) segmenting the fracture type, the fracture level and the fracture strain type, and obtaining the evaluation result of the sliding fracture. Firstly, evaluating the priority levels of different sliding fractures in oil and gas exploration and development practices: the primary oil source sliding fracture is greater than the secondary oil source sliding fracture, the tertiary oil source sliding fracture is greater than the primary non-oil source sliding fracture, the secondary non-oil source sliding fracture is greater than the tertiary non-oil source sliding fracture. Then, evaluating the priority levels of different fracture strain types of the same sliding fracture in the oil and gas exploration and development practice: the stretching section is larger than the translation section and larger than the extrusion section.
Illustratively, the oil and gas exploration and development can be guided in sections according to the fracture types, the fracture grades and the fracture strain types obtained above, and the specific steps are as follows:
the guide trap preferably: the higher the running and sliding fracture priority level of the trap is, the more beneficial the scale of a reservoir layer and the oil gas filling condition in the trap are, for example, the trap controlled by the running and sliding fracture of the primary oil source is more beneficial than the trap controlled by the running and sliding fracture of the secondary oil source; for the same trap controlled by the sliding fracture, the trap controlled by the stretching section is more beneficial than the trap controlled by the translation section, and the trap controlled by the translation section is more beneficial than the extrusion section.
The guiding well position is preferably: the preferred well position is selected according to the size of the carving volume of the single well and the distance from the sliding fracture in the preferred favorable trap. The single well carving volume is more than or equal to 80 square and has high efficiency potential from well positions with the vertical distance from the sliding fracture less than or equal to 1.5 km.
Pilot well type preference and well trajectory optimization: the drilling of the reservoir and the fracture-crack development area should be considered at the same time, and a horizontal well or an inclined well should be optimized to improve the direct drilling rate of the reservoir; before the well track is designed to enter the target layer, a fracture breaking area is avoided, and drilling risks are avoided as much as possible.
According to the quantitative characterization and evaluation method for the skid fracture, which is provided by the embodiment, the state data is obtained according to the drilled well and the three-dimensional post-stack fidelity seismic data, wherein the state data comprises the following steps: constructing three-dimensional seismic data after guide filtering, horizon interpretation results, three-dimensional body attributes, layer root mean square attributes, fracture interpretation results, fracture plane combinations and residual trend surface data; performing conversion processing on the state data to obtain quantitative characterization characteristic parameters of the slip fracture, wherein the quantitative characterization characteristic parameters of the slip fracture comprise: longitudinal break horizon, plane extension length, average fracture zone width, average break distance, fracture density, and remaining trend surface type; quantitatively characterizing the sliding fracture according to the sliding fracture quantitative characterization characteristic parameters to obtain sliding fracture evaluation characteristic parameters, wherein the sliding fracture evaluation characteristic parameters comprise: fracture type, fracture level, fracture strain type segmentation; and the accurate evaluation of the sliding fracture is improved according to the sliding fracture evaluation characteristic parameters. The method comprises the steps of converting state data in oil and gas reservoir exploration and development, and obtaining parameters for quantitatively representing the slip fracture characteristics; the quantitative characterization of the skid fracture is realized by utilizing the quantitative characterization characteristic parameters of the skid fracture, and the skid evaluation characteristic parameters are further obtained on the basis; the characteristic parameters of the sliding fracture evaluation are utilized to finally realize the sliding fracture evaluation, effectively guide the trap implementation and the efficient well position optimization, improve the well drilling success rate and the single well productivity contribution rate, guide the well type and the well drilling track optimization, and avoid the well drilling risk while improving the reservoir drilling meeting rate. The method has good application effect and promotion value to the same industry in the oil and gas reservoir exploration and development of fracture-cavity type carbonate reservoirs.
Referring to fig. 3, a schematic structural diagram of a device for quantitatively characterizing and evaluating a skid fracture according to an embodiment of the present invention is provided, where the device 30 for quantitatively characterizing and evaluating a skid fracture includes:
a status data module 31, configured to obtain status data according to the drilled well and the three-dimensional post-stack fidelity seismic data, where the status data includes: constructing three-dimensional seismic data after guide filtering, horizon interpretation results, three-dimensional body attributes, layer root mean square attributes, fracture interpretation results, fracture plane combinations and residual trend surface data;
the characteristic parameter module 32 is configured to perform a conversion process on the state data, and obtain a quantitative characterization characteristic parameter of the slip break, where the quantitative characterization characteristic parameter of the slip break includes: longitudinal break horizon, plane extension length, average fracture zone width, average break distance, fracture density, and remaining trend surface type;
the quantization characterization module 33 is configured to quantitatively characterize the skid fracture according to the skid fracture quantization characterization characteristic parameter, and obtain a skid fracture evaluation characteristic parameter, where the skid fracture evaluation characteristic parameter includes: fracture type, fracture level, fracture strain type segmentation;
And the evaluation module 34 is used for acquiring the evaluation result of the sliding fracture according to the sliding fracture evaluation characteristic parameter.
The quantitative characterization of the skid break in the embodiment shown in fig. 3, corresponding to the evaluation means, may be used to perform the steps in the method shown in fig. 2, the implementation principle and technical effect of which are similar, and will not be repeated here.
Optionally, the status data module 31 is specifically configured to:
and carrying out fracture enhancement construction guide filtering treatment on the three-dimensional post-stack fidelity seismic data to obtain the three-dimensional seismic data after construction guide filtering.
Optionally, the status data module 31 is specifically configured to:
and carrying out well vibration calibration processing on the drilled well to obtain a horizon interpretation scheme, and carrying out horizon fine interpretation on the three-dimensional seismic data body after construction guide filtering according to the horizon interpretation scheme to obtain a horizon interpretation result.
Optionally, the status data module 31 is specifically configured to:
and carrying out attribute extraction processing on the three-dimensional seismic data after construction guide filtering to obtain the three-dimensional body attribute, wherein the three-dimensional body attribute comprises a third-generation coherent body attribute, a maximum positive curvature body attribute, a structure tensor body attribute and an ant detector attribute.
Optionally, the status data module 31 is specifically configured to:
and acquiring the layer-following root mean square attribute according to the three-dimensional attribute and the horizon interpretation result, wherein the layer-following root mean square attribute comprises: third generation coherent layer-following root mean square attribute, maximum positive curvature layer-following root mean square attribute, ant detection layer-following root mean square attribute, structure tensor layer-following root mean square attribute.
Optionally, the status data module 31 is specifically configured to:
and carrying out fine interpretation processing of the sliding fracture according to the three-dimensional seismic data subjected to the construction guide filtering processing and the isochronous slices of the third-generation coherence attribute, and obtaining the fracture interpretation result.
Optionally, the status data module 31 is specifically configured to:
and acquiring the fracture plane combination according to the fracture interpretation result, the third-generation coherent layer-following root mean square attribute and the maximum positive curvature layer-following root mean square attribute.
Optionally, the status data module 31 is specifically configured to:
smoothing the horizon interpretation result to obtain a smoothed horizon;
and acquiring the residual trend surface data according to the difference value between the horizon interpretation result and the smoothed horizon.
Optionally, the quantitative characterization characteristic parameters of the sliding fracture include a longitudinal fracture horizon, a plane extension length, an average fracture zone width, an average fracture distance, a fracture density and a residual trend surface type;
the characteristic parameter module 32 is specifically configured to:
and converting the state data to obtain a longitudinal breaking horizon, a plane extension length, an average breaking belt width, an average breaking distance, a crack density and a residual trend surface type of the sliding fracture.
Optionally, the characteristic parameter module 32 is specifically configured to:
and acquiring the longitudinal break horizon according to the three-dimensional seismic data after construction guide filtering, the third-generation coherence attribute and the horizon interpretation result.
Optionally, the characteristic parameter module 32 is specifically configured to:
according to the fracture plane combination, the distance between the starting position and the vanishing position of the sliding fracture is obtained;
and determining the plane extension length according to the distance.
Optionally, the characteristic parameter module 32 is specifically configured to:
acquiring the number of cracks of the two sides of the sliding fracture in a preset range according to the ant detection layer-following root mean square attribute;
and obtaining the crack density according to the number of cracks and the plane extension length.
Optionally, the characteristic parameter module 32 is specifically configured to:
and acquiring the average broken belt width according to the structural tensor layer root mean square attribute and the three-dimensional seismic data after construction guide filtering.
Optionally, the characteristic parameter module 32 is specifically configured to:
counting the breaking distance of the sliding breaking according to the breaking interpretation result;
and acquiring the average breaking distance according to the breaking distance data.
Optionally, the characteristic parameter module 32 is specifically configured to:
and determining the type of the residual trend surface according to the residual trend surface data, wherein the residual trend surface data is positive and represents a squeezing type, the residual trend surface data is negative and represents a stretching type, and the residual trend surface data is zero and represents a translation type.
Optionally, the slip fracture evaluation characteristic parameters include fracture type, fracture level, and fracture strain type segmentation.
The quantization characterization module 33 is specifically configured to: and quantitatively characterizing the sliding fracture according to the sliding fracture quantitative characterization characteristic parameters, and acquiring the sliding fracture evaluation characteristic parameters. Wherein, the slip fracture evaluation characteristic parameters comprise: fracture type, fracture level, and fracture strain type segmentation.
Optionally, the quantization characterization module 33 is specifically configured to:
the fracture type is determined based on whether the walk fracture longitudinal fracture horizon breaks into the chills (hydrocarbon development zone).
Optionally, the quantization characterization module 33 is specifically configured to:
and determining the fracture level according to the sliding fracture quantitative characterization result and preset conditions.
The quantitative characterization result of the skid break meets the following two or more preset conditions to determine the skid break as a primary skid break,
the plane extends over a length of greater than 10 kilometers;
the average breaker belt width is greater than 200 meters;
the average break distance is greater than 50 meters;
the crack density is greater than 50.
Optionally, the quantization characterization module 33 is specifically configured to:
the quantitative characterization result of the skid break meets the following two or more preset conditions to determine the skid break as a secondary skid break,
the plane extension length is greater than or equal to 5 kilometers and less than or equal to 10 kilometers;
the average breaker strip width is greater than or equal to 100 meters and less than or equal to 200 meters;
the average break distance is greater than or equal to 30 meters and less than or equal to 50 meters;
the crack density is greater than or equal to 30 and less than or equal to 50.
Optionally, the quantization characterization module 33 is specifically configured to:
and determining the skid break which does not meet the primary skid break preset condition and the secondary skid break preset condition as a tertiary skid break.
Optionally, the quantization characterization module 33 is specifically configured to:
and determining the fracture strain type segmentation according to the residual trend surface type, wherein the extrusion type corresponds to the extrusion segment, the tension type corresponds to the tension segment, and the translation type corresponds to the translation segment.
Optionally, the evaluation module 34 is specifically configured to:
and acquiring an evaluation result of the sliding fracture according to the sliding fracture evaluation characteristic parameters.
Referring to fig. 4, a schematic hardware structure of a quantitative characterization and evaluation device for skid fracture according to an embodiment of the present invention is provided, where the quantitative characterization device 40 for skid fracture includes: a processor 41, a memory 42 and a computer program; wherein the method comprises the steps of
A memory 42 for storing the computer program, which may also be a flash memory (flash). Such as application programs, functional modules, etc. implementing the methods described above.
And a processor 41 for executing the computer program stored in the memory to implement the steps executed by the terminal in the above method. Reference may be made in particular to the description of the embodiments of the method described above.
Alternatively, the memory 42 may be separate or integrated with the processor 41.
When the memory 42 is a device separate from the processor 41, the apparatus may further include:
a bus 43 for connecting the memory 42 and the processor 41.
The present invention also provides a readable storage medium having stored therein a computer program for implementing the methods provided by the various embodiments described above when executed by a processor.
The readable storage medium may be a computer storage medium or a communication medium. Communication media includes any medium that facilitates transfer of a computer program from one place to another. Computer storage media can be any available media that can be accessed by a general purpose or special purpose computer. For example, a readable storage medium is coupled to the processor such that the processor can read information from, and write information to, the readable storage medium. In the alternative, the readable storage medium may be integral to the processor. The processor and the readable storage medium may reside in an application specific integrated circuit (Application Specific Integrated Circuits, ASIC for short). In addition, the ASIC may reside in a user device. The processor and the readable storage medium may reside as discrete components in a communication device. The readable storage medium may be read-only memory (ROM), random-access memory (RAM), CD-ROMs, magnetic tape, floppy disk, optical data storage device, etc.
The present invention also provides a program product comprising execution instructions stored in a readable storage medium. The at least one processor of the device may read the execution instructions from the readable storage medium, the execution instructions being executed by the at least one processor to cause the device to implement the methods provided by the various embodiments described above.
In the above embodiment of the apparatus, it should be understood that the processor may be a central processing unit (english: central Processing Unit, abbreviated as CPU), or may be other general purpose processors, digital signal processors (english: digital Signal Processor, abbreviated as DSP), application specific integrated circuits (english: application Specific Integrated Circuit, abbreviated as ASIC), or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of a method disclosed in connection with the present invention may be embodied directly in a hardware processor for execution, or in a combination of hardware and software modules in a processor for execution.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and not for limiting the same; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some or all of the technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit of the invention.

Claims (6)

1. The quantitative characterization and evaluation method for the skid fracture is characterized by comprising the following steps:
acquiring state data according to the drilled well and the three-dimensional post-stack fidelity seismic data, wherein the state data comprises: constructing three-dimensional seismic data after guide filtering, horizon interpretation results, three-dimensional body attributes, layer root mean square attributes, fracture interpretation results, fracture plane combinations and residual trend surface data;
performing conversion processing on the state data to obtain quantitative characterization characteristic parameters of the slip fracture, wherein the quantitative characterization characteristic parameters of the slip fracture comprise: longitudinal break horizon, plane extension length, average fracture zone width, average break distance, fracture density, and remaining trend surface type;
quantitatively characterizing the sliding fracture according to the sliding fracture quantitative characterization characteristic parameters to obtain sliding fracture evaluation characteristic parameters, wherein the sliding fracture evaluation characteristic parameters comprise: fracture type, fracture level, fracture strain type segmentation, the fracture type comprising: non-oil source walk-slip fracture and oil source walk-slip fracture, the fracture level includes a primary walk-slip fracture, a secondary walk-slip fracture and a tertiary walk-slip fracture, the fracture strain type segment includes: an extrusion section, a stretching section and a translation section;
Acquiring an evaluation result of the sliding fracture according to the sliding fracture evaluation characteristic parameters;
the obtaining the three-dimensional seismic data after the structure-oriented filtering according to the drilled well and the three-dimensional post-stack fidelity seismic data comprises the following steps:
performing fracture-enhanced structure-oriented filtering processing on the three-dimensional post-stack fidelity seismic data to obtain three-dimensional seismic data after structure-oriented filtering;
the obtaining the horizon interpretation result according to the drilled well and the three-dimensional post-stack fidelity seismic data comprises the following steps:
performing well vibration calibration processing on the drilled well to obtain a horizon interpretation scheme;
according to the horizon interpretation scheme, conducting horizon fine interpretation on the three-dimensional seismic data after construction guide filtering, and obtaining a horizon interpretation result;
the obtaining the three-dimensional body attribute according to the drilled well and the three-dimensional post-stack fidelity seismic data comprises:
performing attribute extraction processing on the three-dimensional seismic data after construction guide filtering to obtain the three-dimensional body attribute, wherein the three-dimensional body attribute comprises a third-generation coherent body attribute, a maximum positive curvature body attribute, a structure tensor body attribute and an ant detector attribute;
the obtaining the fracture interpretation result according to the drilled well and the three-dimensional post-stack fidelity seismic data comprises the following steps:
According to the three-dimensional seismic data of the structure-oriented filtering and the isochronous slice of the third-generation coherence attribute, performing fine interpretation of the sliding fracture to obtain a fracture interpretation result;
the obtaining the layer-along root mean square attribute according to the drilled well and the three-dimensional post-stack fidelity seismic data comprises the following steps:
acquiring the layer-following root mean square attribute according to the three-dimensional body attribute and the horizon interpretation result, wherein the layer-following root mean square attribute comprises a third-generation coherent layer-following root mean square attribute, a maximum positive curvature layer-following root mean square attribute, an ant detection layer-following root mean square attribute and a structure tensor layer-following root mean square attribute;
the obtaining the fracture plane combination according to the drilled well and the three-dimensional post-stack fidelity seismic data comprises the following steps:
acquiring the fracture plane combination according to the fracture interpretation result, the third-generation coherent layer-following root mean square attribute and the maximum positive curvature layer-following root mean square attribute;
the obtaining the remaining trend surface data according to the drilled well and the three-dimensional post-stack fidelity seismic data comprises the following steps:
smoothing the horizon interpretation result to obtain a smoothed interpretation result;
acquiring the residual trend surface data according to the difference between the horizon interpretation result and the smooth interpretation result;
The converting the state data to obtain the longitudinal disconnection horizon includes:
acquiring the longitudinal break horizon according to the three-dimensional seismic data after construction guide filtering, the third-generation coherence body attribute and the horizon interpretation result;
the converting the state data to obtain the plane extension length includes:
according to the fracture plane combination, obtaining the distance between the starting position and the disappearing position of the sliding fracture;
determining the plane extension length according to the distance;
the converting the state data to obtain the crack density includes:
acquiring the number of cracks of the two sides of the sliding fracture in a preset range according to the layer-following root mean square attribute;
acquiring the crack density according to the number of cracks and the plane extension length;
the converting the state data to obtain the average broken belt width comprises the following steps:
acquiring the average broken belt width according to the structure tensor layer root mean square attribute and the three-dimensional seismic data after construction guide filtering;
the converting the state data to obtain the average break distance includes:
According to the fracture interpretation result, counting the breaking distance of the sliding fracture, and obtaining breaking distance data of the sliding fracture;
acquiring the average break distance according to the break distance data;
the converting the state data to obtain the remaining trend surface type includes:
determining the type of the residual trend surface according to the residual trend surface data, wherein the positive value of the residual trend surface data represents a squeezing type, the negative value of the residual trend surface data represents a stretching type, and the zero value of the residual trend surface data represents a translation type;
the step of quantitatively characterizing the step of sliding fracture according to the characteristic parameters of quantitatively characterizing the step of sliding fracture, and the step of obtaining the fracture type comprises the following steps:
acquiring the position relationship between the longitudinal breaking horizon and the hydrocarbon source rock development area;
acquiring the fracture type according to the position relation;
the step of quantitatively characterizing the step of sliding fracture according to the characteristic parameters of quantitatively characterizing the step of sliding fracture, and the step of obtaining the fracture level comprises the following steps:
quantitatively representing the sliding fracture according to the plane extension length, the average broken belt width, the average breaking distance and the crack density, and obtaining a sliding fracture quantitative representation result;
Acquiring the fracture level according to the sliding fracture quantitative characterization result and preset conditions;
the step of quantitatively characterizing the step fracture according to the step fracture quantitatively characterizing characteristic parameters to obtain the fracture strain type segment comprises the following steps:
and acquiring the sliding fracture strain type segmentation according to the residual trend surface type, wherein the extrusion type corresponds to the extrusion section, the tension type corresponds to the tension section, and the translation type corresponds to the translation section.
2. The method of claim 1, wherein the obtaining the primary skid break according to the skid break quantitative characterization result and the preset condition comprises:
the quantitative characterization result of the skid break meets the following two or more preset conditions to determine the skid break as the primary skid break,
the plane extends over a length of greater than 10 kilometers;
the average breaker belt width is greater than 200 meters;
the average break distance is greater than 50 meters;
the crack density is greater than 50.
3. The method according to claim 2, wherein the obtaining the secondary skid break according to the skid break quantitative characterization result and the preset condition comprises:
The quantitative characterization result of the skid break meets the following two or more preset conditions to determine the skid break as the secondary skid break,
the plane extension length is greater than or equal to 5 kilometers and less than or equal to 10 kilometers;
the average breaker strip width is greater than or equal to 100 meters and less than or equal to 200 meters;
the average break distance is greater than or equal to 30 meters and less than or equal to 50 meters;
the crack density is greater than or equal to 30 and less than or equal to 50.
4. The method of claim 3, wherein the obtaining the three-level skid break according to the skid break quantitative characterization result and the preset condition comprises:
and determining the sliding fracture which does not meet the preset conditions of the primary sliding fracture and the secondary sliding fracture as the tertiary sliding fracture.
5. The utility model provides a walk slip fracture quantization characterization and evaluation device which characterized in that includes:
the state data module is used for acquiring state data according to the drilled well and the three-dimensional post-stack fidelity seismic data, wherein the state data comprises: constructing three-dimensional seismic data after guide filtering, horizon interpretation results, three-dimensional body attributes, layer root mean square attributes, fracture interpretation results, fracture plane combinations and residual trend surface data;
The characteristic parameter module is used for carrying out conversion processing on the state data to obtain a quantitative characterization characteristic parameter of the sliding fracture, wherein the quantitative characterization characteristic parameter of the sliding fracture comprises: longitudinal break horizon, plane extension length, average fracture zone width, average break distance, fracture density, and remaining trend surface type;
the quantitative characterization module is used for quantitatively characterizing the sliding fracture according to the sliding fracture quantitative characterization characteristic parameter to obtain a sliding fracture evaluation characteristic parameter, wherein the sliding fracture evaluation characteristic parameter comprises: fracture type, fracture level, fracture strain type segmentation, the fracture type comprising: non-oil source walk-slip fracture and oil source walk-slip fracture, the fracture level includes a primary walk-slip fracture, a secondary walk-slip fracture and a tertiary walk-slip fracture, the fracture strain type segment includes: an extrusion section, a stretching section and a translation section;
the evaluation module is used for acquiring an evaluation result of the sliding fracture according to the sliding fracture evaluation characteristic parameters;
the state data module is specifically configured to:
performing fracture-enhanced structure-oriented filtering processing on the three-dimensional post-stack fidelity seismic data to obtain three-dimensional seismic data after structure-oriented filtering;
The state data module is specifically configured to:
performing well vibration calibration processing on the drilled well to obtain a horizon interpretation scheme, and performing horizon fine interpretation on the three-dimensional seismic data after construction guide filtering according to the horizon interpretation scheme to obtain a horizon interpretation result;
the state data module is specifically configured to:
performing attribute extraction processing on the three-dimensional seismic data after construction guide filtering to obtain the three-dimensional body attribute, wherein the three-dimensional body attribute comprises a third-generation coherent body attribute, a maximum positive curvature body attribute, a structure tensor body attribute and an ant detector attribute;
the state data module is specifically configured to:
and acquiring the layer-following root mean square attribute according to the three-dimensional attribute and the horizon interpretation result, wherein the layer-following root mean square attribute comprises: third generation coherent layer-following root mean square attribute, maximum positive curvature layer-following root mean square attribute, ant detection layer-following root mean square attribute, and structure tensor layer-following root mean square attribute;
the state data module is specifically configured to:
according to the three-dimensional seismic data of the structure-oriented filtering and the isochronous slice of the third-generation coherence attribute, carrying out the fine interpretation processing of the sliding fracture to obtain the fracture interpretation result;
The state data module is specifically configured to:
acquiring the fracture plane combination according to the fracture interpretation result, the third-generation coherent layer-following root mean square attribute and the maximum positive curvature layer-following root mean square attribute;
the state data module is specifically configured to:
smoothing the horizon interpretation result to obtain a smoothed horizon;
acquiring the residual trend surface data according to the difference between the horizon interpretation result and the smoothed horizon;
the characteristic parameter module is specifically used for:
acquiring the longitudinal break horizon according to the three-dimensional seismic data after construction guide filtering, the third-generation coherence body attribute and the horizon interpretation result;
the characteristic parameter module is specifically used for:
according to the fracture plane combination, the distance between the starting position and the vanishing position of the sliding fracture is obtained;
determining the plane extension length according to the distance;
the characteristic parameter module is specifically used for:
acquiring the number of cracks of the two sides of the sliding fracture in a preset range according to the ant detection layer-following root mean square attribute;
acquiring the crack density according to the number of cracks and the plane extension length;
The characteristic parameter module is specifically used for:
acquiring the average broken belt width according to the structure tensor layer root mean square attribute and the three-dimensional seismic data after construction guide filtering;
the characteristic parameter module is specifically used for:
counting the breaking distance of the sliding breaking according to the breaking interpretation result;
acquiring the average break distance according to the break distance data;
the characteristic parameter module is specifically used for:
determining the type of the residual trend surface according to the residual trend surface data, wherein the positive value of the residual trend surface data represents a squeezing type, the negative value of the residual trend surface data represents a stretching type, and the zero value of the residual trend surface data represents a translation type;
the quantization characterization module is specifically configured to:
determining the fracture type according to whether the sliding fracture longitudinal fracture horizon is fractured to a hydrocarbon-based rock development zone;
the quantization characterization module is specifically configured to:
quantitatively representing the sliding fracture according to the plane extension length, the average broken belt width, the average breaking distance and the crack density, and obtaining a sliding fracture quantitative representation result;
determining the fracture level according to the sliding fracture quantitative characterization result and preset conditions;
The quantization characterization module is specifically configured to:
and determining the fracture strain type segmentation according to the residual trend surface type, wherein the extrusion type corresponds to the extrusion segment, the tension type corresponds to the tension segment, and the translation type corresponds to the translation segment.
6. A skid break quantitative characterization and evaluation device, comprising: a memory, a processor and a computer program stored in the memory, the processor running the computer program to perform the method of any one of claims 1 to 4.
CN202010265385.XA 2020-04-07 2020-04-07 Quantitative characterization and evaluation method and device for sliding fracture Active CN113495294B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010265385.XA CN113495294B (en) 2020-04-07 2020-04-07 Quantitative characterization and evaluation method and device for sliding fracture

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010265385.XA CN113495294B (en) 2020-04-07 2020-04-07 Quantitative characterization and evaluation method and device for sliding fracture

Publications (2)

Publication Number Publication Date
CN113495294A CN113495294A (en) 2021-10-12
CN113495294B true CN113495294B (en) 2023-12-22

Family

ID=77995573

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010265385.XA Active CN113495294B (en) 2020-04-07 2020-04-07 Quantitative characterization and evaluation method and device for sliding fracture

Country Status (1)

Country Link
CN (1) CN113495294B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105510964A (en) * 2015-11-27 2016-04-20 中国石油大学(华东) Seismic recognition method of low-order strike-slip faults in complex structural areas
CN106569268A (en) * 2015-10-10 2017-04-19 中国石油化工股份有限公司 Method for quantitatively identifying type of fault structure
CN108107476A (en) * 2016-11-25 2018-06-01 中国石油化工股份有限公司 Broom shape strike-slip fault walks the computational methods of sliding displacement
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
CN110658556A (en) * 2019-10-24 2020-01-07 西南石油大学 Seismic technology combination method for identifying and evaluating broken zone of sliding fracture of carbonate rock
CN110850468A (en) * 2019-10-09 2020-02-28 中国石油大学(华东) Method for quantitatively representing pressurization and pressure release strength of bending slip fault

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106569268A (en) * 2015-10-10 2017-04-19 中国石油化工股份有限公司 Method for quantitatively identifying type of fault structure
CN105510964A (en) * 2015-11-27 2016-04-20 中国石油大学(华东) Seismic recognition method of low-order strike-slip faults in complex structural areas
CN108107476A (en) * 2016-11-25 2018-06-01 中国石油化工股份有限公司 Broom shape strike-slip fault walks the computational methods of sliding displacement
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
CN110850468A (en) * 2019-10-09 2020-02-28 中国石油大学(华东) Method for quantitatively representing pressurization and pressure release strength of bending slip fault
CN110658556A (en) * 2019-10-24 2020-01-07 西南石油大学 Seismic technology combination method for identifying and evaluating broken zone of sliding fracture of carbonate rock

Also Published As

Publication number Publication date
CN113495294A (en) 2021-10-12

Similar Documents

Publication Publication Date Title
US10359541B2 (en) Creating virtual production logging tool profiles for improved history matching
CN111458767B (en) Method and system for identifying lithology based on intersection graph method
JP7137889B2 (en) Prediction method of the lower limit of oil pool endowment depth of carbonatite in the deep-ultra-deep ocean facies
Hu et al. Water flooding flowing area identification for oil reservoirs based on the method of streamline clustering artificial intelligence
CN112147684A (en) Method and device for characterizing activity intensity of co-deposition fault
US10627542B2 (en) Hybrid approach to assisted history matching in large reservoirs
Yamamoto Quantification of uncertainty in ore-reserve estimation: Applications to Chapada copper deposit, State of Goiás, Brazil
CN113495294B (en) Quantitative characterization and evaluation method and device for sliding fracture
CN111025391B (en) Quantitative evaluation method for fault activity
CN114488297A (en) Fault identification method and device
CN112288243B (en) Method and device for evaluating associated metal resources in coal
CN110703328B (en) Overpressure interface identification method, device, equipment and storage medium
CN110569512B (en) Quantitative characterization method for development degree of permeability interlayer in shale
Hongwei et al. Fault identification technology of ant attribute optimization
CN108875163A (en) A kind of method and system for evaluating three-dimensional cracking network connectivty
CN114114390B (en) Oil and gas exploration method, device, equipment and storage medium
CN104462792A (en) Log data lithological stratum numerical value reduction method
CN111781661A (en) Method and device for predicting sedimentary microphase plane spread in small-well area
CN116757096B (en) Determination method, device, equipment and medium for deep, thin and compact sandstone gas reservoir development well
US20240110472A1 (en) Optimal drilling and fracturing sequences for placing numerous horizontal wells in tight reservoirs
US11808148B2 (en) Systems and methods for back-allocation of oil produced by waterflooding
CN117823041A (en) Geosteering model building method and device, electronic equipment and storage medium
CN117988714A (en) Method, device, equipment and medium for determining drilling track of horizontal well
CN112989257A (en) Method for measuring gas production rate of marine phase shale oil and gas reservoir
CN114117948A (en) Visual flow field characterization method based on data mining

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