Prediction method, system and equipment for width of fracture zone of thrust fault and medium

Abstract

The invention relates to a method, a system, equipment and a medium for predicting the width of a fracture zone of a thrust fault, and the method comprises the following steps: obtaining original data which comprise geological data, seismic data and logging data; processing the seismic data to obtain a fault displacement and a fault section inclination angle of a thrust fault; processing the logging data to obtain a brittleness index; processing the geological data to obtain a region coefficient; and based on the fault displacement, the fault section inclination angle of the thrust fault, the brittleness index and the region coefficient, calculating to obtain the thrust fault fracture zone width. According to the novel method for quantitatively calculating the width of the thrust fault fracture zone after the two stratum parameters of the brittleness index and the thrust fault section inclination angle are considered, the lithologic brittleness and the thrust fault section angle are considered, the calculation precision is improved, and a good effect is achieved in calculating the width of the thrust fault fracture zone.

Description

Technical Field

[0001] The present invention belongs to the technical field of seismic exploration, and in particular relates to a method, system, equipment and medium for predicting the width of a thrust fault fracture zone. Background Art

[0002] Faults are a crucial component in describing geological models and structures. Faults are often accompanied by fracture zones, which can serve as potential oil and gas migration pathways and reservoirs. The study of faults and fracture zones, particularly their width and extent, is crucial in petroleum exploration. Accurately describing faults and fault-fracture zones is a challenging task, requiring detailed geological data and fault-related information.

[0003] This invention involves the quantitative calculation of fracture zones in reverse faults. The geological meaning of a reverse fault is that regional strata are squeezed from both sides. When the squeezing exceeds the rock-fragmentation limit, the strata move toward each other, forming fracture surfaces and fracture zones near the fracture surfaces. Compared with normal faults and strike-slip faults, reverse faults are generally more likely to form large-scale fracture zones or fracture zones. The width and extent of these fracture zones are affected by a variety of factors, including the fault type, the magnitude of the fault activity, the rock types involved, and local geological conditions.

[0004] From a geological perspective, some key factors that usually affect the width and extent of fault fracture zones are:

[0005] (1) Fault type: Different types of faults exhibit different degrees of slip and motion. For example, a horizontal strike-slip fault (such as the San Andreas Fault in California) typically has a wider fracture zone than a normal fault or reverse fault.

[0006] (2) The magnitude of the fault activity: The magnitude of the fault activity, usually measured by the moment magnitude of the earthquake, affects the extent of the fault fracture zone. Larger earthquake events generally result in more significant damage and a wider affected area.

[0007] (3) Rock type and properties: The composition and properties of the rocks adjacent to the fault play a role in determining the extent of damage. Brittle rocks (such as crystalline rock formations or certain types of sedimentary rocks) are more susceptible to fracture and fault-related damage than ductile rocks (such as shale or clay).

[0008] (4) Stress regime: The stress regime in the region, including pre-existing stresses and tectonic forces, will affect the propagation and extent of fault damage.

[0009] (5) Geological heterogeneity: Changes in lithology, rock strength, and structural complexity within the stratum can affect the extent of the fault fracture zone. These changes may cause fault segments to bifurcate or propagate differently, thus affecting the overall fracture zone.

[0010] Accurately calculating the extent and width of a fault fracture zone requires detailed geological surveys, geophysical data, and a detailed analysis of the fault system. This typically involves field mapping, geotechnical investigations, seismic surveys, and other geological and geophysical techniques. It's important to note that without specific details about the fault and geological setting, it's impossible to quantitatively calculate or estimate the fault damage zone for a specific situation.

[0011] This patent involves the problem of how to quantitatively calculate the width of the reverse fault fracture zone. In the field of seismic oil and gas exploration, the common geological understanding is that the width of the fracture zone is proportional to the relative movement of the two plates of the reverse fault, that is, the fault throw. Based on this understanding, the simplest method to calculate the width of the fracture zone can be formed. fault =α*D fault

[0012] Among them D fault is the break distance, W fault is the crack zone width, and α is the area coefficient.

[0013] In some cases, this method has large calculation errors and does not take into account the influence of formation parameters such as rock brittleness and cross-section angle.

[0014] At present, the methods for predicting brittleness index mainly include prediction based on well logging and prediction based on seismic data. The main well logging prediction methods include the following two:

[0015] ① Array acoustic wave method. Grieser et al. used array acoustic wave logging curves to calculate the brittleness index, and used the longitudinal and shear wave velocities in the logging data to calculate the Poisson's ratio and Young's modulus curves. The Poisson's ratio and Young's modulus curves were normalized and outliers were removed. Finally, the average of the two was obtained to obtain the brittleness curve.

[0016] ② Mineral composition method: The mineral composition method obtains the brittleness index by calculating the ratio of the brittle mineral content or volume to the total content or total volume of rock skeleton minerals. However, the types and contents of brittle minerals in the strata of different target areas are different, so this method requires the establishment of a suitable brittle mineral model based on the rock mineral characteristics of the study area.

[0017] At present, through experimental analysis of mineral components, it is found that the content of quartz minerals is the main factor determining the size of brittleness, so it is proposed to calculate the proportion of quartz content in rock minerals to calculate the brittleness index, such as formula (1): Studies have found that in addition to quartz minerals in rocks, carbonate minerals also have an impact on brittleness and should also be classified as brittle minerals. Therefore, the brittleness index of shale is obtained by calculating the ratio of the total of the two, such as formula (2); the percentages of quartz, carbonate rock, and feldspar minerals are currently calculated to obtain the brittleness index, such as formula (3).

[0018] I b = Qtz / ( Qtz + Carb + Clay ) (1)

[0019] I b =(V Qtz + Carb ) / ( Qtz + Carb + Clay ) (2)

[0020] I b =(V Qtz + Carb + Feld ) / ( Qtz + Carb + Feld + clay ) (3)

[0021] In formula (1) to formula (3), I b is the brittleness index; V Qtz is the relative volume of quartz mineral; V Carb is the relative volume of carbonate minerals; V Feld is the relative volume of feldspar mineral; V Clay is the relative volume of clay.

[0022] The present invention proposes a new method for quantitatively calculating the width of a thrust fault fracture zone by taking into account two formation parameters, namely, a brittleness index and a thrust fault section dip angle. Summary of the Invention

[0023] The purpose of the present invention is to provide a method, system, equipment and medium for predicting the width of a thrust fault fracture zone in order to solve the above problems.

[0024] The present invention achieves the above-mentioned purpose through the following technical solutions:

[0025] A method for predicting the width of a thrust fault fracture zone comprises the following steps:

[0026] Acquiring original data, including geological data, seismic data, and well logging data;

[0027] The seismic data is processed to obtain the fault distance D fault and the fault section dip angle θ of the thrust fault;

[0028] The logging data is processed to obtain the brittleness index I b ;

[0029] Processing the geological data to obtain a regional coefficient α;

[0030] Based on the break distance D fault , the fault section inclination angle θ of the thrust fault, the brittleness index I b The thrust fault fracture zone width W is calculated using the regional coefficient α fault .

[0031] As a further optimization solution of the present invention, the seismic data is processed to obtain the fault distance D fault The specific process of the fault section dip angle θ of the thrust fault is as follows:

[0032] Performing data interpretation on the seismic data to obtain the structure of the thrust fault;

[0033] By analyzing the morphological characteristics of the thrust fault structure and using relevant algorithms, the fault distance D is calculated. fault and the fault section dip angle θ of the thrust fault.

[0034] As a further optimization solution of the present invention, the specific process of interpreting the seismic data to obtain the structure of the thrust fault is as follows:

[0035] Based on the regional geological background and geological knowledge, the seismic data are interpreted structurally and by horizon and fault, to obtain the structural morphological characteristics of the thrust fault and the geological structural morphology and spatial position of the area.

[0036] As a further optimization scheme of the present invention, by analyzing the morphological statistics of the thrust fault structure, the fault distance D is calculated based on the morphological characteristics of the thrust fault using a related algorithm. fault The specific process of the fault section dip angle θ of the thrust fault is as follows:

[0037] By analyzing the structural morphological characteristics of the thrust fault and performing time-depth conversion on the seismic data of the fault plane, sequentially reading the point data of the fault plane, determining the coherence radius value of the current point, and then performing a neighborhood search based on the coherence radius, performing plane fitting and parabola fitting on the many point data obtained from the search, and calculating the inclination, dip and strike of the current point based on the results of the plane fitting, thereby calculating the fault section dip angle θ of the thrust fault;

[0038] Find the intersection of the fault plane and the horizon. Ideally, the intersection is two curves distributed in space, which are the intersection lines of the upper and lower walls of the fault and the horizon. The fault distance D of the thrust fault can be solved based on the intersection data. fault .

[0039] As a further optimization solution of the present invention, the well logging data is processed to obtain the brittleness index I b The specific process is as follows:

[0040] Based on the well logging data, the lithologic composition of the stratum where the thrust fault is located is estimated, and the relative volume contents of quartz minerals, carbonate minerals, feldspar minerals, and clay are calculated. The brittleness index I of the thrust fault is calculated according to the following formula: b :

[0041] b =(V Qtz + Carb + Feld ) / ( Qtz + Carb + Feld + Clay );

[0042] Where, I b is the brittleness index; V Qtz is the relative volume of quartz mineral; V Carb is the relative volume of carbonate minerals; V Feld is the relative volume of feldspar mineral; V Clay is the relative volume of clay.

[0043] As a further optimization solution of the present invention, the specific process of processing the geological data to obtain the regional coefficient α is as follows:

[0044] The regional coefficient α is obtained based on the known fault throw in the thrust fault region and the estimated fracture zone width in the region.

[0045] As a further optimization solution of the present invention, based on the break distance D fault , the fault section inclination angle θ of the thrust fault, the brittleness index I b The thrust fault fracture zone width W is calculated using the regional coefficient α fault The specific process is as follows:

[0046] The thrust fault fracture zone width W is calculated according to the following formula: fault :

[0047] W fault =*I b * fault *sin(θ);

[0048] Where α is the area coefficient, I b is the brittleness index, θ is the thrust fault section inclination angle, 0°≤θ≤90°, D fault is the break distance, W fault is the width of the thrust fault zone.

[0049] A thrust fault fracture zone width prediction system, comprising:

[0050] A data acquisition module is used to acquire raw data, including geological data, seismic data, and well logging data;

[0051] Parameter calculation module, used to process the seismic data to obtain the fault distance D fault and the fault section inclination angle θ of the thrust fault; the logging data is processed to obtain the brittleness index I b ; Processing the geological data to obtain the regional coefficient α;

[0052] Width calculation module, for calculating the width of the fault , the fault section inclination angle θ of the thrust fault, the brittleness index I b The thrust fault fracture zone width W is calculated using the regional coefficient α fault .

[0053] An electronic device comprises a processor, a communication interface, a memory and a communication bus, wherein the processor, the communication interface and the memory communicate with each other via the communication bus;

[0054] Memory for storing computer programs;

[0055] The processor is used to implement a method for predicting the width of a thrust fault fracture zone when executing a program stored in the memory.

[0056] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements a method for predicting the width of a thrust fault fracture zone.

[0057] The beneficial effects of the present invention are:

[0058] In order to overcome the large calculation errors of the original method in certain cases and the failure to consider the influence of two formation parameters, lithologic brittleness and cross-section angle, the present invention proposes a new method for quantitatively calculating the width of the thrust fault fracture zone after considering two formation parameters, brittleness index and thrust fault cross-section angle. This method takes into account lithologic brittleness and reverse fault cross-section angle, improves the calculation accuracy, and achieves good results in calculating the width of the thrust fault fracture zone. BRIEF DESCRIPTION OF THE DRAWINGS

[0059] Figure 1 is a flow chart of the method of the present invention;

[0060] Figure 2 This is a flowchart of a method for predicting the width of a thrust fault fracture zone in an embodiment of the present invention;

[0061] Figure 3 Schematic diagram of the relationship between the width of the thrust fault fracture zone, the thrust fault angle, and the fault throw in an embodiment of the present invention;

[0062] Figure 4 This is a fracture characteristic map of the fault zone of Keshen 207 Well in the Keshen 2 gas reservoir in the Kuqa Depression in an embodiment of the present invention;

[0063] Figure 5 This is a seismic data morphology interpretation diagram of the structural morphology of the thrust fault in an embodiment of the present invention;

[0064] Figure 6 1 is a schematic diagram of calculation results of the width of the thrust fault fracture zone in an embodiment of the present invention;

[0065] Figure 7 is a schematic diagram of the lithologic components of the formation estimated from the well logging data in an embodiment of the present invention;

[0066] Figure 8 is a statistical diagram of the regional coefficient α of the thrust fault in the embodiment of the present invention;

[0067] Figure 9 is a system structure block diagram in an embodiment of the present invention;

[0068] Figure 10 It is a block diagram of the device structure in an embodiment of the present invention. DETAILED DESCRIPTION

[0069] The present application is further described in detail below in conjunction with the accompanying drawings. It is necessary to point out here that the following specific implementation methods are only used to further illustrate the present application and cannot be understood as limiting the scope of protection of the present application. Technical personnel in this field can make some non-essential improvements and adjustments to the present application based on the above application content.

[0070] The main idea of ​​this method is to first interpret the thrust fault using seismic data to interpret the thrust fault and obtain the structural morphology of the thrust fault. The fault throw Dfault is calculated based on the morphological characteristics of the thrust fault and the fault section dip angle θ (0° to 90°) of the thrust fault is obtained; the lithologic components of the thrust fault are estimated by obtaining well logging data from the work area where the thrust fault is located, and the relative volume contents of quartz mineral VQtz, carbonate mineral VCarb, feldspar mineral VFeld, and clay VClay are calculated. The brittleness index Ib of the thrust fault is calculated according to the calculation formula; then, the regional coefficient α is statistically calculated by considering the known fault throws and estimated fracture zone widths in the thrust fault area; finally, the width of the thrust fault fracture zone is calculated using the method proposed in this invention.

[0071] like Figure 1 As shown, a method for predicting the width of a thrust fault fracture zone includes the following steps:

[0072] Acquiring original data, including geological data, seismic data, and well logging data;

[0073] The seismic data is processed to obtain the fault distance D fault and the fault section dip angle θ of the thrust fault;

[0074] The logging data is processed to obtain the brittleness index I b ;

[0075] Processing the geological data to obtain a regional coefficient α;

[0076] Based on the break distance D fault , the fault section inclination angle θ of the thrust fault, the brittleness index I b The thrust fault fracture zone width W is calculated using the regional coefficient α fault .

[0077] Figure 2 As shown, in this embodiment, the steps are as follows:

[0078] (1) Obtain seismic data and interpret the data to obtain the structure of the thrust fault.

[0079] (2) By analyzing the morphological characteristics of the thrust fault and using relevant algorithms, the fault distance D is calculated based on the morphological characteristics of the thrust fault. fault And calculate the fault section inclination angle θ of the thrust fault.

[0080] (3) By collecting well logging data, the lithologic composition of the stratum where the thrust fault is located is estimated, and the V of the quartz mineral is calculated. Qtz 、Carbonate mineral V Carb , feldspar mineral V Feld , Clay V Clay The relative volume content of the thrust fault is calculated according to formula (3): b .

[0081] (4) The coefficient α of the region is calculated based on the known fault throw and estimated fracture zone width in the thrust fault region.

[0082] (5) Finally, the width of the thrust fault fracture zone is calculated according to the method for calculating the width of the thrust fault fracture zone proposed in the present invention.

[0083] Specifically, calculations of reservoir fracture zones in wells such as Keshen 207 and Dibei 5 in the Kuqa Depression of the Tarim Basin show that, combined with test and production trial data and imaging logging data, the reservoir fractures are relatively wide and contribute significantly to the communication between the edge water of the gas reservoir and the matrix pores inside the gas reservoir ( Figure 4 、 Figure 5 ).

[0084] The specific implementation plan is as follows:

[0085] A method for predicting the width of a thrust fault fracture zone provided by the present invention includes the following detailed steps:

[0086] Step 1: First, obtain the seismic data of a certain work area containing thrust faults, and perform structural interpretation, stratigraphic interpretation, and fault interpretation on the seismic data based on the geological background and geological knowledge of the area. Then, obtain the structural morphological characteristics of the thrust faults and the geological structural morphology and spatial position of the area, such as Figure 6 ;

[0087] Step 2: The relationship between the width of the thrust fault fracture zone and the thrust fault angle and fault throw is as follows: Figure 3 As shown in the figure, by analyzing the structural morphological characteristics of the thrust fault and performing time-depth conversion on the seismic data of the fault plane, the point data of the fault plane are sequentially read in, the coherence radius value of the current point is determined, and then a neighborhood search is performed based on the coherence radius. The many point data obtained from the search are subjected to plane fitting and parabola fitting. The dip, inclination and strike of the current point are obtained according to the results of the plane fitting, so as to calculate the angular attributes of the thrust fault. At the same time, the intersection of the fault plane and the horizon is found. Ideally, these intersections are two curves distributed in space, which are the intersection lines of the upper and lower walls of the fault and the horizon. Then, the fault throw attributes of the thrust fault are solved based on the intersection data.

[0088] Step 3: Estimate the lithologic composition of the thrust fault formation by collecting well logging data and calculate the quartz mineral V Qtz 、Carbonate mineral V Carb , feldspar mineral V Feld , Clay V Clay The relative volume content ( Figure 7 ), the brittleness index I of the thrust fault is calculated by formula (3): b ;

[0089] I b =(V Qtz +V Carb +V Feld ) / (V Qtz +V Carb +V Feld +V Clay )

[0090] Where, I b is the brittleness index; V Qtz is the relative volume of quartz mineral; V Carb is the relative volume of carbonate minerals; V Feld is the relative volume of feldspar mineral; V Clay is the relative volume of clay.

[0091] Step 4: Using the existing data, obtain a statistical diagram of known fault throws and estimated fracture zone widths in the thrust fault area ( Figure 8 ), and obtain the α of the regional coefficient;

[0092] Step 5: Finally, the width of the thrust fault fracture zone is calculated according to the formula (4) proposed in the present invention.

[0093] W fault =α*I b *D fault *sin(θ)

[0094] Where α is the area coefficient, I b is the brittleness index, θ is the dip angle of the thrust fault (0°~90°), D fault is the break distance, W fault is the width of the thrust fault zone.

[0095] It will be apparent to those skilled in the art that the present invention is not limited to the details of the exemplary embodiments described above and that the present invention can be implemented in other specific forms without departing from the spirit or essential characteristics of the present invention. Therefore, the embodiments should be considered in all respects as illustrative and non-restrictive, and the scope of the present invention is defined by the appended claims rather than the foregoing description.

[0096] like Figure 8 As shown, an embodiment of the present disclosure provides a system for predicting the width of a thrust fault fracture zone, comprising:

[0097] A data acquisition module is used to acquire raw data, including geological data, seismic data, and well logging data;

[0098] Parameter calculation module, used to process the seismic data to obtain the fault distance D fault and the fault section inclination angle θ of the thrust fault; the logging data is processed to obtain the brittleness index I b ; Processing the geological data to obtain the regional coefficient α;

[0099] Width calculation module, for calculating the width of the fault , the fault section inclination angle θ of the thrust fault, the brittleness index I b The thrust fault fracture zone width W is calculated using the regional coefficient α fault .

[0100] The implementation process of the functions and effects of each module in the above system is specifically described in the implementation process of the corresponding steps in the above method, which will not be repeated here.

[0101] For the system embodiment, since it basically corresponds to the method embodiment, the relevant parts can be referred to the partial description of the method embodiment. The system embodiment described above is only illustrative, wherein the modules described as separate components may or may not be physically separated, and the components displayed as modules may or may not be physical modules, that is, they may be located in one place, or they may be distributed on multiple network modules. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of the present invention. Those of ordinary skill in the art can understand and implement it without paying any creative work.

[0102] In the above embodiment, any number of all modules can be combined into one module for implementation, or any one of the modules can be split into multiple modules. Alternatively, at least some of the functions of one or more of these modules can be combined with at least some of the functions of other modules and implemented in one module. At least one of all modules can be at least partially implemented as a hardware circuit, such as a field programmable gate array (FPGA), a programmable logic array (PLA), a system on a chip, a system on a substrate, a system on a package, an application specific integrated circuit (ASIC), or can be implemented by hardware or firmware such as any other reasonable way of integrating or packaging the circuit, or implemented in any one of the three implementation modes of software, hardware and firmware or in a suitable combination of any of them. Alternatively, at least one of all modules can be at least partially implemented as a computer program module, which can perform the corresponding function when the computer program module is run.

[0103] See also Figure 9 The electronic device provided by an embodiment of the present disclosure includes a processor 1110, a communication interface 1120, a memory 1130 and a communication bus 1140, wherein the processor 1110, the communication interface 1120, and the memory 1130 communicate with each other through the communication bus 1140;

[0104] Memory 1130, for storing computer programs;

[0105] The processor 1110 is configured to implement the following method for predicting the width of a thrust fault fracture zone when executing the program stored in the memory 1130 .

[0106] The communication bus 1140 can be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus. The communication bus 1140 can be divided into an address bus, a data bus, a control bus, etc. For ease of illustration, the figure shows only one thick line, but this does not mean that there is only one bus or only one type of bus.

[0107] The communication interface 1120 is used for communication between the electronic device and other devices.

[0108] The memory 1130 may include a random access memory (RAM) or a non-volatile memory, such as at least one disk storage. Alternatively, the memory 1130 may be at least one storage device located away from the processor 1110.

[0109] The above-mentioned processor 1110 can be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it can also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.

[0110] The embodiments of the present disclosure further provide a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for predicting the width of a thrust fault fracture zone as described above.

[0111] The computer-readable storage medium may be included in the device / apparatus described in the above embodiments, or may exist independently and not be incorporated into the device / apparatus. The computer-readable storage medium carries one or more programs that, when executed, implement the method for predicting the width of a thrust fault fracture zone according to the embodiments of the present disclosure.

[0112] According to an embodiment of the present disclosure, a computer-readable storage medium may be a non-volatile computer-readable storage medium, such as, but not limited to, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination thereof. In the present disclosure, a computer-readable storage medium may be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, apparatus, or device.

[0113] The above-described embodiments merely illustrate several implementations of the present invention. While the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the present invention. It should be noted that a person skilled in the art would be able to make numerous variations and improvements without departing from the spirit of the present invention, and all such variations and improvements fall within the scope of protection of the present invention.

Claims

1. A method for predicting the width of a thrust fault fracture zone, characterized in that: The following steps are involved: Acquiring original data, including geological data, seismic data, and well logging data; The seismic data is processed to obtain the fault distance D fault and the fault section dip angle θ of the thrust fault; The logging data is processed to obtain the brittleness index I b ; Processing the geological data to obtain a regional coefficient α; Based on the break distance D fault , the fault section inclination angle θ of the thrust fault, the brittleness index I b The thrust fault fracture zone width W is calculated using the regional coefficient α fault .

2. The method for predicting the width of a thrust fault fracture zone according to claim 1, wherein: The seismic data is processed to obtain the fault distance D fault The specific process of the fault section dip angle θ of the thrust fault is as follows: Performing data interpretation on the seismic data to obtain the structure of the thrust fault; By analyzing the morphological characteristics of the thrust fault structure and using relevant algorithms, the fault distance D is calculated. fault and the fault section dip angle θ of the thrust fault.

3. The method for predicting the width of a thrust fault fracture zone according to claim 2, wherein: The specific process of interpreting the seismic data to obtain the structure of the thrust fault is as follows: Based on the regional geological background and geological knowledge, the seismic data are interpreted structurally and by horizon and fault, to obtain the structural morphological characteristics of the thrust fault and the geological structural morphology and spatial position of the area.

4. The method for predicting the width of a thrust fault fracture zone according to claim 3, wherein: By analyzing the morphological characteristics of the thrust fault structure and using relevant algorithms, the fault distance D is calculated. fault The specific process of the fault section dip angle θ of the thrust fault is as follows: By analyzing the structural morphological characteristics of the thrust fault and performing time-depth conversion on the seismic data of the fault plane, sequentially reading the point data of the fault plane, determining the coherence radius value of the current point, and then performing a neighborhood search based on the coherence radius, performing plane fitting and parabola fitting on the many point data obtained from the search, and calculating the inclination, dip and strike of the current point based on the results of the plane fitting, thereby calculating the fault section dip angle θ of the thrust fault; Find the intersection of the fault plane and the horizon. Ideally, the intersection is two curves distributed in space, which are the intersection lines of the upper and lower walls of the fault and the horizon. The fault distance D of the thrust fault can be solved based on the intersection data. fault .

5. The method for predicting the width of a thrust fault fracture zone according to claim 1, wherein: The logging data is processed to obtain the brittleness index I b The specific process is as follows: Based on the well logging data, the lithologic composition of the stratum where the thrust fault is located is estimated, and the relative volume contents of quartz minerals, carbonate minerals, feldspar minerals, and clay are calculated. The brittleness index I of the thrust fault is calculated according to the following formula: b : I b =(V Qtz +V Carb +V Feld ) / (V Qtz +V Carb +V Feld +V Clay ); Where, I b is the brittleness index; V Qtz is the relative volume of quartz mineral; V Carb is the relative volume of carbonate minerals; V Feld is the relative volume of feldspar mineral; V Clay is the relative volume of clay.

6. The method for predicting the width of a thrust fault fracture zone according to claim 1, wherein: The specific process of processing the geological data to obtain the regional coefficient α is as follows: The regional coefficient α is obtained based on the known fault throw in the thrust fault region and the estimated fracture zone width in the region.

7. The method for predicting the width of a thrust fault fracture zone according to claim 1, wherein: Based on the break distance D fault , the fault section inclination angle θ of the thrust fault, the brittleness index I b The thrust fault fracture zone width W is calculated using the regional coefficient α fault The specific process is as follows: The thrust fault fracture zone width W is calculated according to the following formula: fault : W fault =α*I b *D fault *sin(θ); Where α is the area coefficient, I b is the brittleness index, θ is the thrust fault section inclination angle, 0°≤θ≤90°, D fault is the break distance, W fault is the width of the thrust fault zone.

8. A system for predicting the width of a thrust fault fracture zone, characterized in that: include: A data acquisition module is used to acquire raw data, including geological data, seismic data, and well logging data; Parameter calculation module, used to process the seismic data to obtain the fault distance D fault and the fault section inclination angle θ of the thrust fault; the logging data is processed to obtain the brittleness index I b ; Processing the geological data to obtain a regional coefficient α; Width calculation module, for calculating the width of the fault , the fault section inclination angle θ of the thrust fault, the brittleness index I b The thrust fault fracture zone width W is calculated using the regional coefficient α fault .

9. An electronic device, characterized in that: The processor, the communication interface, the memory and the communication bus are connected to each other via the communication bus. Memory for storing computer programs; The processor is configured to implement the method for predicting the width of a thrust fault fracture zone according to any one of claims 1 to 7 when executing the program stored in the memory.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the method for predicting the width of a thrust fault fracture zone according to any one of claims 1 to 7 is implemented.

Images