Earthquake assessment method and device for fault stability

Abstract

The embodiment of the invention provides an earthquake assessment method and device for fault stability, and belongs to the technical field of oil and gas seismic exploration. The method comprises the following steps: acquiring seismic basic data of a drilling section; fracture interpretation is conducted based on actual drilling data, and the fracture development position of the drilling section is recognized; based on the seismic data, seismic attributes capable of predicting fractures are extracted; determining a sensitive attribute body, a fracture quantitative characterization result, a crustal stress prediction result and a crustal stress direction prediction result; on the basis of the Moire-Coulomb theory and on the basis of the fracture quantitative characterization result, the shear stress and the effective normal stress of the fracture surface are calculated in combination with the ground stress prediction result and the ground stress direction prediction result, and a fracture stability prediction plane graph is obtained on the basis of the shear stress and the effective normal stress of the fracture surface. According to the method, fracture surface shear stress and effective normal stress are calculated according to earthquake fracture and crustal stress field prediction results, so that the stable state of the fracture surface is evaluated, and the stability of the fault can be quantitatively evaluated.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas seismic exploration, and in particular to a seismic evaluation method for fault stability, a seismic evaluation device for fault stability, a machine-readable storage medium and a processor. Background Art

[0002] In unconventional exploration and development, casing change is an important factor that restricts the development of scale benefits. In the current shale gas exploration and development, casing change is common in the drilling and fracturing stages. There are many factors that cause casing change, such as temperature, asymmetric fracturing, cementing quality, and fault slip. At present, the research on casing change caused by fault slip focuses on single wellbore analysis, including cementing stage, casing load, cement sheath, and ground stress analysis. By establishing a multi-stage fracturing numerical model, the ground stress changes near the wellbore are calculated. There is no literature mentioning the study of the mechanical stability of the fracture itself. The distribution of underground stress, the spatial angle relationship between the maximum principal stress and the fracture surface, and the relationship between the friction force and shear stress of the fracture surface are all important factors affecting the mechanical stability of the fracture.

[0003] Existing geophysical prediction methods include: (1) prediction of faults of different scales, using post-stack data to extract coherence, curvature, texture, symmetry, ant volume, maximum likelihood attributes, etc. for qualitative prediction; (2) seismic prediction of ground stress, using pre-stack inversion of P-wave impedance, S-wave impedance, density and P-S wave velocity ratio, and then calculating rock mechanical parameters such as Young's modulus, Poisson's ratio, brittleness index, etc.

[0004] However, existing geophysical prediction methods cannot quantitatively assess the stability of faults, making it difficult to guide well deployment and avoid drilling risks. Summary of the invention

[0005] The purpose of the embodiments of the present application is to provide a seismic assessment method and device for fault stability, which can quantitatively assess the stability of the fault to guide well site deployment and avoid drilling risks.

[0006] In order to achieve the above-mentioned object, the present application provides, in a first aspect, a method for seismic assessment of fault stability, the method comprising: Acquire basic seismic data of the drilling section; wherein the basic seismic data of the drilling section includes seismic data and actual drilling data; Performing fracture interpretation based on the actual drilling data to identify fracture development locations in the drilling section; Based on the seismic data, extracting seismic attributes capable of predicting rupture; Calibrate the fracture development position of the drilling section and the seismic attributes to determine a sensitive attribute body, and determine a fracture quantitative characterization result based on the sensitive attribute body; Performing elastic parameter inversion based on the seismic data to obtain an elastic parameter body, and determining a rock mechanics parameter body based on the elastic parameter body; Determining a ground stress prediction result based on the elastic parameter body and the rock mechanics parameter body; Based on the thin plate theory and using the rock mechanics parameter body, geological, mechanical and mathematical models are established, and the three-dimensional finite element difference simulation method is used to determine the prediction results of the ground stress direction; Based on the Mohr-Coulomb theory, on the basis of the quantitative characterization results of the fracture, combined with the ground stress prediction results and the ground stress direction prediction results, the fracture surface shear stress and effective normal stress are calculated, and the fracture stability prediction plane diagram is obtained based on the fracture surface shear stress and effective normal stress.

[0007] Optionally, performing fracture interpretation based on the actual drilling data to identify fracture development locations in the drilling section includes: The location of fault development in the drilling section is interpreted through the actual drilling engineering recovery map, FMI imaging logging data, conventional logging data, horizontal section logging data and microseismic monitoring data.

[0008] Optionally, the seismic data includes a time domain prestack migration gather; and performing elastic parameter inversion based on the seismic data to obtain an elastic parameter volume includes: The elastic parameter inversion is completed using the time domain prestack migration gathers to obtain the elastic parameter volume.

[0009] Optionally, the geostress prediction result includes overlying formation pressure; and the method includes: Based on the elastic parameter body and the rock mechanics parameter body, a normal compaction curve of density is calculated using a seismic travel time equation; Create a density curve using the normal compaction curve of density and the actual measured density curve; The overburden pressure is calculated using the integral of the density curve.

[0010] Optionally, the geostress prediction result also includes formation pore pressure; and the method includes: Based on the elastic parameter body and the rock mechanics parameter body, the formation pore pressure is calculated using Eaton's empirical formula.

[0011] Optionally, the geostress prediction result also includes horizontal principal stress; and the method includes: Based on the elastic parameter body and the rock mechanics parameter body, the horizontal principal stress is calculated using the Huang model.

[0012] Optionally, a fracture stability prediction plane diagram is obtained based on the fracture surface shear stress and effective normal stress, including: The magnitudes of the shear stress and the effective normal stress of the fracture surface are compared to obtain a comparison result; and a fracture stability prediction plane diagram is drawn according to the comparison result.

[0013] A second aspect of the present application provides a seismic assessment device for fault stability, the device comprising: A basic data acquisition module is used to acquire basic seismic data of the drilling section; wherein the basic seismic data of the drilling section includes seismic data and actual drilling data; A fracture interpretation module, used to interpret fractures based on the actual drilling data and identify fracture development locations in the drilling section; A seismic attribute extraction module, used for extracting seismic attributes capable of predicting fractures based on the seismic data; A calibration module, used to calibrate the fracture development position of the drilling section and the seismic attribute, determine a sensitive attribute body, and determine a fracture quantitative characterization result based on the sensitive attribute body; An inversion module, used for performing elastic parameter inversion based on the seismic data to obtain an elastic parameter body, and determining a rock mechanics parameter body according to the elastic parameter body; A geostress prediction module, used for determining a geostress prediction result based on the elastic parameter body and the rock mechanics parameter body; The module for predicting the direction of geostress is used to establish geological, mechanical and mathematical models based on the thin plate theory and the rock mechanics parameter body, and to determine the prediction results of the direction of geostress using the three-dimensional finite element difference simulation method; The fracture stability prediction module is used to calculate the fracture surface shear stress and effective normal stress based on the Mohr-Coulomb theory and the fracture quantitative characterization results, combined with the ground stress prediction results and the ground stress direction prediction results, and obtain the fracture stability prediction plane diagram based on the fracture surface shear stress and the effective normal stress.

[0014] A third aspect of the present application provides a processor configured to execute the above-mentioned seismic assessment method of fault stability.

[0015] A fourth aspect of the present application provides a machine-readable storage medium having instructions stored thereon, which, when executed by a processor, configure the processor to execute the above-mentioned seismic assessment method for fault stability.

[0016] Compared with the prior art, the above technical solution of the present invention has the following beneficial effects: The present application provides a seismic assessment method and device for fault stability. The method uses the prediction results of earthquake fractures and geostress fields to calculate the shear stress and effective normal stress of the fracture surface to evaluate the stability of the fracture surface. The method can quantitatively assess the stability of the fault to guide well site deployment and avoid drilling risks.

[0017] Other features and advantages of the embodiments of the present application will be described in detail in the subsequent specific implementation section. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] The accompanying drawings are used to provide a further understanding of the embodiments of the present invention and constitute a part of the specification. Together with the following specific implementations, they are used to explain the embodiments of the present invention, but do not constitute a limitation on the embodiments of the present invention. In the accompanying drawings: Figure 1 The overall process diagram of the seismic evaluation method of fault stability according to the embodiment of the present application is schematically shown; Figure 2 A schematic diagram of a specific process of a seismic evaluation method for fault stability according to an embodiment of the present application is schematically shown; Figure 3 A schematic diagram of the predicted plane distribution of overburden formation pressure according to an embodiment of the present application is shown schematically; Figure 4 A schematic diagram of the horizontal principal stress difference superimposed on the maximum horizontal principal stress direction according to an embodiment of the present application is shown; Figure 5 Schematically shows a fracture space channel and a quantitative characterization schematic diagram according to an embodiment of the present application; Figure 6 A schematic diagram of a fracture-slip risk prediction plan view according to an embodiment of the present application is shown; Figure 7 The structural block diagram of a seismic assessment device for fault stability according to an embodiment of the present application is schematically shown. DETAILED DESCRIPTION

[0019] In order to make the purpose, technical scheme and advantages of the embodiments of the present application clearer, the technical scheme in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. It should be understood that the specific implementation methods described herein are only used to illustrate and explain the embodiments of the present application, and are not used to limit the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of this application.

[0020] Figure 1 The overall process diagram of the seismic assessment method of fault stability according to an embodiment of the present application is schematically shown. Figure 2 The specific flow chart of the seismic evaluation method of fault stability according to the embodiment of the present application is schematically shown. Figure 1-2As shown, in one embodiment of the present application, a seismic assessment method for fault stability is provided to quantitatively assess the stability of the fault and solve the engineering prediction problem of avoiding high-slip risk fractures during the drilling process, including the following steps: Step 110, obtaining basic seismic data of the drilling section; wherein the basic seismic data of the drilling section includes seismic data and actual drilling data.

[0021] Specifically, it is necessary to first obtain basic data for seismic assessment of fracture mechanics stability. Seismic data include pre-stack migration gathers in the time domain and full stack migration data in the time domain. Actual drilling data include drilling coordinates, drilling layers, drilling displays, drilling well trajectories, logging data from the surface to the underground target layer (including but not limited to longitudinal waves, density, shear waves, natural gamma, total hydrocarbon curves, gamma while drilling, etc.), element logging curves, actual drilling data of geostress (maximum horizontal principal stress and direction, minimum horizontal principal stress, overburden pressure, pore pressure, etc.), FMI imaging logging data, actual drilling engineering stratum recovery profiles, and microseismic monitoring data.

[0022] Step 120, performing fracture interpretation based on the actual drilling data to identify fracture development locations in the drilling section.

[0023] Specifically, the fracture development locations in the drilling sections (including the fracture development locations in the vertical and horizontal sections of the drilling) are identified through actual drilling engineering formation recovery profiles, FMI imaging logging data, conventional logging data (natural gamma, deep and shallow resistivity, acoustic time difference, density and neutron), horizontal section logging (total hydrocarbon curve, gamma while drilling) and microseismic monitoring results.

[0024] Step 130: extracting seismic attributes capable of predicting fractures based on the seismic data.

[0025] Specifically, based on the pre-stack migration gather in the time domain, the drilling trajectory in the depth domain is converted to the time domain by extracting seismic wavelets and making synthetic seismic records, so as to facilitate matching with the full stack migration data volume in the time domain. In addition, faults with more than half a phase mismatch are manually interpreted on the time domain seismic data. A variety of seismic attributes that can currently predict faults are extracted, including structural attributes (similarity, curvature, texture, maximum likelihood body, etc.), geometric attributes (body curvature, dip, azimuth, etc.), instantaneous attributes (amplitude, frequency, phase, envelope, energy, ant body anomalies, etc.), and fault prediction based on artificial intelligence.

[0026] Step 140, calibrating the fracture development position of the drilling section and the seismic attributes, determining a sensitive attribute body, and determining a fracture quantitative characterization result based on the sensitive attribute body.

[0027] Figure 5The schematic diagram of the fracture space channel and quantitative characterization according to the embodiment of the present application is schematically shown. Specifically, the fracture development position of the drilling section is calibrated with the seismic attributes, and 3-4 sensitive attribute bodies (usually including but not limited to body curvature, ant body, artificial intelligence fracture prediction) of possible fractures are selected. The fracture prediction attribute body (dip, azimuth, similarity or artificial intelligence fracture prediction or curvature) is selected, and the attribute values ​​are extracted along the main target layer segment. The extracted attribute values ​​(dip, similarity or curvature) are intersected, and points with high dip, low similarity or high body curvature values ​​are filtered out. The collection of these points in space is the fracture channel. The positions and characteristic values ​​(dip, dip) of these points are determined, and each fracture piece in the fracture channel can be simulated based on these values, such as Figure 5 shown.

[0028] Step 150, performing elastic parameter inversion based on the seismic data to obtain an elastic parameter body, and determining a rock mechanics parameter body based on the elastic parameter body.

[0029] Specifically, the elastic parameter inversion is completed using the time domain prestack migration gathers to obtain elastic parameter bodies such as longitudinal waves, shear waves, and density. Based on these elastic parameter bodies, rock mechanics parameter bodies such as Young's modulus, Poisson's ratio, and Lame constant are calculated.

[0030] Step 160: Determine a ground stress prediction result based on the elastic parameter body and the rock mechanics parameter body.

[0031] Figure 3 The schematic diagram of the planar distribution of the overburden formation pressure prediction according to the embodiment of the present application is schematically shown. Specifically, the above-mentioned elastic parameters and rock mechanics parameter bodies are used, and the relationship between porosity and depth (based on the Athy equation of seismic travel time) is used to calculate the normal compaction curve of the acoustic wave and density curve; the resistivity or gamma curve is used to calculate the mudstone indicator trend line and the clay content; for overpressure formations, a normal density curve is created using the normal compaction curve of density and the actual measured density curve; the overburden formation pressure is calculated using the integral of the density curve; the formation pore pressure is calculated using Eaton's empirical formula; and the horizontal principal stress is calculated using the Huang model. Based on the above process, the overburden formation pressure body, the formation pore pressure body and the horizontal stress body can be obtained.

[0032] Step 170, based on the thin plate theory and using the rock mechanics parameter body, geological, mechanical and mathematical models are established, and the three-dimensional finite element difference simulation method is used to determine the prediction result of the ground stress direction.

[0033] Figure 4The schematic diagram of the horizontal principal stress difference superimposed on the maximum horizontal principal stress direction according to the embodiment of the present application is schematically shown. Specifically, mainly based on the thin plate theory, using the rock mechanics parameter body, the geological, mechanical and mathematical models are established, and the maximum horizontal principal stress direction is simulated by using the three-dimensional finite element difference simulation method. The result is calibrated with the measured maximum principal stress direction, and the model parameters are continuously optimized and adjusted according to the calibration result until the consistency with the actual maximum principal stress direction is greater than 70%, such as Figure 4 shown.

[0034] Step 180, based on the Mohr-Coulomb theory, on the basis of the quantitative characterization results of the fracture, combined with the ground stress prediction results and the ground stress direction prediction results, the fracture surface shear stress and the effective normal stress are calculated, and the fracture stability prediction plane diagram is obtained based on the fracture surface shear stress and the effective normal stress.

[0035] Figure 6 The fracture slip risk prediction plane diagram according to the embodiment of the present application is schematically shown. Specifically, guided by the Mohr-Coulomb criterion (a classical theory of fracture surface sliding in mechanics), the results obtained by the above process, i.e., the maximum horizontal principal stress and direction, the minimum horizontal principal stress, and the cross-section inclination, are used to calculate the fracture surface shear stress plane distribution diagram and the effective normal stress plane distribution diagram, and the shear stress and friction force (effective normal stress) are compared to obtain the fracture stability prediction plane diagram, as shown in FIG. Figure 6 shown.

[0036] Figure 2 FIG. 1 is a flow chart of a method for seismic assessment of fault stability in one embodiment. Figure 2 The steps in the flowchart are shown in sequence as indicated by the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order restriction for the execution of these steps, and these steps can be executed in other orders. Moreover, Figure 2 At least part of the steps may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least part of the sub-steps or stages of other steps.

[0037] The method described in the embodiment of the present application borrows the mature seismic inversion technology, seismic attribute extraction technology, rock mechanics Eaton formation pore pressure empirical formula, three-dimensional finite element differential simulation method, ground stress calculation based on Huang's model, and the classical theory of fracture surface sliding dynamics Mohr-Coulomb criterion, and proposes an original technology for fracture surface seismic assessment. Mohr-Coulomb is often used to assess the sliding of wellbore fracture surfaces, and can only predict the sliding risk for drilling locations in the region. The embodiment of the present application takes advantage of the continuous spatial distribution of seismic data, and can predict the sliding risk of regional fracture surfaces, assess the sliding risk of different positions of the fracture surface, and assess the risk of different stress locations. Compared with the method of using digital random modeling to assess fracture risk in the later stage, the data between wellbores are all from deterministic seismic data, rather than the mathematical interpolation method of random space modeling, the data calculation amount is small, and the prediction accuracy and stability are higher.

[0038] The method described in the embodiment of this application has been effectively used in the exploration and development of deep shale gas in the southern Sichuan Basin. The method is used to evaluate fracture stability, which is consistent with actual drilling. It has currently supported the deep shale gas optimization pilot test plan in Luzhou area. The number of well positions in Lu 203 well area has increased by 16%, and more than 20 wells are being drilled. The drilling adjustment plan has increased the utilization rate of underground resources and effectively guided production, significantly reduced the probability of drilling casing changes, saved drilling time, and improved economic benefits.

[0039] In one embodiment, Figure 7 As shown, a seismic assessment device for fault stability is provided, the device comprising: The basic data acquisition module 210 is used to acquire the basic seismic data of the drilling section; wherein the basic seismic data of the drilling section includes seismic data and actual drilling data; A fracture interpretation module 220 is used to interpret fractures based on the actual drilling data and identify fracture development locations in the drilling section; A seismic attribute extraction module 230, for extracting seismic attributes capable of predicting fractures based on the seismic data; The calibration module 240 is used to calibrate the fracture development position of the drilling section and the seismic attribute, determine the sensitive attribute body, and determine the fracture quantitative characterization result according to the sensitive attribute body; An inversion module 250 is used to perform elastic parameter inversion based on the seismic data to obtain an elastic parameter body, and determine a rock mechanics parameter body based on the elastic parameter body; A geostress prediction module 260, configured to determine a geostress prediction result based on the elastic parameter body and the rock mechanics parameter body; The ground stress direction prediction module 270 is used to establish geological, mechanical and mathematical models based on the thin plate theory and the rock mechanics parameter body, and to determine the ground stress direction prediction result by using the three-dimensional finite element difference simulation method; The fracture stability prediction module 280 is used to calculate the fracture surface shear stress and effective normal stress based on the Mohr-Coulomb theory and the fracture quantitative characterization results, combined with the ground stress prediction results and the ground stress direction prediction results, and obtain the fracture stability prediction plane diagram based on the fracture surface shear stress and the effective normal stress.

[0040] The seismic assessment device for fault stability includes a processor and a memory. The basic data acquisition module, fracture interpretation module, seismic attribute extraction module, calibration module, inversion module, geostress prediction module, geostress direction prediction module and fracture stability prediction module are all stored in the memory as program units, and the processor executes the program units stored in the memory to realize corresponding functions.

[0041] The processor includes a kernel, and the kernel calls the corresponding program unit from the memory. One or more kernels can be set, and the seismic assessment method of fault stability is realized by adjusting kernel parameters.

[0042] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM, and the memory includes at least one memory chip.

[0043] An embodiment of the present invention provides a storage medium on which a program is stored. When the program is executed by a processor, the seismic assessment method of fault stability is implemented.

[0044] An embodiment of the present invention provides a processor, which is used to run a program, wherein the seismic assessment method of fault stability is executed when the program is run.

[0045] An embodiment of the present invention provides a device, the device including a processor, a memory, and a program stored in the memory and executable on the processor, and when the processor executes the program, the following steps are implemented: Step 110, obtaining basic seismic data of the drilling section; wherein the basic seismic data of the drilling section includes seismic data and actual drilling data.

[0046] Step 120, performing fracture interpretation based on the actual drilling data to identify fracture development locations in the drilling section.

[0047] Step 130: extracting seismic attributes capable of predicting fractures based on the seismic data.

[0048] Step 140, calibrating the fracture development position of the drilling section and the seismic attributes, determining a sensitive attribute body, and determining a fracture quantitative characterization result based on the sensitive attribute body.

[0049] Step 150, performing elastic parameter inversion based on the seismic data to obtain an elastic parameter body, and determining a rock mechanics parameter body based on the elastic parameter body.

[0050] Step 160: Determine a ground stress prediction result based on the elastic parameter body and the rock mechanics parameter body.

[0051] Step 170, based on the thin plate theory and using the rock mechanics parameter body, geological, mechanical and mathematical models are established, and the three-dimensional finite element difference simulation method is used to determine the prediction result of the ground stress direction.

[0052] Step 180, based on the Mohr-Coulomb theory, on the basis of the quantitative characterization results of the fracture, combined with the ground stress prediction results and the ground stress direction prediction results, the fracture surface shear stress and the effective normal stress are calculated, and the fracture stability prediction plane diagram is obtained based on the fracture surface shear stress and the effective normal stress.

[0053] Optionally, fracture interpretation is performed based on the actual drilling data to identify the fracture development location of the drilling section, including: interpreting the fracture development location of the drilling section through the actual drilling engineering recovery map of the drilling section, FMI imaging logging data, conventional logging data, horizontal section logging data and microseismic monitoring data.

[0054] Optionally, the seismic data includes a time-domain prestack migration gather; and performing elastic parameter inversion based on the seismic data to obtain an elastic parameter body includes: performing elastic parameter inversion using a time-domain prestack migration gather to obtain an elastic parameter body.

[0055] Optionally, the geostress prediction result includes overlying formation pressure; and the method includes: Based on the elastic parameter body and the rock mechanics parameter body, a normal compaction curve of density is calculated using a seismic travel time equation; Create a density curve using the normal compaction curve of density and the actual measured density curve; The overburden pressure is calculated using the integral of the density curve.

[0056] Optionally, the geostress prediction result also includes formation pore pressure; and the method includes: Based on the elastic parameter body and the rock mechanics parameter body, the formation pore pressure is calculated using Eaton's empirical formula.

[0057] Optionally, the geostress prediction result also includes horizontal principal stress; the method includes: calculating the horizontal principal stress using the Huang model based on the elastic parameter body and the rock mechanics parameter body.

[0058] Optionally, obtaining a fracture stability prediction plane diagram based on the fracture surface shear stress and the effective normal stress includes: comparing the fracture surface shear stress and the effective normal stress to obtain a comparison result; and drawing the fracture stability prediction plane diagram according to the comparison result.

[0059] The method described in the embodiment of the present application can guide drilling to avoid areas with high risk of fault sliding, improve the guiding role of seismic results in shale gas exploration and development, support overall platform deployment and optimization, and improve the utilization rate of underground resources.

[0060] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present application may adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.

[0061] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0062] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0063] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0064] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0065] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0066] Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

[0067] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

[0068] The above are only embodiments of the present application and are not intended to limit the present application. For those skilled in the art, the present application may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included within the scope of the claims of the present application.

Claims

1. A method for seismic assessment of fault stability, characterized in that: The method comprises: Acquire basic seismic data of the drilling section; wherein the basic seismic data of the drilling section includes seismic data and actual drilling data; Performing fracture interpretation based on the actual drilling data to identify fracture development locations in the drilling section; Based on the seismic data, extracting seismic attributes capable of predicting rupture; Calibrate the fracture development position of the drilling section and the seismic attributes to determine a sensitive attribute body, and determine a fracture quantitative characterization result based on the sensitive attribute body; Performing elastic parameter inversion based on the seismic data to obtain an elastic parameter body, and determining a rock mechanics parameter body based on the elastic parameter body; Determining a ground stress prediction result based on the elastic parameter body and the rock mechanics parameter body; Based on the thin plate theory and using the rock mechanics parameter body, geological, mechanical and mathematical models are established, and the three-dimensional finite element difference simulation method is used to determine the prediction results of the ground stress direction; Based on the Mohr-Coulomb theory, on the basis of the quantitative characterization results of the fracture, combined with the ground stress prediction results and the ground stress direction prediction results, the fracture surface shear stress and effective normal stress are calculated, and the fracture stability prediction plane diagram is obtained based on the fracture surface shear stress and effective normal stress.

2. The seismic assessment method of fault stability according to claim 1, characterized in that: Based on the actual drilling data, the fracture interpretation is performed to identify the fracture development location of the drilling section, including: The location of fault development in the drilling section is interpreted through the actual drilling engineering recovery map, FMI imaging logging data, conventional logging data, horizontal section logging data and microseismic monitoring data.

3. The seismic assessment method of fault stability according to claim 1, characterized in that: The seismic data includes a time domain pre-stack migration gather; performing elastic parameter inversion based on the seismic data to obtain an elastic parameter body includes: The elastic parameter inversion is completed using the time domain prestack migration gathers to obtain the elastic parameter volume.

4. The seismic assessment method of fault stability according to claim 1, characterized in that: The ground stress prediction result includes the overlying formation pressure; the method includes: Based on the elastic parameter body and the rock mechanics parameter body, a normal compaction curve of density is calculated using a seismic travel time equation; Create a density curve using the normal compaction curve of density and the actual measured density curve; The overburden pressure is calculated using the integral of the density curve.

5. The seismic evaluation method of fault stability according to claim 4, characterized in that: The geostress prediction result also includes formation pore pressure; the method includes: Based on the elastic parameter body and the rock mechanics parameter body, the formation pore pressure is calculated using Eaton's empirical formula.

6. The seismic assessment method of fault stability according to claim 5, characterized in that: The geostress prediction result also includes horizontal principal stress; the method includes: Based on the elastic parameter body and the rock mechanics parameter body, the horizontal principal stress is calculated using the Huang model.

7. The seismic assessment method of fault stability according to claim 1, characterized in that: The fracture stability prediction plane diagram is obtained based on the fracture surface shear stress and effective normal stress, including: Compare the magnitudes of fracture surface shear stress and effective normal stress to obtain comparative results; A fracture stability prediction plane diagram is drawn according to the comparison result.

8. A seismic assessment device for fault stability, characterized in that: The device comprises: A basic data acquisition module is used to acquire basic seismic data of the drilling section; wherein the basic seismic data of the drilling section includes seismic data and actual drilling data; A fracture interpretation module, used to interpret fractures based on the actual drilling data and identify fracture development locations in the drilling section; A seismic attribute extraction module, used for extracting seismic attributes capable of predicting fractures based on the seismic data; A calibration module, used to calibrate the fracture development position of the drilling section and the seismic attribute, determine a sensitive attribute body, and determine a fracture quantitative characterization result based on the sensitive attribute body; An inversion module, used for performing elastic parameter inversion based on the seismic data to obtain an elastic parameter body, and determining a rock mechanics parameter body according to the elastic parameter body; A geostress prediction module, used for determining a geostress prediction result based on the elastic parameter body and the rock mechanics parameter body; The module for predicting the direction of geostress is used to establish geological, mechanical and mathematical models based on the thin plate theory and the rock mechanics parameter body, and to determine the prediction results of the direction of geostress using the three-dimensional finite element difference simulation method; The fracture stability prediction module is used to calculate the fracture surface shear stress and effective normal stress based on the Mohr-Coulomb theory and the fracture quantitative characterization results, combined with the ground stress prediction results and the ground stress direction prediction results, and obtain the fracture stability prediction plane diagram based on the fracture surface shear stress and the effective normal stress.

9. A processor, characterized in that: A method for seismic assessment of fault stability configured to perform any one of claims 1 to 7.

10. A machine-readable storage medium having instructions stored thereon, characterized in that: When the instructions are executed by a processor, the processor is configured to perform the seismic assessment method of fault stability according to any one of claims 1 to 7.

Images