Method and system for identifying facies belt boundary of volcanic institution

Abstract

The invention provides a method and system for identifying a volcanic institution facies belt boundary, and belongs to the technical field of oil exploration and development. The method comprises the following steps: dividing a drilled facies belt according to three-dimensional seismic reflection characteristics and a preset type of lithofacies proportion of a drilled well; constructing a near-source facies belt pseudo well and a far-source facies belt pseudo well based on the division result of the drilled facies belt and the multiple types of drilled lithofacies; based on the near-source facies zone pseudo well and the far-source facies zone pseudo well, determining a value range corresponding to a preset type of lithofacies proportion of the volcanic institution; based on the three-dimensional seismic data volume of the volcanic institution, a planar graph representing a facies belt is obtained; according to the value range corresponding to the preset type lithofacies proportion of the volcanic institution, volcanic institution facies belt division is carried out on the planar graph representing the facies belt, and the boundary between the near-source facies belt and the far-source facies belt is obtained. Therefore, the boundary between the near-source phase belt and the far-source phase belt is accurately and effectively defined, the near-source favorable phase belt of the volcanic institution is quickly focused, and a foundation is laid for an exploration target.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of oil exploration and development, in particular to a method for identifying boundaries of a volcanic mechanism facies belt and a system for identifying boundaries of a volcanic mechanism facies belt. BACKGROUND

[0002] Due to the differences in lithology, rock fabric and genesis, volcanic facies are usually divided into five types, namely, explosive facies, overflow facies, extrusion facies, volcanic channel facies and volcanic sedimentary facies. It is difficult to track and identify these five types of facies on actual seismic data. In the practice of volcanic rock exploration and development, volcanic mechanisms are usually divided into three facies belts, namely, crater-near crater, near-source and far-source facies belts, according to the distance from the crater. The distribution of the gas reservoirs of the basic volcanic rocks in the Shishanling Formation of the Songnan fault depression is mainly controlled by the effective reservoirs. The primary porosity is more developed in the near crater-near source facies belt, the proportion of the high porosity and permeability zone is large, and the gas reservoirs are more enriched. Therefore, finding the boundaries between the near-source facies belt and the far-source facies belt is the key to finding the effective reservoirs. At present, the near-source facies belt is mainly depicted by analyzing the seismic reflection characteristics and judging the position of the crater through the geometric attributes such as the stratigraphic dip and the coherence, and then qualitatively judging the distribution of the near-source facies belt. There is no method for quantitatively predicting the boundaries between the near-source facies belt and the far-source facies belt of the volcanic mechanism.

[0003] The energy half-life attribute mainly reflects the change relationship of the seismic wave energy, which is generally used to indicate the change of the sedimentary environment and the lithology and facies, and has been applied in the identification of the dense sandstone such as turbidite. However, there is no data to prove that it can be applied in the identification of the volcanic rock facies and facies belt. Due to the high volcanic velocity, the top envelope of the volcanic mechanism has strong seismic reflection energy, and the volcanic eruption causes the explosive facies to have chaotic and weak reflection characteristics. Therefore, there is a certain difference in the energy half-life. In the process of volcanic eruption, the volcanic materials are accumulated near the crater, and the stratum thickness near the crater is large. Therefore, the combination of the two can be used as the main basis for judging the boundaries between the near-source facies belt and the far-source facies belt. The existing method can only find the position of the crater by extrapolating the seismic reflection characteristics and the seismic geometric attributes of the actual drilled well, and then qualitatively judge the distribution, which cannot effectively define the planar boundaries between the near-source facies belt and the far-source facies belt of the volcanic mechanism. SUMMARY

[0004] The purpose of the embodiments of the present application is to provide a method and system for identifying the boundaries of the facies belt of the volcanic mechanism, so as to at least solve the problem that the boundaries between the near-source facies belt and the far-source facies belt of the volcanic mechanism cannot be accurately and effectively defined.

[0005] In order to achieve the above-mentioned purpose, the first aspect of the present application provides a method for identifying the boundaries of the facies belt of the volcanic mechanism, comprising:

[0006] According to the three-dimensional seismic reflection characteristics and the preset type of rock facies proportion in a single volcanic eruption period of the drilled well, the drilled well facies belt is divided;

[0007] construct a near-source facies belt pseudo well and a far-source facies belt pseudo well based on the drilled well facies belt division result and the characteristic curves corresponding to the multiple types of lithofacies of the drilled well in a single volcanic eruption period;

[0008] determine a value range corresponding to a preset type of lithofacies proportion of the volcanic mechanism based on the characteristic curve of the near-source facies belt pseudo well and the characteristic curve of the far-source facies belt pseudo well;

[0009] obtain a plane map representing facies belts based on a three-dimensional seismic data volume of the volcanic mechanism;

[0010] perform facies belt division of the volcanic mechanism on the plane map representing facies belts according to the value range corresponding to the preset type of lithofacies proportion of the volcanic mechanism, to obtain a boundary between a near-source facies belt and a far-source facies belt in a single volcanic eruption period.

[0011] Optionally, the method for identifying the facies belt boundary of the volcanic mechanism further includes:

[0012] divide the volcanic eruption periods according to drilled well drilling data;

[0013] perform lithofacies division according to the lithofacies in each volcanic eruption period to obtain a lithofacies division result;

[0014] The lithofacies division result includes an explosive facies, an overflow facies, an extrusion facies, a volcanic channel facies, and a volcanic sedimentary facies.

[0015] Optionally, the characteristic curves include acoustic curves and density curves.

[0016] The near-source facies belt pseudo well and the far-source facies belt pseudo well are constructed based on the drilled well facies belt division result and the characteristic curves corresponding to the multiple types of lithofacies of the drilled well in a single volcanic eruption period, and include:

[0017] determine acoustic curves and density curves of the explosive facies, the overflow facies, and the volcanic sedimentary facies based on the acoustic curves and the density curves of the drilled well and the lithofacies division result;

[0018] construct the near-source facies belt pseudo well using the acoustic curves and the density curves of the explosive facies and the overflow facies of the drilled well and construct the far-source facies belt pseudo well using the acoustic curves and the density curves of the volcanic sedimentary facies of the drilled well based on the drilled well facies belt division result.

[0019] Optionally, the drilled well facies belt division according to the three-dimensional seismic reflection characteristics and the preset type of lithofacies proportion of the drilled well in a single volcanic eruption period includes:

[0020] perform drilled well facies belt division according to the three-dimensional seismic reflection characteristics and the volcanic sedimentary facies proportion of the drilled well in a single volcanic eruption period;

[0021] S = volcanic sedimentary facies ÷ sum of all lithofacies;

[0022] Wherein, S represents the proportion of volcanic sedimentary facies, S is less than 50% is a proximal facies belt, S is higher than 50% is a distal facies belt.

[0023] Optionally, the characteristic curve includes an acoustic curve and a density curve.

[0024] The characteristic curve based on the proximal facies belt pseudo well and the characteristic curve based on the distal facies belt pseudo well determine the value range corresponding to the preset type lithofacies proportion of the volcanic mechanism, including:

[0025] Based on the acoustic curve and the density curve of the proximal facies belt pseudo well and the acoustic curve and the density curve of the distal facies belt pseudo well, a two-dimensional seismic forward model is established by using a convolution model, ray tracing and wave equation method.

[0026] Based on the two-dimensional seismic forward model, an energy half-life time T curve and a thickness H curve corresponding to the volcanic eruption period are calculated.

[0027] Based on the energy half-life time T curve and the thickness H curve corresponding to the volcanic eruption period, an R curve is calculated by using the formula R = (T / H)*100, the value A corresponding to the volcanic sedimentary facies proportion in the R curve is analyzed, and the value range of the value A is obtained.

[0028] Optionally, according to the value range corresponding to the preset type lithofacies proportion of the volcanic mechanism, the volcanic mechanism facies belt is divided on the plane map representing the facies belt, and the boundary between the proximal facies belt and the distal facies belt corresponding to the volcanic eruption period is obtained, including:

[0029] Based on the value A corresponding to the volcanic sedimentary facies proportion in the R curve, the volcanic mechanism facies belt is divided on the plane map representing the facies belt, and the distal facies belt position and the proximal facies belt position of the volcanic mechanism are obtained.

[0030] Based on the distal facies belt position and the proximal facies belt position, the boundary between the proximal facies belt and the distal facies belt is obtained.

[0031] Optionally, the volcanic mechanism facies belt is divided on the plane map representing the facies belt based on the value A corresponding to the volcanic sedimentary facies proportion in the R curve, including:

[0032] The area greater than the value range of the value A on the plane map is taken as the distal facies belt of the volcanic mechanism, and the area less than the value range of the value A is taken as the proximal facies belt of the volcanic mechanism.

[0033] Optionally, the three-dimensional seismic data volume of the volcanic mechanism is obtained to obtain the plane map representing the facies belt, including:

[0034] based on the three-dimensional seismic data volume of the volcanic edifice, a planar diagram t of energy half-decay time of a single volcanic eruption episode, a bottom surface and a top surface of the volcanic edifice are obtained;

[0035] based on the bottom surface and the top surface of the single volcanic eruption episode, a planar diagram h of thickness of the single volcanic eruption episode is calculated;

[0036] a planar diagram of facies belts is obtained according to a formula r=(t / h)*100; wherein r is a planar diagram of a parameter corresponding to the facies belts.

[0037] The second aspect of the present application provides a system for identifying facies belt boundaries of a volcanic edifice, comprising:

[0038] a drilled well facies belt division module configured to divide drilled well facies belts according to three-dimensional seismic reflection characteristics and a preset type lithofacies proportion of drilled wells in a single volcanic eruption episode;

[0039] a pseudo-well construction module configured to construct a near-source facies belt pseudo-well and a far-source facies belt pseudo-well based on a drilled well facies belt division result and characteristic curves corresponding to multiple types of lithofacies of drilled wells in a single volcanic eruption episode;

[0040] a value range determination module configured to determine a value range corresponding to a preset type lithofacies proportion of the volcanic edifice based on a characteristic curve of the near-source facies belt pseudo-well and a characteristic curve of the far-source facies belt pseudo-well;

[0041] a planar diagram obtaining module configured to obtain a planar diagram of facies belts based on a three-dimensional seismic data volume of the volcanic edifice;

[0042] a volcanic edifice facies belt division module configured to divide volcanic edifice facies belts on the planar diagram of facies belts according to the value range corresponding to the preset type lithofacies proportion of the volcanic edifice, to obtain boundaries between near-source facies belts and far-source facies belts in a single volcanic eruption episode.

[0043] The third aspect of the present application provides a machine-readable storage medium, the machine-readable storage medium stores instructions, the instructions, when executed by a processor, cause the processor to be configured to execute the above-mentioned method for identifying facies belt boundaries of a volcanic edifice.

[0044] The fourth aspect of the present application provides an electronic device, the electronic device comprises a memory, a processor and a computer program stored in the memory and executable on the processor, and the processor executes the computer program to realize the above-mentioned method for identifying facies belt boundaries of a volcanic edifice.

[0045] According to the technical scheme, a method and system for identifying boundaries of volcanic mechanism facies belts are provided, and the method and system determine facies belt division of drilled wells in a single volcanic eruption period according to three-dimensional seismic reflection characteristics and proportions of volcanic sedimentary facies of the drilled wells in the single volcanic eruption period. Based on the facies belt division results of the drilled wells and characteristic curves corresponding to various types of lithofacies of the drilled wells in the single volcanic eruption period, a near-source facies belt pseudo-well and a far-source facies belt pseudo-well are constructed. According to the characteristic curve of the near-source facies belt pseudo-well and the characteristic curve of the far-source facies belt pseudo-well, a value range corresponding to a proportion of volcanic sedimentary facies of the volcanic mechanism is determined. According to the value range corresponding to the proportion of the volcanic sedimentary facies of the volcanic mechanism, facies belt division of the volcanic mechanism is performed on a planar graph representing facies belts obtained by interpreting three-dimensional seismic data of the volcanic mechanism, so as to obtain boundaries between near-source facies belts and far-source facies belts in the single volcanic eruption period. Thus, the method and system can quantitatively, accurately and effectively define the boundaries between the near-source facies belts and the far-source facies belts, quickly focus on near-source favorable facies belts of the volcanic mechanism, and are beneficial to further searching for effective volcanic reservoirs, and thus lay a foundation for exploration targets.

[0046] Other features and advantages of the present application will be explained in detail in the following detailed description. BRIEF DESCRIPTION OF DRAWINGS

[0047] The accompanying drawings are included to provide a further understanding of the present application and are incorporated in and constitute a part of this specification, illustrate embodiments of the present application and serve to explain the principles of the present application, but are not intended to limit the present application. In the drawings:

[0048] Figure 1 is a flowchart of a method for identifying boundaries of volcanic mechanism facies belts provided by an embodiment of the present application;

[0049] Figure 2 is a flowchart of another method for identifying boundaries of volcanic mechanism facies belts provided by an embodiment of the present application;

[0050] Figure 3 is a planar graph of energy half-life time of a volcanic mechanism provided by an embodiment of the present application;

[0051] Figure 4 is a thickness graph of a volcanic eruption period provided by an embodiment of the present application;

[0052] Figure 5 is a curve schematic diagram in a volcanic eruption period provided by an embodiment of the present application;

[0053] Figure 6 is a planar distribution graph of facies belt division of a volcanic mechanism provided by an embodiment of the present application;

[0054] Figure 7is a block diagram of a system for identifying a volcanic edifice facies belt boundary according to an embodiment of the present application;

[0055] Figure 8 is a schematic diagram of an electronic device structure according to a preferred embodiment of the present application.

[0056] Legend of reference signs

[0057] 10 - electronic device, 100 - processor, 101 - memory, 102 - computer program. DETAILED DESCRIPTION

[0058] The specific embodiments of the present application are described in detail below with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are merely intended to illustrate and explain the present application, and are not intended to limit the present application.

[0059] Figure 1 is a flowchart of a method for identifying a volcanic edifice facies belt boundary according to an embodiment of the present application, Figure 2 is a flowchart of another method for identifying a volcanic edifice facies belt boundary according to an embodiment of the present application. As shown in Figure 1 and Figure 2 The method for identifying a volcanic edifice facies belt boundary according to an embodiment of the present application comprises:

[0060] S110: According to the three-dimensional seismic reflection characteristics and the preset type lithofacies proportion of the drilled well in a single volcanic eruption period, the drilled well facies belt is divided.

[0061] The preset type lithofacies proportion is the volcanic sedimentary facies proportion S.

[0062] In some embodiments of the present embodiment, the above-mentioned drilled well facies belt division according to the three-dimensional seismic reflection characteristics and the preset type lithofacies proportion of the drilled well in a single volcanic eruption period comprises: drilled well facies belt division according to the volcanic sedimentary facies proportion of the drilled well in a single volcanic eruption period: S = volcanic sedimentary facies ÷ sum of all lithofacies; wherein S represents the volcanic sedimentary facies proportion, and S less than 50% is a proximal facies belt and S higher than 50% is a distal facies belt.

[0063] Specifically, the five lithofacies proportions of the drilled well in a single volcanic eruption period are used, combined with the three-dimensional seismic reflection characteristics, to judge the drilled well facies belt division according to the volcanic sedimentary facies proportion S, taking S = 50% as the benchmark, the volcanic sedimentary facies proportion S less than 50% is a proximal facies belt, and the volcanic sedimentary facies proportion S higher than 50% is a distal facies belt. Thus, the purpose of judging the drilled well facies belt is achieved.

[0064] S120: constructing a near-source facies belt pseudo well and a far-source facies belt pseudo well based on the drilled well facies belt division result and the characteristic curves corresponding to the multiple types of lithofacies of the drilled well in a single volcanic eruption period;

[0065] In some embodiments of the present embodiment, the method of identifying the boundaries of the volcanic mechanism facies belt further comprises: dividing the volcanic eruption period according to the drilled well drilling data; and performing lithofacies division according to the lithofacies in each volcanic eruption period to obtain a lithofacies division result; wherein the lithofacies division result includes explosive facies, overflow facies, extrusion facies, volcanic channel facies, and volcanic sedimentary facies.

[0066] For example, the drilled well drilling data includes well logging and coring data of the drilled well.

[0067] Specifically, the volcanic eruption period is divided according to the well logging and coring data of the drilled well, and the five types of lithofacies, i.e., explosive facies, overflow facies, extrusion facies, volcanic channel facies, and volcanic sedimentary facies, are divided according to the lithofacies in each volcanic eruption period.

[0068] In some embodiments of the present embodiment, the characteristic curves include acoustic curves and density curves; and the constructing of the near-source facies belt pseudo well and the far-source facies belt pseudo well based on the drilled well facies belt division result and the characteristic curves corresponding to the multiple types of lithofacies of the drilled well in a single volcanic eruption period comprises: determining the acoustic curves and the density curves of the explosive facies, the overflow facies, and the volcanic sedimentary facies based on the acoustic curves and the density curves of the drilled well and the lithofacies division result; and constructing the near-source facies belt pseudo well by using the acoustic curves and the density curves of the explosive facies and the overflow facies of the drilled well based on the drilled well facies belt division result, and constructing the far-source facies belt pseudo well by using the acoustic curves and the density curves of the volcanic sedimentary facies of the drilled well.

[0069] Specifically, the acoustic curves and the density curves of the explosive facies, the overflow facies, and the volcanic sedimentary facies are determined by using the acoustic curves and the density curves of the well logging of the drilled well and the lithofacies division result, and the pseudo well characteristics are spliced by using the acoustic curves and the density curves of the explosive facies, the overflow facies, and the volcanic sedimentary facies and the drilled well facies belt division result, so as to construct the a pseudo well and the b pseudo well. The a pseudo well and the b pseudo well have only one volcanic eruption period, the acoustic curves and the density curves of the a pseudo well in a single volcanic eruption period are represented by the acoustic curves and the density curves of the explosive facies and the overflow facies of the drilled well, and the a pseudo well is a near-source facies belt. The acoustic curves and the density curves of the b pseudo well in a single volcanic eruption period are represented by the acoustic curves and the density curves of the volcanic sedimentary facies of the drilled well, and the b pseudo well is a far-source facies belt. Thus, the purpose of constructing the near-source facies belt pseudo well and the far-source facies belt pseudo well is achieved.

[0070] S130: determining a value range corresponding to a preset type of lithofacies proportion of the volcanic mechanism based on the characteristic curves of the near-source facies belt pseudo well and the characteristic curves of the far-source facies belt pseudo well;

[0071] wherein the preset type lithofacies proportion is a volcanic sedimentary facies proportion.

[0072] In some embodiments of the present embodiment, the characteristic curve includes an acoustic curve and a density curve; and the determining of the value range corresponding to the preset type lithofacies proportion of the volcanic mechanism based on the characteristic curve of the near-source facies zone pseudo well and the characteristic curve of the far-source facies zone pseudo well includes: based on the acoustic curve and the density curve of the near-source facies zone pseudo well and the acoustic curve and the density curve of the far-source facies zone pseudo well, a two-dimensional seismic forward model is established by using a convolution model, ray tracing and wave equation method; based on the two-dimensional seismic forward model, an energy half-life time T curve and a thickness H curve in a corresponding volcanic eruption period are calculated; based on the energy half-life time T curve and the thickness H curve in the corresponding volcanic eruption period, an R curve is calculated by using a formula R = (T / H)*100, a value A corresponding to the volcanic sedimentary facies proportion in the R curve is analyzed, and a value range of the value A is obtained. Figure 5 is a curve diagram in a volcanic eruption period provided by an embodiment of the present application, as shown in Figure 5 the upper half of the figure is a two-dimensional seismic forward model of two pseudo wells, the left well in the figure is an actual drilled eruption facies well, the right well is a designed drilled volcanic sedimentary facies well, and the lower half of the figure is a curve extracted in a volcanic eruption period, light gray is an extracted energy half-life time T curve, and dark gray is an R curve obtained by dividing thickness by energy half-life time.

[0073] Specifically, by using the constructed acoustic curve and density curve of the a pseudo well and the acoustic curve and density curve of the b pseudo well, a two-dimensional seismic forward model is established by using a convolution model, ray tracing and wave equation method, and an energy half-life time T curve and a thickness H curve in a corresponding volcanic eruption period are calculated by using the two-dimensional seismic forward model (wherein the thickness H is obtained by subtracting a top surface from a bottom surface of the volcanic mechanism). An R curve is calculated by using a formula R = (T / H)*100, a value A corresponding to the volcanic sedimentary facies proportion in the R curve is analyzed, and a value range of the value A is the value range corresponding to the volcanic sedimentary facies proportion of the volcanic mechanism. The method considers the energy half-life time T curve and the thickness H curve, finds out the value A matching the volcanic sedimentary facies proportion on the R curve by constructing the R curve, and can more accurately reflect the value range corresponding to the volcanic sedimentary facies proportion of the volcanic mechanism.

[0074] It should be noted that the energy half-life time attribute is calculated in a given analysis time window to obtain a relative time position when the seismic reflection wave energy reaches half; the basic algorithm is to first calculate the cumulative amplitude sum of each sampling point in the given time window, then divide the obtained amplitude sum E by 2 to obtain the half energy Ef = E / 2 of the time window; the position of the half energy in the time window is analyzed, the number of sampling points n corresponding to the position is divided by the total number of sampling points m in the entire time window to obtain the energy half-life time T of the position point.

[0075] S140: obtaining a plane map representing facies belts based on the three-dimensional seismic data volume of the volcanic edifice;

[0076] Please refer to Figure 3 and Figure 4 , Figure 3 is a plane map of energy half-life time of a volcanic edifice provided by an embodiment of the present application. Figure 4 is a thickness map of a volcanic eruption period provided by an embodiment of the present application. In some embodiments of the present embodiment, the above-mentioned obtaining a plane map representing facies belts based on the three-dimensional seismic data volume of the volcanic edifice comprises: obtaining a plane map of energy half-life time t of a single volcanic eruption period based on the three-dimensional seismic data volume of the volcanic edifice, a bottom surface and a top surface of the volcanic edifice; calculating a thickness h plane map of the single volcanic eruption period based on the bottom surface and the top surface of the volcanic edifice of the single volcanic eruption period; and obtaining a plane map representing facies belts according to the formula r = (t / h)*100; wherein r is a plane map of a parameter representing facies belts.

[0077] Specifically, the bottom surface and the top surface of each volcanic eruption period are interpreted on the three-dimensional seismic data volume of the volcanic edifice, a plane map of energy half-life time t of a single volcanic eruption period is extracted, and a thickness plane h of the single volcanic eruption period is calculated. The plane map of energy half-life time t is divided by the thickness h plane map of the corresponding volcanic eruption period to obtain a plane map of a parameter r representing facies belts, r = (t / h)*100.

[0078] S150: performing volcanic edifice facies belt division on the plane map representing facies belts according to a value range corresponding to a preset type of lithofacies proportion of the volcanic edifice to obtain a boundary between a proximal facies belt and a distal facies belt within the corresponding volcanic eruption period.

[0079] In some embodiments of the present embodiment, the above-mentioned performing volcanic edifice facies belt division on the plane map representing facies belts according to a value range corresponding to a preset type of lithofacies proportion of the volcanic edifice to obtain a boundary between a proximal facies belt and a distal facies belt within the corresponding volcanic eruption period comprises: performing volcanic edifice facies belt division on the plane map representing facies belts based on a value A corresponding to the R curve of the volcanic sedimentary facies proportion to obtain a distal facies belt position and a proximal facies belt position of the volcanic edifice; wherein a region with a value greater than the value A on the plane map is taken as a distal facies belt of the volcanic edifice, and a region with a value less than the value A on the plane map is taken as a proximal facies belt of the volcanic edifice. The boundary between the proximal facies belt and the distal facies belt is obtained based on the distal facies belt position and the proximal facies belt position. As shown in Figure 6 Figure 6 ​is a plane distribution map of a volcanic mechanism facies belt provided by an embodiment of the present application, wherein the black boundary is the boundary between the near-source facies belt and the far-source facies belt of the volcanic mechanism. The thick black line is the range of A, and the inside of the thick black line (i.e. the range enclosed by the thick black line) is the near-source facies belt, and the outside of the thick black line (i.e. the range not enclosed by the thick black line) is the far-source facies belt.

[0080] In the implementation process, the method determines the drilled well facies belt division according to the three-dimensional seismic reflection characteristics and the volcanic sedimentary facies proportion of the drilled well in a single volcanic eruption period. Based on the drilled well facies belt division result and the characteristic curves corresponding to multiple types of lithofacies of the drilled well in a single volcanic eruption period, the near-source facies belt pseudo well and the far-source facies belt pseudo well are constructed. According to the characteristic curve of the near-source facies belt pseudo well and the characteristic curve of the far-source facies belt pseudo well, the value range corresponding to the volcanic sedimentary facies proportion of the volcanic mechanism is determined. According to the value range corresponding to the volcanic sedimentary facies proportion of the volcanic mechanism, the volcanic mechanism facies belt division is performed on the plane map representing the facies belt obtained by interpreting the three-dimensional seismic data volume of the volcanic mechanism, so as to obtain the boundary between the near-source facies belt and the far-source facies belt in the corresponding volcanic eruption period. Thus, the method can quantitatively, accurately and effectively define the boundary between the near-source facies belt and the far-source facies belt, quickly focus on the near-source favorable facies belt of the volcanic mechanism, and is beneficial to further searching for effective volcanic rock reservoirs, thereby laying a foundation for exploration targets.

[0081] Figure 7 is a block diagram of a system for identifying the facies belt boundary of a volcanic mechanism provided by an embodiment of the present application. As shown in Figure 7 , the embodiment of the present application provides a system for identifying the facies belt boundary of a volcanic mechanism, which comprises:

[0082] The drilled well facies belt division module is configured to perform drilled well facies belt division according to the three-dimensional seismic reflection characteristics and the preset type of lithofacies proportion of the drilled well in a single volcanic eruption period.

[0083] The pseudo well construction module is configured to construct the near-source facies belt pseudo well and the far-source facies belt pseudo well based on the drilled well facies belt division result and the characteristic curves corresponding to multiple types of lithofacies of the drilled well in a single volcanic eruption period.

[0084] The value range determination module is configured to determine the value range corresponding to the preset type of lithofacies proportion of the volcanic mechanism based on the characteristic curve of the near-source facies belt pseudo well and the characteristic curve of the far-source facies belt pseudo well.

[0085] The plane map obtaining module is configured to obtain a plane map representing the facies belt based on the three-dimensional seismic data volume of the volcanic mechanism.

[0086] The volcanic edifice facies belt division module is configured to divide the volcanic edifice facies belt on the planar graph representing the facies belt according to the value range corresponding to the volcanic sedimentary facies proportion of the volcanic edifice of the preset type, so as to obtain the boundary between the proximal facies belt and the distal facies belt in a single volcanic eruption period.

[0087] Specifically, the system determines the facies belt division of the drilled well according to the three-dimensional seismic reflection characteristics and the volcanic sedimentary facies proportion of the drilled well in a single volcanic eruption period. Based on the facies belt division result of the drilled well and the characteristic curves corresponding to the multiple types of lithofacies of the drilled well in a single volcanic eruption period, the proximal facies belt pseudo well and the distal facies belt pseudo well are constructed. According to the characteristic curve of the proximal facies belt pseudo well and the characteristic curve of the distal facies belt pseudo well, the value range corresponding to the volcanic sedimentary facies proportion of the volcanic edifice is determined. According to the value range corresponding to the volcanic sedimentary facies proportion of the volcanic edifice, the volcanic edifice facies belt division is performed on the planar graph representing the facies belt obtained by interpreting the three-dimensional seismic data volume of the volcanic edifice, so as to obtain the boundary between the proximal facies belt and the distal facies belt in a single volcanic eruption period. Therefore, the system can quantitatively, accurately and effectively define the boundary between the proximal facies belt and the distal facies belt, quickly focus on the proximal favorable facies belt of the volcanic edifice, and is beneficial to further search for the effective volcanic reservoir, thereby laying a foundation for exploration targets.

[0088] The embodiment of the present application also provides a machine readable storage medium, which stores instructions, and the instructions make the processor 100 be configured to execute the above-mentioned method for identifying the boundary of the volcanic edifice facies belt when the processor 100 executes the instructions.

[0089] The machine readable storage medium includes permanent and non-permanent, removable and non-removable media, and can be realized by any method or technology. The information can be computer readable instructions, data structures, program modules or other data. Examples of the storage medium of the computer 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 technologies, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassette tape, magnetic tape magnetic disk storage or other magnetic storage devices or any other non-transmission medium that can be used to store information accessible to a computing device. According to the definition in this paper, the computer readable medium does not include transitory computer readable media such as modulated data signals and carriers.

[0090] The embodiment of the present application further provides an electronic device 10, which comprises a memory 101, a processor 100, and a computer program 102 stored in the memory 101 and capable of running on the processor 100, and the processor 100 implements the method for identifying the boundary of the volcanic mechanism facies belt when executing the computer program 102.

[0091] As shown in Figure 8 Fig. 1 is a schematic diagram of an electronic device according to an embodiment of the present application. As shown in Figure 8 Fig. 1, the electronic device 10 of the embodiment comprises a processor 100, a memory 101, and a computer program 102 stored in the memory 101 and capable of running on the processor 100. The processor 100 implements the steps in the method embodiment when executing the computer program 102. Alternatively, the processor 100 implements the functions of the modules / units in the device embodiment when executing the computer program 102.

[0092] For example, the computer program 102 can be divided into one or more modules / units, which are stored in the memory 101 and executed by the processor 100 to complete the present application. The one or more modules / units can be a series of computer program instruction segments capable of completing a specific function, which are used to describe the execution process of the computer program 102 in the electronic device 10. For example, the computer program 102 can be divided into a drilled well facies belt division module, a pseudo well construction module, a value range determination module, a plan view obtaining module, and a volcanic mechanism facies belt division module.

[0093] The electronic device 10 can be a desktop computer, a notebook computer, a palm computer, a cloud server, and the like. The electronic device 10 can comprise, but is not limited to, the processor 100 and the memory 101. Those skilled in the art can understand that Figure 8 The electronic device 10 is only an example and does not constitute a limitation on the electronic device 10, and can comprise more or fewer components than those shown, or combine certain components, or different components, for example, the electronic device can further comprise an input / output device, a network access device, a bus, and the like.

[0094] The processor 100 can be a central processing unit (CPU), and can also be other general-purpose processors, a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or the like. The general-purpose processor can be a microprocessor, or the processor can also be any conventional processor.

[0095] The memory 101 can be an internal storage unit of the electronic device 10, such as a hard disk or a memory of the electronic device 10. The memory 101 can also be an external storage device of the electronic device 10, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, or the like equipped on the electronic device 10. Further, the memory 101 can include both the internal storage unit and the external storage device of the electronic device 10. The memory 101 is used to store computer programs and other programs and data required by the electronic device 10. The memory 101 can also be used to temporarily store data that has been output or will be output.

[0096] Those skilled in the art can clearly understand that, for the convenience and brevity of description, only the division of the above functional units and modules is taken as an example for illustration, and in actual application, the above functions can be completed by different functional units and modules according to needs, that is, the internal structure of the device is divided into different functional units or modules to complete all or part of the functions described above. Each functional unit and module in the embodiment can be integrated in one processing unit, or each unit can exist physically, or two or more units can be integrated in one unit, and the integrated unit can be realized in the form of hardware or in the form of software functional unit. In addition, the specific names of each functional unit and module are only for easy distinction, and do not limit the protection scope of the present application. The specific working process of the units and modules in the above system can refer to the corresponding process in the foregoing method embodiments, which will not be described here.

[0097] Those skilled in the art will appreciate that embodiments of the application can be readily used as a method, a system or an article of manufacture, such as a computer program product 102. Accordingly, embodiments of the application can take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, embodiments of the application can take the form of a computer program product 102 on one or more computer-usable storage media (including, but not limited to, magnetic disks, CD-ROMs, optical storage media such as DVD s, etc.) having computer usable program code embodied thereon.

[0098] The present application is described in reference to the flow diagrams and / or block diagrams of the methods, apparatus (systems) and computer program products 102 according to embodiments of the application. It will be understood that each block of the flow diagrams and / or block diagrams, and combinations of blocks in the flow diagrams and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor 100 of a general purpose computer, special purpose computer, embedded processing machine, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor 100 of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flow diagrams and / or block diagrams block or blocks. Figure 1 one or more functions specified in one or more of the flow diagram or block diagram blocks. Figure 1 one or more functions specified in one or more of the flow diagram or block diagram blocks.

[0099] These computer program instructions can also be stored in a computer- readable memory 101 that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory 101 produce an article of manufacture including instructions which implement the flow diagrams and / or block diagrams block or blocks. Figure 1 one or more functions specified in one or more of the flow diagram or block diagram blocks. Figure 1 one or more functions specified in one or more of the flow diagram or block diagram blocks.

[0100] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the flow diagrams and / or block diagrams block or blocks. Figure 1 one or more functions specified in one or more of the flow diagram or block diagram blocks. Figure 1 one or more functions specified in one or more of the flow diagram or block diagram blocks.

[0101] It should also be noted that the terms "comprising", "comprises" or other variations thereof 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 can also include other elements not expressly listed or inherent to such process, method, article, or apparatus.

[0102] The above embodiments are only used to illustrate the present application, but not to limit it. Instead of the above, various modifications and changes can be made to the application by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the application shall fall into the scope of the claims of the application.

Claims

1. A method of identifying a volcanic edifice facies belt boundary, characterized by, The method comprises the following steps: According to the three-dimensional seismic reflection characteristics and the preset type lithofacies proportion of the drilled well in a single volcanic eruption period, the drilled well facies belt is divided; Based on the drilled well facies belt division result and the characteristic curve corresponding to the multiple type lithofacies of the drilled well in a single volcanic eruption period, the near-source facies belt pseudo well and the far-source facies belt pseudo well are constructed; Based on the characteristic curve of the near-source facies belt pseudo well and the characteristic curve of the far-source facies belt pseudo well, the value range corresponding to the preset type lithofacies proportion of the volcanic mechanism is determined; Based on the three-dimensional seismic data volume of the volcanic mechanism, the plane map representing the facies belt is obtained; According to the value range corresponding to the preset type lithofacies proportion of the volcanic mechanism, the volcanic mechanism facies belt division is carried out on the plane map representing the facies belt, and the boundary between the near-source facies belt and the far-source facies belt in the corresponding volcanic eruption period is obtained.

2. The method of identifying a volcanic edifice facies belt boundary according to claim 1, wherein, Further comprising: According to the drilling data of the drilled well, the volcanic eruption period is divided; According to the lithofacies in each volcanic eruption period, the lithofacies division result is obtained; The lithofacies division result includes explosive facies, overflow facies, outburst facies, volcanic channel facies and volcanic sedimentary facies.

3. The method of identifying a volcanic edifice facies belt boundary according to claim 2, wherein, The characteristic curve includes acoustic curve and density curve; The construction of the near-source facies belt pseudo well and the far-source facies belt pseudo well based on the drilled well facies belt division result and the characteristic curve corresponding to the multiple type lithofacies of the drilled well in a single volcanic eruption period comprises: Based on the acoustic curve and the density curve of the drilled well and the lithofacies division result, the acoustic curve and the density curve of the explosive facies, the overflow facies and the volcanic sedimentary facies are determined; Based on the drilled well facies belt division result, the near-source facies belt pseudo well is constructed by using the acoustic curve and the density curve of the explosive facies and the overflow facies of the drilled well, and the far-source facies belt pseudo well is constructed by using the acoustic curve and the density curve of the volcanic sedimentary facies of the drilled well.

4. The method of identifying a volcanic edifice facies belt boundary according to claim 2, wherein, The drilled well facies belt division according to the three-dimensional seismic reflection characteristics and the preset type lithofacies proportion of the drilled well in a single volcanic eruption period comprises: According to the three-dimensional seismic reflection characteristics and the volcanic sedimentary facies proportion of the drilled well in a single volcanic eruption period, the drilled well facies belt is divided: S=volcanic sedimentary facies÷sum of all lithofacies; Wherein, S represents the volcanic sedimentary facies proportion, S lower than 50% is the near-source facies belt, and S higher than 50% is the far-source facies belt.

5. The method of identifying a volcanic edifice facies belt boundary of claim 1, wherein, The characteristic curve includes acoustic curve and density curve; The determination of the value range corresponding to the preset type lithofacies proportion of the volcanic mechanism based on the characteristic curve of the near-source facies belt pseudo well and the characteristic curve of the far-source facies belt pseudo well comprises: Based on the acoustic curve and the density curve of the near-source facies belt pseudo well and the acoustic curve and the density curve of the far-source facies belt pseudo well, a two-dimensional seismic forward model is established by using convolution model, ray tracing and wave equation method; Based on the two-dimensional seismic forward model, the energy half-life time T curve and the thickness H curve corresponding to the volcanic eruption period are calculated; Based on the energy half-life time T curve and the thickness H curve corresponding to the volcanic eruption period, the R curve is calculated by using the formula R=(T / H)*100, the value A corresponding to the volcanic sedimentary facies proportion in the R curve is analyzed, and the value range of the value A is obtained.

6. The method of identifying a volcanic edifice facies belt boundary according to claim 5, wherein, The value range corresponding to the preset type of lithofacies proportion of the volcanic edifice is used for volcanic edifice facies belt division on the plane map representing the facies belt, to obtain a boundary between a proximal facies belt and a distal facies belt in a corresponding volcanic eruption period, and the boundary includes: Based on the value A corresponding to the volcanic sedimentary facies proportion on the R curve, the volcanic edifice facies belt is divided on the plane map representing the facies belt, to obtain a distal facies belt position and a proximal facies belt position of the volcanic edifice; Based on the distal facies belt position and the proximal facies belt position, a boundary between the proximal facies belt and the distal facies belt is obtained.

7. The method of identifying a volcanic edifice facies belt boundary according to claim 6, wherein, The volcanic edifice facies belt is divided on the plane map representing the facies belt based on the value A corresponding to the volcanic sedimentary facies proportion on the R curve, and the division includes: The area on the plane map with a value greater than the value A is taken as a distal facies belt of the volcanic edifice, and the area with a value less than the value A is taken as a proximal facies belt of the volcanic edifice.

8. The method of identifying a volcanic edifice facies belt boundary of claim 1, wherein, The plane map representing the facies belt is obtained based on the three-dimensional seismic data volume of the volcanic edifice, and the obtaining includes: Based on the three-dimensional seismic data volume of the volcanic edifice, an energy half-life time plane t of a single volcanic eruption period, a bottom surface and a top surface of the volcanic edifice are obtained; Based on the bottom surface and the top surface of the single volcanic eruption period, a thickness h plane of the single volcanic eruption period is calculated; According to a formula r = (t / h)*100, the plane map representing the facies belt is obtained; wherein r is a plane map of a parameter corresponding to the facies belt.

9. A system for identifying a volcanic edifice facies belt boundary, characterized by, The method includes: A drilled well facies belt division module is configured to divide drilled well facies belts according to three-dimensional seismic reflection characteristics and a preset type of lithofacies proportion of a single volcanic eruption period of drilled wells; A pseudo-well construction module is configured to construct a proximal facies belt pseudo-well and a distal facies belt pseudo-well based on drilled well facies belt division results and characteristic curves of multiple types of lithofacies of a single volcanic eruption period of drilled wells; A value range determination module is configured to determine a value range corresponding to a preset type of lithofacies proportion of a volcanic edifice based on characteristic curves of the proximal facies belt pseudo-well and the distal facies belt pseudo-well; A plane map obtaining module is configured to obtain a plane map representing the facies belt based on a three-dimensional seismic data volume of the volcanic edifice; A volcanic edifice facies belt division module is configured to divide volcanic edifice facies belts on the plane map representing the facies belt according to the value range corresponding to the preset type of lithofacies proportion of the volcanic edifice, to obtain a boundary between a proximal facies belt and a distal facies belt in a corresponding volcanic eruption period.

10. A machine-readable storage medium having stored thereon instructions, the instructions being executable by a machine to cause the machine to: The instructions, when executed by the processor, cause the processor to be configured to perform the method of identifying a volcanic edifice facies belt boundary according to any one of claims 1 to 8.

11. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, The processor performs the computer program to implement the method of identifying a volcanic edifice facies belt boundary according to any one of claims 1 to 8.

Images