CN-121978754-A - Fracture basin sand body distribution prediction method, device, equipment, storage medium and product based on fault activity quantitative characterization

CN121978754ACN 121978754 ACN121978754 ACN 121978754ACN-121978754-A

Abstract

The application provides a fracture basin sand body distribution prediction method, device, equipment, storage medium and product based on fault activity quantitative characterization, and relates to the technical field of sand body prediction, comprising the steps of acquiring three-dimensional seismic data and logging data of a target area, and performing construction interpretation on the three-dimensional seismic data to obtain construction interpretation data; the method comprises the steps of calculating a thickness difference grid between a fault upper disc and a fault lower disc of a target fault in a target interval based on construction interpretation data, calculating a fault activity index of the target fault in a target interval deposition period based on the construction interpretation data, determining a tolerant space and supply ratio of the target interval based on the construction interpretation data and logging data, and generating a sand thickness distribution map of the target interval based on the thickness difference grid, the fault activity index and the tolerant space and supply ratio. According to the application, more accurate sand body distribution prediction results can be obtained.

Inventors

LIU YICHEN
XU QINGYAN
LIU YUMING
BAO LEI
LIU PEIPEI
CHEN QI
LIU HAOCHEN
ZHANG FAN
MA XUAN
LIU HENGZHI

Assignees

中国石油大学(北京)

Dates

Publication Date: 20260505
Application Date: 20251223

Claims (10)

1. A fracture basin sand body distribution prediction method based on fault activity quantitative characterization is characterized by comprising the following steps: acquiring three-dimensional seismic data and logging data of a target area, and performing construction interpretation on the three-dimensional seismic data to obtain construction interpretation data; Calculating a thickness difference grid between an upper fault disc and a lower fault disc of the target fault within the target interval based on the construction interpretation data; Calculating a fault activity index of the target fault at a target interval deposition period based on the formation interpretation data; determining a ratio of available space to supply for the target interval based on the formation interpretation data and the logging data; And generating a sand thickness profile of the target interval based on the thickness difference grid, the fault activity index, and the space-to-supply ratio.
2. The method for predicting the distribution of fractured basin sand based on quantitative characterization of fault activity according to claim 1, wherein before the step of obtaining three-dimensional seismic data and logging data of the target area, further comprises: acquiring initial logging data of the target area; Performing depth correction processing, same-layer alignment processing and curve standard processing on the initial logging data to obtain standardized logging data; lithology calibration is carried out on the standardized logging data by using core data and sheet analysis data to obtain calibrated logging data; and carrying out standardization processing on the calibrated logging data to obtain logging data.
3. The fault basin sand distribution prediction method based on quantitative characterization of fault activity of claim 1, wherein the calculating a thickness difference grid between a fault upper plate and a fault lower plate of a target fault within a target interval based on the formation interpretation data comprises: Extracting discrete deposition thickness data on two sides of a fault in the target interval from the construction interpretation data; Calculating a thickness difference value of discrete distribution between the fault upper disc and the fault lower disc based on the discrete deposition thickness data; taking the trend of the target fault as a main change path direction, and adopting an anisotropic kriging interpolation algorithm to perform spatial interpolation on the thickness difference values of the discrete distribution to generate a continuous thickness difference initial grid; And applying structural consistency space smoothness constraint to the continuous thickness difference initial grid to generate a thickness difference grid.
4. The fault basin sand distribution prediction method based on quantitative characterization of fault activity of claim 1, wherein the calculating a fault activity index of the target fault during a target interval deposition period based on the formation interpretation data comprises: Based on the construction interpretation data, acquiring the vertical sedimentation rate of the target fault in the target interval sedimentation period through a construction inversion technology and a stratum recovery technology, and acquiring the stratum thickness of the target interval; identifying different cutting layers of the target fault from the construction interpretation data, and determining the number of activities of the target fault according to the number of the cutting layers; Extracting late displacement amplitude of the target fault in late construction activities from the construction interpretation data; a fault activity index is calculated based on the vertical settlement rate, the formation thickness, the number of activities, and the late displacement amplitude.
5. The fault basin sand distribution prediction method based on quantitative characterization of fault activity of claim 1, wherein the generating a sand thickness profile for the target interval based on the thickness difference grid, the fault activity index, and the space-to-supply ratio comprises: constructing a leachable mapping relation with the thickness of the sand body as a baseline prediction model by taking the thickness difference grid, the fault activity index and the ratio of the allowable space to the supply as input variables; training the baseline prediction model by adopting a weighted least square algorithm to obtain a parameter estimation value of the baseline prediction model; And generating an initial prediction result based on the parameter estimation value, applying Laplace regularization constraint to the initial prediction result, and generating a sand thickness distribution map by minimizing the spatial second-order differential of the prediction result.
6. The method for predicting a fault-activity quantitative characterization-based fractured-basin sand body according to claim 5, wherein the steps of generating an initial prediction result based on the parameter estimation value, applying a laplace regularization constraint to the initial prediction result, and generating a sand body thickness distribution map by minimizing a spatial second order differential of the prediction result, further comprise: Based on a Monte Carlo and a self-help sampling method, obtaining a plurality of new input variables after performing a plurality of disturbance on the input variables; Generating a plurality of groups of sand thickness distribution graphs based on the plurality of new input variables; Calculating a sand thickness confidence interval on each spatial position based on the multiple groups of sand thickness distribution diagrams, and generating a confidence level plan corresponding to each group of sand thickness distribution diagrams; and generating a comprehensive prediction report based on the multiple groups of sand thickness distribution diagrams and the corresponding confidence level plan.
7. The utility model provides a fracture basin sand body distribution prediction unit based on fault activity ration characterization which characterized in that includes: The acquisition module is used for acquiring three-dimensional seismic data and logging data of the target area, and performing construction interpretation on the three-dimensional seismic data to obtain construction interpretation data; A calculation module for calculating a thickness difference grid between an upper fault disc and a lower fault disc of a target fault within a target interval based on the formation interpretation data; the calculation module is used for calculating a fault activity index of the target fault in a target interval deposition period based on the construction interpretation data; A determination module for determining a ratio of the volume of the target interval to the supply based on the formation interpretation data and the logging data; and the generation module is used for generating a sand thickness distribution map of the target interval based on the thickness difference grid, the fault activity index and the allowable space and supply ratio.
8. An electronic device comprising a memory, a processor and a computer program stored on the memory and running on the processor, wherein the processor when executing the computer program implements the fault-activity quantitative characterization-based fracture basin sand body distribution prediction method of any one of claims 1 to 6.
9. A non-transitory computer readable storage medium having stored thereon a computer program, which when executed by a processor, implements the fracture basin sand body distribution prediction method based on quantitative characterization of fault activity according to any one of claims 1 to 6.
10. A computer program product comprising a computer program, characterized in that the computer program, when executed by a processor, implements a fracture basin sand body distribution prediction method based on quantitative characterization of fault activity according to any one of claims 1 to 6.

Description

Fracture basin sand body distribution prediction method, device, equipment, storage medium and product based on fault activity quantitative characterization Technical Field The invention relates to the technical field of sand body prediction, in particular to a fracture basin sand body distribution prediction method, device, equipment, storage medium and product based on fault activity quantitative characterization. Background Fault activity is a key factor in controlling deposit basin architecture evolution and deposit system development. The formation, evolution and differential activities of faults not only control the spatial spread and thickness change of stratum, but also directly influence the migration path, the material source supply direction and the geometric form and spatial connectivity of sand body of a deposition system. Especially in land-phase fracture basin (such as Bohai Bay basin, songlao basin and other newly-generated fracture areas), faults develop densely and have strong activities, and the differential activities often cause migration of a deposition center and periodical change of a receivable space, so that a multi-level and overlapped deposition system is formed, the sand body is discontinuous in the vertical direction and asymmetric in the plane, and the heterogeneity of a reservoir is obviously enhanced. In recent years, with the development of high-resolution three-dimensional earthquake, logging interpretation and geologic modeling technologies, the accuracy of fault identification and sedimentary system analysis is remarkably improved. The existing method mainly discloses the correlation between fault activity and sand body spreading qualitatively or semi-quantitatively through means of construction explanation, seismic attribute analysis, logging phase identification and the like. However, these methods are mostly focused on static descriptions or empirical statistics, and fail to establish a systematic and quantitative response relationship between fault activity parameters and sand distribution characteristics, which results in difficulty in meeting the requirements of fine exploration in areas with low well control or complicated geological conditions. Disclosure of Invention The invention provides a fracture basin sand body distribution prediction method, device, equipment, storage medium and product based on fault activity quantitative characterization, which are used for solving the defect of low accuracy of sand body distribution prediction in the prior art and realizing high-accuracy sand body distribution prediction. The invention provides a fracture basin sand body distribution prediction method based on fault activity quantitative characterization, which comprises the following steps: acquiring three-dimensional seismic data and logging data of a target area, and performing construction interpretation on the three-dimensional seismic data to obtain construction interpretation data; Calculating a thickness difference grid between an upper fault disc and a lower fault disc of the target fault within the target interval based on the construction interpretation data; Calculating a fault activity index of the target fault at a target interval deposition period based on the formation interpretation data; determining a ratio of available space to supply for the target interval based on the formation interpretation data and the logging data; And generating a sand thickness profile of the target interval based on the thickness difference grid, the fault activity index, and the space-to-supply ratio. According to the fracture basin sand body distribution prediction method based on fault activity quantitative characterization provided by the invention, before the step of obtaining the three-dimensional seismic data and the logging data of the target area, the method further comprises the following steps: acquiring initial logging data of the target area; Performing depth correction processing, same-layer alignment processing and curve standard processing on the initial logging data to obtain standardized logging data; lithology calibration is carried out on the standardized logging data by using core data and sheet analysis data to obtain calibrated logging data; and carrying out standardization processing on the calibrated logging data to obtain logging data. According to the fracture basin sand body distribution prediction method based on quantitative fault activity characterization, the thickness difference grid between the fault upper disc and the fault lower disc of the target fault is calculated in the target interval based on the construction interpretation data, and the method comprises the following steps: Extracting discrete deposition thickness data on two sides of a fault in the target interval from the construction interpretation data; Calculating a thickness difference value of discrete distribution between the fault upper disc and the fault lower disc based o