阅读说明：本技术 基于地质-地球物理的沉积相演化定量表征方法及系统 (Sedimentary facies evolution quantitative characterization method and system based on geology-geophysical ) 是由 张旺 梁裳恣 单小彩 于 2021-09-03 设计创作，主要内容包括：本发明属于沉积相定量表征领域,具体涉及了一种基于地质-地球物理的沉积相演化定量表征方法及系统,旨在解决在地质模型约束下,现有技术无法系统化、体系化地利用地球物理技术实现对沉积相演化的定量表征的问题。本发明包括：利用岩心-测井-地震,建立高分辨率层序地层格架,并获取地震数据体上的三维层序界面；结合表征沉积期地势高低的旅行时差,获取沉积期近似古地貌图；以高分辨率层序格架为模型,开展测井-地震联合波阻抗反演；波阻抗体地层切片属性分析；建立砂地比-波阻抗映射关系；基于砂地比-波阻抗映射关系,动态、定量表征沉积相演化过程。本发明可以实现高分辨率层序格架内沉积相演化过程的动态、定量表征。(The invention belongs to the field of sedimentary facies quantitative characterization, and particularly relates to a sedimentary facies evolution quantitative characterization method and system based on geology-geophysical, aiming at solving the problem that the prior art cannot systematically and systematically utilize the geophysical technology to realize quantitative characterization of sedimentary facies evolution under the constraint of a geological model. The invention comprises the following steps: establishing a high-resolution sequence stratigraphic framework by using rock core-logging-earthquake, and acquiring a three-dimensional sequence interface on an earthquake data body; acquiring an approximate ancient apparent map of the sedimentary period by combining the travel time difference representing the terrain height of the sedimentary period; carrying out logging-seismic combined wave impedance inversion by taking the high-resolution sequence grid as a model; analyzing the attribute of the wave impedance body stratigraphic slice; establishing a sand-ground ratio-wave impedance mapping relation; and dynamically and quantitatively representing the sedimentary phase evolution process based on the sand-ground ratio-wave impedance mapping relation. The invention can realize dynamic and quantitative characterization of the evolution process of the sedimentary phase in the high-resolution sequence grid.)

1. A geologic-geophysical based depositional phase evolution quantitative characterization method is characterized by comprising the following steps:

step S10, based on the obtained core sample, well logging and seismic data volume, establishing a high-resolution stratigraphic framework on the seismic data volume by using a well logging-seismic calibration and combining seismic sequence boundary identification through a core and well logging sequence boundary identification technology;

step S20, acquiring a three-dimensional sequence interface on the seismic data body through core-logging-seismic comprehensive calibration based on the high-resolution stratigraphic framework on the seismic data body;

step S30, obtaining an approximate ancient apparent map of the depositional period based on a three-dimensional sequence interface on the seismic data volume and by combining with a travel time difference representing the terrain height of the depositional period;

s40, performing logging-seismic combined wave impedance inversion by taking the high-resolution stratum lattice frame as a model to obtain a wave impedance data volume;

step S50, acquiring root mean square wave impedance of the wave impedance body based on the high-resolution stratum lattice and the wave impedance data body, and acquiring sand-ground ratio of the sequence lattice on a single logging well; the sand-to-ground ratio is the ratio of the thickness of a sand layer of a target layer to the thickness of a stratum;

step S60, establishing a sand-ground ratio-wave impedance mapping relation based on the interlayer attributes of the wave impedance body and the sand-ground ratio of the sequence grid;

and step S70, realizing the dynamic evolution quantitative characterization of the sedimentary facies based on the sand-to-wave impedance mapping relation and the sedimentary period approximate ancient apparent map, the direction of an analyte source, the sand body distribution range and the accumulation advancing or accumulation retreating process.

2. The geologic-geophysical based depositional phase evolution quantitative characterization method of claim 1, wherein the high resolution stratigraphic framework, its high resolution sequence cyclic partitioning method comprises:

the stratum unit section is in the sequence level and corresponds to long-term convolution;

the stratum unit sand layer group is a quasi-stratum sequence group, corresponds to a middle period cycle and is a sedimentary facies characterization framework;

depositing a subphase development framework as a quasi-sequence corresponding to a short-term cycle.

3. The geologic-geophysical based depositional phase evolution quantitative characterization method according to claim 1, wherein said depositional phase is approximately ancient apparent map as a three-dimensional display stratigraphic map corresponding to a seismic two-way trip;

the calculation method of the earthquake two-way travel time comprises the following steps:

wherein the content of the first and second substances,when the earthquake is represented as a two-way trip,when the travel of the bottom interface of the layer sequence is represented,representing the time of travel of the top interface of the sequence.

4. The geologic-geophysical based depositional phase evolution quantitative characterization method of claim 1, wherein the root mean square wave impedance is calculated by:

wherein the content of the first and second substances,the root mean square wave impedance is represented,in-layer first representing wave impedanceThe value of the wave impedance at each sampling point,the number of sampling points in the wave resistor.

5. The geologic-geophysical based quantitative characterization method of dephase evolution according to claim 4, wherein the sand-to-ground ratio of the sequence trellis is calculated by:

wherein the content of the first and second substances,represents the sand-ground ratio of the sequence lattice,represents the thickness of the sand layer of the target layer,representing the thickness of the formation.

6. The geologic-geophysical based depositional phase evolution quantitative characterization method of claim 5, wherein the sand-to-ground ratio-wave impedance mapping relationship is:

wherein the content of the first and second substances,the mapping function is fitted to preset statistics.

7. The geologic-geophysical based depositional facies evolution quantitative characterization method of claim 6, wherein the sand-to-ground ratios of the sequence trellis are used to characterize the geomorphic depositional facies;

if the sand-ground ratio of the sequence grid is higher than a set first threshold value, the sand content of the stratum is high, the hydrodynamic force in the sedimentation period is strong, and the stratum is deposited near a source;

if the sand-ground ratio of the sequence grid is lower than a set second threshold value, the sand content of the stratum is low, the hydrodynamic force in the sedimentation period is weak, and the stratum is sedimentated in a basin of a remote source;

wherein the first threshold is greater than the second threshold.

8. A geologic-geophysical based depositional phase evolution quantitative characterization system, the system comprising:

the high-resolution stratigraphic framework establishing module is configured to establish a high-resolution stratigraphic framework on the seismic data volume by utilizing logging-seismic calibration and combining seismic sequence boundary identification through a rock core and logging sequence boundary identification technology based on the obtained rock core sample, the logging and the seismic data volume;

the three-dimensional sequence interface acquisition module is configured to acquire a three-dimensional sequence interface on the seismic data body through core-logging-seismic comprehensive calibration based on the high-resolution stratigraphic framework on the seismic data body;

the sedimentary period approximate ancient apparent map acquisition module is configured to acquire a sedimentary period approximate ancient apparent map based on a three-dimensional sequence of layers interface on the seismic data body in combination with a travel time difference representing the height of a sedimentary period topography;

the wave impedance data volume acquisition module is configured to perform logging-seismic combined wave impedance inversion by taking the high-resolution stratum lattice frame as a model to obtain a wave impedance data volume;

a sand-to-ground ratio acquisition module configured to acquire root mean square wave impedance of the wave impedance data volume based on the high resolution stratigraphic framework and the wave impedance data volume, and acquire sand-to-ground ratio of the sequence framework on a single log; the sand-to-ground ratio is the ratio of the thickness of a sand layer of a target layer to the thickness of a stratum;

a sand-to-ground ratio-wave impedance mapping establishing module configured to establish a sand-to-ground ratio-wave impedance mapping relation based on the interlayer properties of the wave impedance and the sand-to-ground ratio of the sequence trellis;

and the sedimentary phase dynamic evolution quantitative characterization module is configured to realize sedimentary phase dynamic evolution quantitative characterization based on the sand-to-wave impedance mapping relation and the sedimentary phase approximate ancient apparent map, the direction of an analyte source, the sand body distribution range and the accumulation advancing or accumulation retreating process.

9. An electronic device, comprising:

at least one processor; and

a memory communicatively coupled to at least one of the processors; wherein the content of the first and second substances,

the memory stores instructions executable by the processor for execution by the processor to implement the geologic-geophysical based depositional phase evolution quantitative characterization method of any of claims 1-7.

10. A computer-readable storage medium storing computer instructions for execution by the computer to implement the geologic-geophysical based depositional phase evolution quantitative characterization method of any one of claims 1-7.

Technical Field

The invention belongs to the field of sedimentary facies quantitative characterization, and particularly relates to a sedimentary facies evolution quantitative characterization method and system based on geology-geophysical characteristics.

Background

Sedimentary facies research is one of important links in selected area belt fixing and trap forecasting research in oil and gas exploration. Sedimentary facies beneficial to oil and gas accumulation, such as conventional oil and gas accumulation delta front edge branch river sand bodies, biological reefs, unconventional oil and gas accumulation deep sea-semi-deep sea terracotta facies or high-quality shale of deep lakes, are important exploration targets. Therefore, how to quantitatively predict and characterize depositional facies has been a key technical problem in the field of oil and gas exploration.

The main sedimentary facies characterization techniques in the prior art mainly include two types of geological techniques and geophysical techniques: the geological technology comprises outcrop, core and slice observation, the determination of sedimentary facies is more intuitive and strong, but most of the sedimentary facies are descriptive or semi-quantitative characterization and are difficult to predict by lateral extrapolation; the geophysical technology comprises well logging, earthquake and the like, can quantitatively represent and transversely predict geological abnormal bodies, but uncertainty exists in geological meaning analysis.

In general, the comprehensive geologic-geophysical characterization of dephasis is currently one of the most heavily developed exploration techniques in the field of oil and gas exploration. How to realize quantitative characterization of sedimentary facies evolution by using a geophysical technology under the constraint of a geological model is an exploration practice problem which needs to be solved urgently at present, and a systematized solution aiming at the problem is not developed in the field.

Disclosure of Invention

In order to solve the above problems in the prior art, that is, under the constraint of a geological model, the prior art cannot systematically and systematically utilize a geophysical technology to realize quantitative characterization of dephasic evolution, the invention provides a geological-geophysical-based method for quantitatively characterizing dephasic evolution, which comprises the following steps:

step S10, based on the obtained core sample, well logging and seismic data volume, establishing a high-resolution stratigraphic framework on the seismic data volume by using a well logging-seismic calibration and combining seismic sequence boundary identification through a core and well logging sequence boundary identification technology;

step S20, acquiring a three-dimensional sequence interface on the seismic data body through core-logging-seismic comprehensive calibration based on the high-resolution stratigraphic framework on the seismic data body;

step S30, obtaining an approximate ancient apparent map of the depositional period based on a three-dimensional sequence interface on the seismic data volume and by combining with a travel time difference representing the terrain height of the depositional period;

s40, performing logging-seismic combined wave impedance inversion by taking the high-resolution stratum lattice frame as a model to obtain a wave impedance data volume;

step S50, acquiring root mean square wave impedance of the wave impedance body based on the high-resolution stratum lattice and the wave impedance data body, and acquiring sand-ground ratio of the sequence lattice on a single logging well; the sand-to-ground ratio is the ratio of the thickness of a sand layer of a target layer to the thickness of a stratum;

step S60, establishing a sand-ground ratio-wave impedance mapping relation based on the interlayer attributes of the wave impedance body and the sand-ground ratio of the sequence grid;

and step S70, realizing the dynamic evolution quantitative characterization of the sedimentary facies based on the sand-to-wave impedance mapping relation and the sedimentary period approximate ancient apparent map, the direction of an analyte source, the sand body distribution range and the accumulation advancing or accumulation retreating process.

In some preferred embodiments, the high resolution stratigraphic framework, the high resolution sequence cyclic division method thereof comprises:

the stratum unit section is in the sequence level and corresponds to long-term convolution;

the stratum unit sand layer group is a quasi-stratum sequence group, corresponds to a middle period cycle and is a sedimentary facies characterization framework;

depositing a subphase development framework as a quasi-sequence corresponding to a short-term cycle.

In some preferred embodiments, the dephasing profile is approximately ancient in a three-dimensional display horizon corresponding to a two-way seismic trip;

the calculation method of the earthquake two-way travel time comprises the following steps:

wherein the content of the first and second substances,when the earthquake is represented as a two-way trip,when the travel of the bottom interface of the layer sequence is represented,representing the time of travel of the top interface of the sequence.

In some preferred embodiments, the root mean square wave impedance is calculated by:

wherein the content of the first and second substances,the root mean square wave impedance is represented,in-layer first representing wave impedanceThe value of the wave impedance at each sampling point,the number of sampling points in the wave resistor.

In some preferred embodiments, the sand-ground ratio of the sequence grid is calculated by:

wherein the content of the first and second substances,represents the sand-ground ratio of the sequence lattice,represents the thickness of the sand layer of the target layer,representing the thickness of the formation.

In some preferred embodiments, the sand-to-ground-wave impedance mapping relationship is:

wherein the content of the first and second substances,the mapping function is fitted to preset statistics.

In some preferred embodiments, the sand-to-ground ratio of the sequence lattice is used to characterize the relief dephasing;

if the sand-ground ratio of the sequence grid is higher than a set first threshold value, the sand content of the stratum is high, the hydrodynamic force in the sedimentation period is strong, and the stratum is deposited near a source;

if the sand-ground ratio of the sequence grid is lower than a set second threshold value, the sand content of the stratum is low, the hydrodynamic force in the sedimentation period is weak, and the stratum is sedimentated in a basin of a remote source;

wherein the first threshold is greater than the second threshold.

In another aspect of the present invention, a geologic-geophysical based depositional phase evolution quantitative characterization system is provided, which includes:

the high-resolution stratigraphic framework establishing module is configured to establish a high-resolution stratigraphic framework on the seismic data volume by utilizing logging-seismic calibration and combining seismic sequence boundary identification through a rock core and logging sequence boundary identification technology based on the obtained rock core sample, the logging and the seismic data volume;

the three-dimensional sequence interface acquisition module is configured to acquire a three-dimensional sequence interface on the seismic data body through core-logging-seismic comprehensive calibration based on the high-resolution stratigraphic framework on the seismic data body;

the sedimentary period approximate ancient apparent map acquisition module is configured to acquire a sedimentary period approximate ancient apparent map based on a three-dimensional sequence of layers interface on the seismic data body in combination with a travel time difference representing the height of a sedimentary period topography;

the wave impedance data volume acquisition module is configured to perform logging-seismic combined wave impedance inversion by taking the high-resolution stratum lattice frame as a model to obtain a wave impedance data volume;

a sand-to-ground ratio acquisition module configured to acquire root mean square wave impedance of the wave impedance data volume based on the high resolution stratigraphic framework and the wave impedance data volume, and acquire sand-to-ground ratio of the sequence framework on a single log; the sand-to-ground ratio is the ratio of the thickness of a sand layer of a target layer to the thickness of a stratum;

a sand-to-ground ratio-wave impedance mapping establishing module configured to establish a sand-to-ground ratio-wave impedance mapping relation based on the interlayer properties of the wave impedance and the sand-to-ground ratio of the sequence trellis;

and the sedimentary phase dynamic evolution quantitative characterization module is configured to realize sedimentary phase dynamic evolution quantitative characterization based on the sand-to-wave impedance mapping relation and the sedimentary phase approximate ancient apparent map, the direction of an analyte source, the sand body distribution range and the accumulation advancing or accumulation retreating process.

In a third aspect of the present invention, an electronic device is provided, including:

at least one processor; and

a memory communicatively coupled to at least one of the processors; wherein the content of the first and second substances,

the memory stores instructions executable by the processor for execution by the processor to implement the geologic-geophysical based depositional phase evolution quantitative characterization method described above.

In a fourth aspect of the present invention, a computer-readable storage medium is provided, which stores computer instructions for being executed by the computer to implement the geologic-geophysical based depositional phase evolution quantitative characterization method.

The invention has the beneficial effects that:

(1) the sedimentary facies evolution quantitative characterization method based on geology-geophysics endows wave impedance attributes with definite geological significance through comprehensive analysis of geology and geophysics, and can quantitatively characterize sedimentary facies.

(2) The sedimentary facies evolution quantitative characterization method based on geology-geophysics can realize dynamic characterization of the sedimentary facies evolution process in the high-resolution sequence grid.

(3) The sedimentary facies evolution quantitative characterization method based on geology-geophysics can realize dynamic and quantitative characterization of sedimentary facies and provide technical support for oil-gas exploration.

Drawings

Other features, objects and advantages of the present application will become more apparent upon reading of the following detailed description of non-limiting embodiments thereof, made with reference to the accompanying drawings in which:

FIG. 1 is a schematic flow chart of a geologic-geophysical based depositional phase evolution quantitative characterization method of the present invention;

FIG. 2 is an ancient apparent map of a dephasing period of an embodiment of the geologic-geophysical based quantitative characterization method for dephasing evolution of the invention;

FIG. 3 is a root mean square wave impedance plane diagram of an embodiment of the geologic-geophysical based depositional phase evolution quantitative characterization method of the present invention;

FIG. 4 is a sand-to-earth ratio-root-mean-square relationship diagram of an embodiment of the geologic-geophysical based depositional facies evolution quantitative characterization method of the present invention;

FIG. 5 is a sedimentary facies evolution quantitative characterization and evolution analysis diagram of an embodiment of the sedimentary facies evolution quantitative characterization method based on geology-geophysics.

Detailed Description

The present application will be described in further detail with reference to the following drawings and examples. It is to be understood that the specific embodiments described herein are merely illustrative of the relevant invention and not restrictive of the invention. It should be noted that, for convenience of description, only the portions related to the related invention are shown in the drawings.

It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict. The present application will be described in detail below with reference to the embodiments with reference to the attached drawings.

The invention relates to a geologic-geophysical based sedimentary facies evolution quantitative characterization method, which comprises the following steps:

step S10, based on the obtained core sample, well logging and seismic data volume, establishing a high-resolution stratigraphic framework on the seismic data volume by using a well logging-seismic calibration and combining seismic sequence boundary identification through a core and well logging sequence boundary identification technology;

step S20, acquiring a three-dimensional sequence interface on the seismic data body through core-logging-seismic comprehensive calibration based on the high-resolution stratigraphic framework on the seismic data body;

step S30, obtaining an approximate ancient apparent map of the depositional period based on a three-dimensional sequence interface on the seismic data volume and by combining with a travel time difference representing the terrain height of the depositional period;

s40, performing logging-seismic combined wave impedance inversion by taking the high-resolution stratum lattice frame as a model to obtain a wave impedance data volume;

step S50, acquiring root mean square wave impedance of the wave impedance body based on the high-resolution stratum lattice and the wave impedance data body, and acquiring sand-ground ratio of the sequence lattice on a single logging well; the sand-to-ground ratio is the ratio of the thickness of a sand layer of a target layer to the thickness of a stratum;

step S60, establishing a sand-ground ratio-wave impedance mapping relation based on the interlayer attributes of the wave impedance body and the sand-ground ratio of the sequence grid;

and step S70, realizing the dynamic evolution quantitative characterization of the sedimentary facies based on the sand-to-wave impedance mapping relation and the sedimentary period approximate ancient apparent map, the direction of an analyte source, the sand body distribution range and the accumulation advancing or accumulation retreating process.

In order to more clearly illustrate the geologic-geophysical-based depositional facies evolution quantitative characterization method of the present invention, the following describes in detail the steps in the embodiment of the present invention with reference to fig. 1.

The geologic-geophysical-based depositional facies evolution quantitative characterization method of the first embodiment of the invention comprises the steps of S10-S70, wherein the steps are described in detail as follows:

and step S10, based on the obtained core sample, well logging and seismic data volume, establishing a high-resolution stratum framework on the seismic data volume by using a well logging-seismic calibration and combining seismic sequence boundary identification through a core and well logging sequence boundary identification technology.

And step S20, acquiring a three-dimensional sequence interface on the seismic data body through core-logging-seismic comprehensive calibration based on the high-resolution stratigraphic framework on the seismic data body.

The high-resolution stratum lattice has high-resolution sequence cyclic division method including:

the stratum unit section is in the sequence level and corresponds to long-term convolution;

the stratum unit sand layer group is a quasi-stratum sequence group, corresponds to a middle period cycle and is a sedimentary facies characterization framework;

depositing a subphase development framework as a quasi-sequence, and correspondingly revolving in a short period;

sedimentary microfacies are generally not easily identified from seismic sequence grids due to limited seismic resolution.

And step S30, obtaining an approximate ancient apparent map of the depositional period based on the three-dimensional sequence interface on the seismic data volume and by combining the travel time difference representing the terrain height of the depositional period.

The sedimentary phase is similar to an ancient apparent map and is a three-dimensional display layer map corresponding to earthquake double-travel, and the calculation method of the earthquake double-travel is shown as the formula (1):

wherein the content of the first and second substances,when the earthquake is represented as a two-way trip,representing the bottom interface of the layer sequenceWhen the user travels, the user can select the time,representing the time of travel of the top interface of the sequence.

Approximating an ancient map, it is possible to indicate the height of the terrain during the deposition phase, where the time difference is small and the terrain is high; the time difference is large and the terrain is low. The accurate two-way travel time difference graph can reflect landform phenomena such as ancient river valleys, ancient sedimentary slope folds and the like.

As shown in FIG. 2, the depositional phase ancient apparent map of an embodiment of the geologic-geophysical based depositional phase evolution quantitative characterization method of the invention is shown, wherein light colors represent highland, and dark colors represent lowland.

And step S40, performing logging-seismic combined wave impedance inversion by taking the high-resolution stratum lattice as a model to obtain a wave impedance data volume.

After the step is carried out on the basis of rock physical modeling analysis, wave impedance is determined and can be used for distinguishing lithology of sandstone and mudstone. The wave impedance data volume obtained by the steps has definite lithology indication meaning and higher resolution.

The root mean square wave impedance is calculated according to the following formula (2):

wherein the content of the first and second substances,the root mean square wave impedance is represented,in-layer first representing wave impedanceThe value of the wave impedance at each sampling point,the number of sampling points in the wave resistor.

As shown in fig. 3, which is a root mean square wave impedance plane diagram of an embodiment of the geologic-geophysical based depositional phase evolution quantitative characterization method of the present invention, a light color is a sector range.

Step S50, acquiring root mean square wave impedance of the wave impedance body based on the high-resolution stratum lattice and the wave impedance data body, and acquiring sand-ground ratio of the sequence lattice on a single logging well; the sand ground ratio is the ratio of the sand layer thickness of a target layer to the formation thickness, and has sedimentary facies indication significance, such as high sand ground ratio, indication of high sand content of the formation, strong hydrodynamic force during sedimentary period and deposition of a near source, and low sand ground ratio, indication of low sand content of the formation, weak hydrodynamic force during sedimentary period and deposition of a basin of a far source.

The sand-ground ratio of the sequence grillwork is calculated according to the formula (3):

wherein the content of the first and second substances,represents the sand-ground ratio of the sequence lattice,represents the thickness of the sand layer of the target layer,representing the thickness of the formation.

The sand-ground ratio of the sequence lattice is used for representing landform sedimentary facies;

if the sand-ground ratio of the sequence grid is higher than a set first threshold value, the sand content of the stratum is high, the hydrodynamic force in the sedimentation period is strong, and the stratum is deposited near a source;

if the sand-ground ratio of the sequence grid is lower than a set second threshold value, the sand content of the stratum is low, the hydrodynamic force in the sedimentation period is weak, and the stratum is sedimentated in a basin of a remote source;

wherein the first threshold is greater than the second threshold.

As shown in fig. 4, a sand-ground ratio-root mean square relationship diagram of an embodiment of the geologic-geophysical-based depositional phase evolution quantitative characterization method of the present invention is shown, the left table is numerical statistics of the sand-ground ratio and the root mean square wave impedance, and the right table is a linear fitting and function of the sand-ground ratio and the root mean square wave impedance.

Step S60, establishing a sand-ground ratio-wave impedance mapping relation based on the interlayer attributes of the wave impedance body and the sand-ground ratio of the sequence grid, as shown in formula (4):

wherein the content of the first and second substances,the mapping function is fitted to the preset statistics, typically a positive correlation function.

In one embodiment of the present invention, the first and second substrates are,

and step S70, realizing the dynamic evolution quantitative characterization of the sedimentary facies based on the sand-to-wave impedance mapping relation and the sedimentary period approximate ancient apparent map, the direction of an analyte source, the sand body distribution range and the accumulation advancing or accumulation retreating process.

The deposition period of the step S30 is similar to an ancient apparent map, and the direction of the object source can be judged, and the deposit is converged from a raised position to a low position; the plane distribution of the interlayer attributes (many interlayer attributes, such as root mean square wave impedance, maximum wave impedance, minimum wave impedance and the like, and root mean square wave impedance obtained by calculating root mean square value through the wave impedance value of the interlayer sampling points) in the step S50, and the displayed fan shape can also judge the source direction and is unfolded from the source direction by the fan; the dynamic representation of the sand body distribution comprises plane position and area change. The accumulation advancing or accumulation retreating process can be judged according to the position of the fan body, and the movement of the fan body towards the object source is the accumulation retreating process, otherwise, the movement is the accumulation advancing process.

As shown in fig. 5, a sedimentary facies evolution quantitative characterization and evolution analysis chart of an embodiment of the sedimentary facies evolution quantitative characterization method based on geology-geophysical of the present invention can quantitatively characterize the fan body evolution sector delta in the depopulation process in the northwest direction, and the sand-ground ratio and the root mean square wave impedance in this embodiment have good linesRelationship of (A) to (B)R ²= 0.7313), as shown in formula (5):

wherein the content of the first and second substances,represents the sand-ground ratio of the sequence lattice,representing the rms wave impedance.

The range with large root mean square wave impedance is an underwater diversion river channel development area at the front edge of the fan delta, and the range with small root mean square wave impedance is a estuary dam and a mat sand development area at the front edge of the fan delta. The dynamic evolution of the dephasing can also be revealed: from the first sector to the third sector deposition period, the sand development area in the well area of gram 82 is favorable to moving towards the source. In the first stage and the second stage of fan body deposition period, the gram 82 well area mainly develops an underwater diversion river channel, and in the third stage of fan body deposition period, the gram 82 well area is favorable for sand body to not develop.

Although the foregoing embodiments describe the steps in the above sequential order, those skilled in the art will understand that, in order to achieve the effect of the present embodiments, the steps may not be executed in such an order, and may be executed simultaneously (in parallel) or in an inverse order, and these simple variations are within the scope of the present invention.

The geologic-geophysical based depositional phase evolution quantitative characterization system of the second embodiment of the present invention comprises:

the high-resolution stratigraphic framework establishing module is configured to establish a high-resolution stratigraphic framework by utilizing a core and logging sequence boundary identification technology, utilizing logging-seismic calibration and combining seismic sequence boundary identification;

the three-dimensional sequence interface acquisition module is configured to acquire a three-dimensional sequence interface on a seismic data body through core-logging-seismic comprehensive calibration;

the sedimentary period approximate ancient apparent map acquisition module is configured to acquire a sedimentary period approximate ancient apparent map based on a three-dimensional sequence of layers interface on the seismic data body in combination with a travel time difference representing the height of a sedimentary period topography;

the wave impedance data volume acquisition module is configured to perform logging-seismic combined wave impedance inversion by taking the high-resolution stratum lattice frame as a model to obtain a wave impedance data volume;

a sand-to-ground ratio acquisition module configured to acquire interlayer properties of the wave impedance data volume based on the high resolution stratigraphic framework and the wave impedance data volume and acquire a sand-to-ground ratio of the sequence framework on a single log; the sand-to-ground ratio is the ratio of the thickness of a sand layer of a target layer to the thickness of a stratum;

a sand-to-ground ratio-wave impedance mapping establishing module configured to establish a sand-to-ground ratio-wave impedance mapping relation based on the interlayer properties of the wave impedance and the sand-to-ground ratio of the sequence trellis;

and the sedimentary facies dynamic evolution quantitative characterization module is configured to realize the sedimentary facies dynamic evolution quantitative characterization based on the sand-to-ground ratio-wave impedance mapping relation, the analyte source direction, the sand body distribution range and the accumulation entering or accumulation retreating process.

It can be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working process and related description of the system described above may refer to the corresponding process in the foregoing method embodiments, and will not be described herein again.

It should be noted that the geologic-geophysical-based depositional facies evolution quantitative characterization system provided in the above embodiment is only illustrated by the division of the above functional modules, and in practical applications, the above functions may be allocated to different functional modules according to needs, that is, the modules or steps in the embodiment of the present invention are further decomposed or combined, for example, the modules in the above embodiment may be combined into one module, or may be further split into multiple sub-modules, so as to complete all or part of the above described functions. The names of the modules and steps involved in the embodiments of the present invention are only for distinguishing the modules or steps, and are not to be construed as unduly limiting the present invention.

An electronic apparatus according to a third embodiment of the present invention includes:

at least one processor; and

a memory communicatively coupled to at least one of the processors; wherein the content of the first and second substances,

the memory stores instructions executable by the processor for execution by the processor to implement the geologic-geophysical based depositional phase evolution quantitative characterization method described above.

A computer-readable storage medium of a fourth embodiment of the present invention stores computer instructions for execution by the computer to implement the geologic-geophysical based depositional phase evolution quantitative characterization method described above.

It can be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working processes and related descriptions of the storage device and the processing device described above may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.

Those of skill in the art would appreciate that the various illustrative modules, method steps, and modules described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that programs corresponding to the software modules, method steps may be located in Random Access Memory (RAM), memory, Read Only Memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art. To clearly illustrate this interchangeability of electronic hardware and software, various illustrative components and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as electronic hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.

The terms "first," "second," and the like are used for distinguishing between similar elements and not necessarily for describing or implying a particular order or sequence.

The terms "comprises," "comprising," or any other similar term are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus.

So far, the technical solutions of the present invention have been described in connection with the preferred embodiments shown in the drawings, but it is easily understood by those skilled in the art that the scope of the present invention is obviously not limited to these specific embodiments. Equivalent changes or substitutions of related technical features can be made by those skilled in the art without departing from the principle of the invention, and the technical scheme after the changes or substitutions can fall into the protection scope of the invention.