Method and system for identifying volcanic edifice facies belt boundaries
By constructing pseudo-well and forward modeling models using three-dimensional seismic reflection characteristics and drilled well data, the accuracy problem of identifying volcanic structure facies boundaries was solved, and the quantitative identification and focusing of effective reservoirs was achieved.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- CHINA PETROLEUM & CHEMICAL CORP
- Filing Date
- 2024-04-22
- Publication Date
- 2026-07-21
AI Technical Summary
Existing technologies cannot accurately and effectively define the boundary between the near-source and far-source facies zones of volcanic structures, making it difficult to quantitatively identify effective reservoirs.
By using three-dimensional seismic reflection characteristics and the lithofacies proportion of drilled wells within a single volcanic eruption period, pseudo-wells of near-source and far-source facies zones are constructed. A two-dimensional seismic forward model is established using sonic curves and density curves to calculate energy half-decay time and thickness, determine the range of volcanic sedimentary facies proportion, and delineate facies zone boundaries.
It enables quantitative and accurate definition of the boundary between near-source and far-source facies zones, and quickly focuses on favorable near-source facies zones of volcanic structures, helping to find effective reservoirs.
Smart Images
Figure CN120831708B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of petroleum exploration and development technology, and specifically to a method and a system for identifying the boundaries of volcanic facies zones. Background Technology
[0002] Volcanic facies are generally classified into five major facies based on their lithology, rock fabric, and genesis: eruptive facies, effusive facies, volcanic conduit facies, and volcanic sedimentary facies. Identifying these five facies in actual seismic events presents significant challenges. In volcanic rock exploration and development, volcanic structures are typically divided into three facies zones based on their distance from the crater: crater-near-crater, near-source, and far-source. In the Songnan fault depression, the distribution of intermediate-basic volcanic gas reservoirs in the Huoshiling Formation is primarily controlled by effective reservoirs. The near-crater-near-source facies zone exhibits more developed primary porosity, a larger proportion of high-porosity and permeability zones, and higher gas concentrations. Therefore, identifying the boundary between the near-source and far-source facies zones is crucial for finding effective reservoirs. Currently, the near-source facies zone is characterized mainly by analyzing seismic reflection characteristics and geometric properties such as stratigraphic dip and coherence to determine the crater's location, thus qualitatively approximating the distribution of the near-source facies zone. However, there is still no quantitative method for predicting the boundary between the near-source and far-source facies zones of a volcanic structure.
[0003] Energy half-decay properties primarily reflect changes in seismic wave energy and are generally used to indicate variations in sedimentary environments and lithological facies. They have been applied in the identification of tight sandstones such as turbidites, but there is no evidence to prove their application in the identification of volcanic rock facies and facies zones. Due to the high volcanic velocity resulting in strong seismic reflections enveloping the volcanic structure, and the chaotic and weak reflections of eruptive facies, there are certain differences in energy half-decay times. Furthermore, during volcanic eruptions, volcanic material accumulates near the crater, resulting in thick strata near the crater. Therefore, combining these two factors can serve as a primary basis for determining the boundaries between near-source and far-source facies zones. However, existing methods can only extrapolate seismic reflection characteristics from actual drilling and combine seismic geometric properties to locate the crater and determine the approximate distribution, failing to effectively define the planar boundaries between near-source and far-source facies zones of volcanic structures. Summary of the Invention
[0004] The purpose of this invention is to provide a method and system for identifying the boundaries of facies zones in volcanic structures, so as to at least solve the problem mentioned above of the inability to accurately and effectively define the boundaries between the near-source and far-source facies zones of volcanic structures.
[0005] To achieve the above objectives, a first aspect of the present invention provides a method for identifying the boundaries of volcanic facies zones, comprising:
[0006] Based on the three-dimensional seismic reflection characteristics and the proportion of pre-defined lithofacies types in drilled wells within a single volcanic eruption period, the drilled well facies zones are divided.
[0007] Based on the results of the facies zone division of drilled wells and the characteristic curves corresponding to various rock types in a single volcanic eruption, pseudo-wells of near-source facies zone and pseudo-wells of far-source facies zone are constructed.
[0008] Based on the characteristic curves of pseudo-wells in the near-source facies zone and pseudo-wells in the far-source facies zone, the range of values corresponding to the pre-defined type of lithofacies proportion of volcanic structures is determined.
[0009] Based on the 3D seismic data volume of the volcanic structure, a planar map representing the facies zone is obtained;
[0010] Based on the range of values corresponding to the pre-defined type of facies proportion of the volcanic structure, the facies zones of the volcanic structure are divided on the plane map representing the facies zones, and the boundary between the near-source facies zone and the far-source facies zone within the corresponding volcanic eruption period is obtained.
[0011] Optionally, the above methods for identifying the boundaries of volcanic facies zones also include:
[0012] Based on the drilling data, the volcanic eruption phases are divided;
[0013] The lithofacies were classified according to the lithofacies of each volcanic eruption period, and the lithofacies classification results were obtained.
[0014] The lithofacies classification results include explosive facies, overflow facies, intrusive facies, volcanic conduit facies, and volcanic sedimentary facies.
[0015] Optionally, the aforementioned characteristic curves include acoustic curves and density curves;
[0016] Based on the facies zone division results of drilled wells and the characteristic curves corresponding to various rock types in a single volcanic eruption, pseudo-wells of near-source and far-source facies zones are constructed, including:
[0017] Based on the sonic and density curves of the drilled wells and the lithofacies classification results, the sonic and density curves of the eruptive facies, the overflow facies, and the volcanic sedimentary facies were determined.
[0018] Based on the facies zone division results of the drilled wells, pseudo-wells of the near-source facies zone are constructed using the sonic curves and density curves of the eruptive and overflow facies zones of the drilled wells, and pseudo-wells of the distant-source facies zone are constructed using the sonic curves and density curves of the volcanic sedimentary facies zones of the drilled wells.
[0019] Optionally, the above-mentioned well facies zone division is based on three-dimensional seismic reflection characteristics and the proportion of pre-defined lithofacies types in drilled wells within a single volcanic eruption period, including:
[0020] Based on three-dimensional seismic reflection characteristics and the proportion of volcanic sedimentary facies in drilled wells within a single volcanic eruption period, the drilled well facies zones are divided as follows:
[0021] S = volcanic sedimentary facies ÷ sum of all lithofacies;
[0022] Where S represents the proportion of volcanic sedimentary facies; S below 50% indicates a near-source facies zone, and S above 50% indicates a far-source facies zone.
[0023] Optionally, the aforementioned characteristic curves include acoustic curves and density curves;
[0024] The aforementioned characteristic curves of near-source facies pseudo-wells and far-source facies pseudo-wells determine the value range corresponding to the preset type lithofacies proportion of volcanic structures, including:
[0025] Based on the acoustic and density curves of pseudo-wells in the near-source facies zone and the far-source facies zone, a two-dimensional seismic forward model was established using the convolution model, ray tracing, and wave equation methods.
[0026] Based on the two-dimensional seismic forward model, the energy half-life T curve and thickness H curve within the corresponding volcanic eruption period are calculated.
[0027] Based on the energy half-life T curve and thickness H curve within the corresponding volcanic eruption period, the R curve is calculated using the formula R=(T / H)*100. The value A corresponding to the proportion of volcanic sedimentary phases in the R curve is analyzed, and the range of values for A is obtained.
[0028] Optionally, based on the value range corresponding to the pre-defined type of lithofacies proportion of the volcanic structure, the volcanic structure facies zones are divided on the characterizing facies map to obtain the boundary between the proximal and distal facies zones within the corresponding volcanic eruption period, including:
[0029] Based on the value A corresponding to the proportion of volcanic sedimentary facies on the R curve, the volcanic structure facies zones are divided on the planar map characterizing the facies zones to obtain the location of the distal and proximal facies zones of the volcanic structure.
[0030] Based on the location of the distant source phasor and the near source phasor, the boundary between the near source phasor and the distant source phasor is obtained.
[0031] Optionally, the above-mentioned volcanic facies zone division on the planar map characterizing facies zones, based on the value A corresponding to the proportion of volcanic sedimentary facies in the R curve, includes:
[0032] The region on the planar map with values greater than A is designated as the distal phalanx of the volcanic structure, and the region with values less than A is designated as the proximal phalanx of the volcanic structure.
[0033] Optionally, the above-mentioned 3D seismic data volume based on volcanic structure yields a planar map characterizing the facies zone, including:
[0034] Based on the three-dimensional seismic data volume of the volcanic structure, the energy half-decay planar diagram t, the bottom and top surfaces of the volcanic structure for a single volcanic eruption are obtained;
[0035] Based on the bottom and top surfaces of the volcanic structure of a single eruption, calculate the thickness h planar diagram of a single eruption.
[0036] According to the formula r=(t / h)*100, the planar diagram representing the phase band is obtained; where r is the planar diagram representing the corresponding parameters of the phase band.
[0037] A second aspect of the present invention provides a system for identifying the boundaries of volcanic facies zones, comprising:
[0038] The drilled facies zone delineation module is used to delineate drilled facies zones based on three-dimensional seismic reflection characteristics and the proportion of pre-defined lithofacies types in drilled wells within a single volcanic eruption period.
[0039] The pseudo-well construction module is used to construct near-source facies zone pseudo-wells and far-source facies zone pseudo-wells based on the facies zone division results of drilled wells and the characteristic curves corresponding to various rock types in a single volcanic eruption period of drilled wells.
[0040] The value range determination module is used to determine the value range corresponding to the pre-defined type facies proportion of volcanic structures based on the characteristic curves of near-source facies pseudo-wells and far-source facies pseudo-wells.
[0041] The planar plot module is used to obtain planar plots representing facies zones from 3D seismic data volumes based on volcanic mechanisms.
[0042] The volcanic structure facies zone division module is used to divide the volcanic structure facies zones on a planar map representing facies zones according to the value range corresponding to the preset type of rock facies proportion of the volcanic structure, and obtain the boundary between the near-source facies zone and the far-source facies zone within the corresponding volcanic eruption period.
[0043] In a third aspect, the present invention provides a machine-readable storage medium storing instructions that, when executed by a processor, configure the processor to perform the method for identifying the boundary of a volcanic structure phase zone.
[0044] In a fourth aspect, an electronic device is provided, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the method for identifying the boundary of a volcanic phalanx.
[0045] The above technical solution provides a method and system for identifying facies boundaries of volcanic structures. Based on three-dimensional seismic reflection characteristics and the proportion of volcanic sedimentary facies in drilled wells within a single volcanic eruption period, the facies zone division of drilled wells is determined. Based on the facies zone division results of drilled wells and the corresponding characteristic curves of various rock types within a single volcanic eruption period, pseudo-wells of near-source and far-source facies zones are constructed. According to the characteristic curves of the near-source and far-source facies zone pseudo-wells, the range of values corresponding to the proportion of volcanic sedimentary facies in the volcanic structure is determined. Based on the range of values corresponding to the proportion of volcanic sedimentary facies in the volcanic structure, facies zone division of the volcanic structure is performed on a planar map representing the facies zones obtained by interpreting the three-dimensional seismic data of the volcanic structure, thus obtaining the boundary between the near-source and far-source facies zones within the corresponding volcanic eruption period. Therefore, this method and system can quantitatively, accurately, and effectively define the boundary between near-source and far-source facies zones, quickly focusing on favorable near-source facies zones of volcanic structures, which is beneficial for further searching for effective volcanic reservoirs, thereby laying the foundation for exploration targets.
[0046] Other features and advantages of the embodiments of the present invention will be described in detail in the following detailed description section. Attached Figure Description
[0047] The accompanying drawings are provided to further illustrate embodiments of the present invention and form part of the specification. They are used together with the following detailed description to explain the embodiments of the present invention, but do not constitute a limitation thereof. In the drawings:
[0048] Figure 1 This is a flowchart of a method for identifying the boundary of a volcanic facies zone according to one embodiment of the present invention;
[0049] Figure 2 This is a flowchart of another method for identifying the boundary of a volcanic structure phase zone provided by one embodiment of the present invention;
[0050] Figure 3 This is a planar diagram of the energy half-decay of a volcanic mechanism provided in one embodiment of the present invention;
[0051] Figure 4 This is a thickness map of volcanic eruption phases provided by one embodiment of the present invention;
[0052] Figure 5 This is a schematic diagram of the curves during a volcanic eruption period provided by one embodiment of the present invention;
[0053] Figure 6 This is a planar distribution diagram of volcanic facies zone division provided by one embodiment of the present invention;
[0054] Figure 7This is a block diagram of a system for identifying the boundary of a volcanic facies zone according to one embodiment of the present invention;
[0055] Figure 8 This is a schematic diagram of an electronic device structure provided by a preferred embodiment of the present invention.
[0056] Explanation of reference numerals in the attached figures
[0057] 10 - Electronic device, 100 - Processor, 101 - Memory, 102 - Computer program. Detailed Implementation
[0058] The specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are for illustration and explanation only and are not intended to limit the present invention.
[0059] Figure 1 This is a flowchart of a method for identifying the boundary of a volcanic facies zone according to one embodiment of the present invention. Figure 2 This is a flowchart of another method for identifying volcanic phagocytic zone boundaries provided by one embodiment of the present invention. Figure 1 and Figure 2 As shown, an embodiment of the present invention provides a method for identifying the boundary of a volcanic facies zone, comprising:
[0060] S110: Based on the three-dimensional seismic reflection characteristics and the proportion of pre-defined lithofacies types in drilled wells within a single volcanic eruption period, the drilled well facies zones are divided.
[0061] Among them, the proportion of the preset type of lithofacies is the proportion of volcanic sedimentary facies S.
[0062] In some embodiments of this example, the above-mentioned division of drilled facies zones based on three-dimensional seismic reflection characteristics and the proportion of pre-defined lithofacies types in drilled wells within a single volcanic eruption period includes: dividing drilled facies zones based on three-dimensional seismic reflection characteristics and the proportion of volcanic sedimentary facies in drilled wells within a single volcanic eruption period: S = volcanic sedimentary facies ÷ sum of all lithofacies; where S represents the proportion of volcanic sedimentary facies, S less than 50% is a near-source facies zone, and S greater than 50% is a far-source facies zone.
[0063] Specifically, by utilizing the proportions of five lithofacies in drilled wells during a single volcanic eruption and combining this with three-dimensional seismic reflection characteristics, the facies zones of drilled wells are determined based on the proportion S of volcanic sedimentary facies. Using S = 50% as a benchmark, a proportion S below 50% indicates a near-source facies zone, while a proportion S above 50% indicates a distant-source facies zone. This achieves the purpose of determining the facies zones of drilled wells.
[0064] S120: Based on the results of the facies zone division of drilled wells and the characteristic curves corresponding to various rock types in a single volcanic eruption, pseudo-wells of near-source facies zone and pseudo-wells of far-source facies zone are constructed.
[0065] In some embodiments of this example, the method for identifying the boundaries of volcanic facies zones further includes: dividing volcanic eruption periods based on drilled well data; dividing lithofacies based on the lithofacies within each eruption period to obtain lithofacies division results; wherein the lithofacies division results include eruptive facies, effusive facies, intrusive facies, volcanic conduit facies, and volcanic sedimentary facies.
[0066] For example, the aforementioned drilled well data includes drilled well logging and core data.
[0067] Specifically, based on the drilling, logging and coring data, the volcanic eruption periods are divided into five lithofacies: eruptive facies, effusive facies, intrusive facies, volcanic conduit facies, and volcanic sedimentary facies.
[0068] In some embodiments of this example, the aforementioned characteristic curves include acoustic curves and density curves; the construction of near-source facies zone pseudo-wells and far-source facies zone pseudo-wells based on the facies zone division results of drilled wells and the characteristic curves corresponding to various rock types in a single volcanic eruption period of drilled wells includes: determining the acoustic curves and density curves of eruptive facies, effusive facies, and volcanic sedimentary facies based on the acoustic curves and density curves of drilled wells and the facies division results; constructing near-source facies zone pseudo-wells using the acoustic curves and density curves of eruptive facies and effusive facies of drilled wells based on the facies zone division results of drilled wells, and constructing far-source facies zone pseudo-wells using the acoustic curves and density curves of volcanic sedimentary facies of drilled wells.
[0069] Specifically, using the sonic and density curves from drilled wells, along with lithofacies classification results, the sonic and density curves of the eruptive, effusive, and volcanic sedimentary facies were identified. Then, using these curves, along with the existing well facies zone classification results, pseudo-well characteristics were pieced together to construct pseudo-wells a and b. Pseudo-wells a and b each represent only one volcanic eruption period. Within this single eruption period, the sonic and density curves of pseudo-well a are represented by the sonic and density curves of the drilled eruptive and effusive facies, indicating pseudo-well a is a near-source facies zone. Similarly, the sonic and density curves of pseudo-well b within the same eruption period are represented by the sonic and density curves of the drilled volcanic sedimentary facies, indicating pseudo-well b is a far-source facies zone. This achieves the goal of constructing near-source and far-source facies zone pseudo-wells.
[0070] S130: Based on the characteristic curves of near-source facies zone pseudo-wells and far-source facies zone pseudo-wells, determine the value range corresponding to the preset type lithofacies proportion of volcanic structures.
[0071] The pre-defined type of lithofacies proportion is the proportion of volcanic sedimentary facies.
[0072] In some embodiments of this example, the aforementioned characteristic curves include acoustic curves and density curves; the determination of the value range corresponding to the pre-defined type facies proportion of the volcanic structure based on the characteristic curves of near-source facies pseudo-wells and far-source facies pseudo-wells includes: establishing a two-dimensional seismic forward model using convolution model, ray tracing, and wave equation methods based on the acoustic and density curves of near-source facies pseudo-wells and far-source facies pseudo-wells; calculating the energy half-decay time T curve and thickness H curve within the corresponding volcanic eruption period based on the two-dimensional seismic forward model; calculating the R curve using the formula R = (T / H) * 100 based on the energy half-decay time T curve and thickness H curve within the corresponding volcanic eruption period, analyzing the value A corresponding to the volcanic sedimentary facies proportion on the R curve, and obtaining the value range of value A. Figure 5 This is a schematic diagram of the curves during a volcanic eruption period provided by one embodiment of the present invention, such as... Figure 5 As shown, the upper half of the figure is a two-dimensional seismic forward model of two pseudo-wells. The well on the left is the well that actually encountered the eruptive facies, and the well on the right is the well that was designed to encounter the volcanic sedimentary facies. The lower half of the figure shows the curves extracted during a volcanic eruption period. The light gray curve is the T curve with the extracted energy half-decay, and the dark gray curve is the R curve obtained by dividing the thickness by the energy half-decay.
[0073] Specifically, using the sonic and density curves of the constructed pseudo-wells a and b, a two-dimensional seismic forward model was established using convolution, ray tracing, and wave equation methods. This model was then used to calculate the energy half-decay time (T) curve and thickness (H) curve for the corresponding volcanic eruption period (where thickness H is obtained by subtracting the top surface from the bottom surface of the volcanic structure). The R curve was calculated using the formula R = (T / H) * 100. The value A corresponding to the proportion of volcanic sedimentary facies on the R curve was analyzed; the range of values for A represents the range of values corresponding to the proportion of volcanic sedimentary facies in the volcanic structure. This method considers both the energy half-decay time (T) curve and the thickness (H) curve. By constructing the R curve and finding the value A on it that matches the proportion of volcanic sedimentary facies, it can more accurately reflect the range of values corresponding to the proportion of volcanic sedimentary facies in the volcanic structure.
[0074] It should be noted that the energy half-decay attribute is obtained by calculating the relative time position when the seismic reflected wave energy reaches half within a given analysis time window. The basic algorithm is to first calculate the cumulative amplitude of each sampling point within the given time window, then divide the obtained amplitude sum E by 2 to obtain the half-energy of the time window Ef = E / 2; analyze the time window position corresponding to the half-energy, and divide the number of sampling points n corresponding to that position by the total number of sampling points m in the entire time window to obtain the energy half-decay time T = n / m at that position.
[0075] S140: A three-dimensional seismic data volume based on volcanic structure, yielding a planar map characterizing the facies zone;
[0076] Please refer to Figure 3 and Figure 4 , Figure 3 This is a planar diagram of the energy half-decay of a volcanic mechanism provided in one embodiment of the present invention. Figure 4 This is a thickness map of a volcanic eruption phase provided by one embodiment of the present invention. In some embodiments of this embodiment, the above-mentioned three-dimensional seismic data volume based on the volcanic structure is used to obtain a planar map representing the phase zone, including: obtaining a planar map t of the energy half-decay time of a single volcanic eruption phase, the bottom surface and the top surface of the volcanic structure based on the three-dimensional seismic data volume based on the volcanic structure; calculating a planar map of the thickness h of a single volcanic eruption phase based on the bottom surface and the top surface of the volcanic structure based on the single volcanic eruption phase; and obtaining a planar map representing the phase zone according to the formula r = (t / h) * 100; where r is a planar map representing the parameters corresponding to the phase zone.
[0077] Specifically, the bottom and top surfaces of each volcanic eruption are interpreted on the three-dimensional seismic data volume of the volcanic structure. The energy half-decay plane diagram t of a single volcanic eruption is extracted, and the thickness plane h of a single volcanic eruption is calculated. The energy half-decay plane diagram t is divided by the thickness plane diagram h of the corresponding volcanic eruption to obtain the plane diagram of the parameter r representing the phase zone, r = (t / h)*100.
[0078] S150: Based on the range of values corresponding to the pre-defined type of facies proportion of the volcanic structure, the facies zones of the volcanic structure are divided on the plane map representing the facies zones, and the boundary between the near-source facies zone and the far-source facies zone within the corresponding volcanic eruption period is obtained.
[0079] In some embodiments of this example, the above-mentioned method of dividing volcanic facies zones on a planar map representing facies zones based on the value range corresponding to the preset type facies proportion of the volcanic structure, to obtain the boundary between the near-source and far-source facies zones within the corresponding volcanic eruption period, includes: dividing volcanic facies zones on a planar map representing facies zones based on the value A corresponding to the volcanic sedimentary facies proportion on the R curve, to obtain the location of the far-source and near-source facies zones of the volcanic structure; wherein, the area on the planar map with a value greater than A is taken as the far-source facies zone of the volcanic structure, and the area with a value less than A is taken as the near-source facies zone of the volcanic structure. Based on the location of the far-source and near-source facies zones, the boundary between the near-source and far-source facies zones is obtained. Figure 6 As shown, Figure 6This is a planar distribution diagram of volcanic facies zones provided by one embodiment of the present invention, wherein the black boundary is the boundary between the near-source facies zone and the far-source facies zone of the volcanic structure. The thick black line represents the range A, the area inside the thick black line (i.e., the area enclosed by the thick black line) is the near-source facies zone, and the area outside the thick black line (i.e., the area not enclosed by the thick black line) is the far-source facies zone.
[0080] In the above implementation process, this method determines the facies zone division of drilled wells based on 3D seismic reflection characteristics and the proportion of volcanic sedimentary facies in a single volcanic eruption period. Based on the facies zone division results of drilled wells and the corresponding characteristic curves of various rock types within a single volcanic eruption period, pseudo-wells of near-source and far-source facies zones are constructed. According to the characteristic curves of the near-source and far-source facies zone pseudo-wells, the range of values corresponding to the proportion of volcanic sedimentary facies in the volcanic structure is determined. Based on the range of values corresponding to the proportion of volcanic sedimentary facies in the volcanic structure, facies zone division of the volcanic structure is performed on a planar map representing the facies zones obtained by interpreting the 3D seismic data of the volcanic structure, thus obtaining the boundary between the near-source and far-source facies zones within the corresponding volcanic eruption period. Therefore, this method can quantitatively, accurately, and effectively define the boundary between near-source and far-source facies zones, quickly focusing on favorable near-source facies zones of the volcanic structure, which is beneficial for further searching for effective volcanic reservoirs, thereby laying the foundation for exploration targets.
[0081] Figure 7 This is a block diagram of a system for identifying the boundaries of volcanic facies zones according to one embodiment of the present invention. Figure 7 As shown, an embodiment of the present invention provides a system for identifying the boundary of a volcanic phagocytic structure, comprising:
[0082] The drilled facies zone delineation module is used to delineate drilled facies zones based on three-dimensional seismic reflection characteristics and the proportion of pre-defined lithofacies types in drilled wells within a single volcanic eruption period.
[0083] The pseudo-well construction module is used to construct near-source facies zone pseudo-wells and far-source facies zone pseudo-wells based on the facies zone division results of drilled wells and the characteristic curves corresponding to various types of rocks in a single volcanic eruption period of drilled wells.
[0084] The value range determination module is used to determine the value range corresponding to the pre-defined type facies proportion of volcanic structures based on the characteristic curves of near-source facies pseudo-wells and far-source facies pseudo-wells.
[0085] The planar plot module is used to obtain planar plots representing facies zones from 3D seismic data volumes based on volcanic mechanisms.
[0086] The volcanic structure facies zone division module is used to divide the volcanic structure facies zones on a planar map representing facies zones according to the value range corresponding to the preset type of rock facies proportion of the volcanic structure, and obtain the boundary between the near-source facies zone and the far-source facies zone within the corresponding volcanic eruption period.
[0087] Specifically, this system determines the facies zone division of drilled wells based on 3D seismic reflection characteristics and the proportion of volcanic sedimentary facies in a single volcanic eruption. Based on the facies zone division results and the characteristic curves corresponding to various rock types within a single volcanic eruption, pseudo-wells of near-source and far-source facies zones are constructed. The characteristic curves of these pseudo-wells are used to determine the range of values corresponding to the proportion of volcanic sedimentary facies in the volcanic structure. Based on this range, facies zones are divided on a planar map representing the facies zones obtained from the 3D seismic data volume interpretation of the volcanic structure, thus identifying the boundaries between near-source and far-source facies zones within the corresponding volcanic eruption period. This system enables the quantitative, accurate, and effective definition of the boundaries between near-source and far-source facies zones, quickly focusing on favorable near-source facies zones within the volcanic structure, facilitating the further search for effective volcanic reservoirs and laying the foundation for exploration targets.
[0088] This invention also provides a machine-readable storage medium storing instructions that, when executed by a processor 100, configure the processor 100 to perform the above-described method for identifying the boundary of a volcanic structure phase zone.
[0089] Machine-readable storage media include both permanent and non-permanent, removable and non-removable media, which can store information by any method or technology. Information can be computer-readable instructions, data structures, modules of programs, or other data. Examples of computer storage media include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, CD-ROM, digital versatile optical disc (DVD) or other optical storage, magnetic tape, magnetic magnetic disk storage or other magnetic storage devices, or any other non-transferable medium that can be used to store information accessible by a computing device. As defined herein, computer-readable media does not include transient computer-readable media, such as modulated data signals and carrier waves.
[0090] The present invention also provides an electronic device 10, which includes a memory 101, a processor 100, and a computer program 102 stored in the memory 101 and executable on the processor 100. When the processor 100 executes the computer program 102, it implements the above-described method for identifying the boundary of a volcanic structure phase zone.
[0091] like Figure 8 The diagram shown is a schematic representation of an electronic device according to an embodiment of the present invention. Figure 8 As shown, the electronic device 10 of this embodiment includes a processor 100, a memory 101, and a computer program 102 stored in the memory 101 and executable on the processor 100. When the processor 100 executes the computer program 102, it implements the steps in the method embodiment described above. Alternatively, when the processor 100 executes the computer program 102, it implements the functions of each module / unit in the device embodiment described above.
[0092] For example, computer program 102 can be divided into one or more modules / units, one or more of which are stored in memory 101 and executed by processor 100 to complete the present invention. One or more modules / units can be a series of computer program instruction segments capable of performing specific functions, which describe the execution process of computer program 102 in electronic device 10. For example, computer program 102 can be divided into a drilled well facies zone delineation module, a pseudo-well construction module, a value range determination module, a planar map acquisition module, and a volcanic structure facies zone delineation module.
[0093] Electronic device 10 can be a desktop computer, laptop, handheld computer, cloud server, or other computing device. Electronic device 10 may include, but is not limited to, processor 100 and memory 101. Those skilled in the art will understand that... Figure 8 This is merely an example of electronic device 10 and does not constitute a limitation on electronic device 10. It may include more or fewer components than shown, or combine certain components, or different components. For example, electronic device may also include input / output devices, network access devices, buses, etc.
[0094] The processor 100 can be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor can be a microprocessor or any conventional processor.
[0095] The memory 101 can be an internal storage unit of the electronic device 10, such as a hard disk or RAM 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, Smart Media Card (SMC), Secure Digital (SD) card, or Flash Card equipped on the electronic device 10. Furthermore, the memory 101 can include both internal and external storage units 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 will clearly understand that, for the sake of convenience and brevity, the above-described division of functional units and modules is merely an example. In practical applications, the above functions can be assigned to different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiments can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit. Furthermore, the specific names of the functional units and modules are only for easy differentiation and are not intended to limit the scope of protection of this application. The specific working process of the units and modules in the above system can be referred to the corresponding process in the foregoing method embodiments, and will not be repeated here.
[0097] Those skilled in the art will understand that embodiments of this application can be provided as a method, system, or computer program 102 product. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program 102 product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.
[0098] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program 102 products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program 102 instructions. These computer program 102 instructions can be provided to a processor 100 of a general-purpose computer, special-purpose computer, embedded processor, 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, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.
[0099] These computer program 102 instructions may also be stored in a computer-readable storage medium 101 that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium 101 produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.
[0100] These computer program 102 instructions may 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, thereby providing instructions that execute on the computer or other programmable apparatus for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.
[0101] It should also be noted that the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article, or apparatus. Unless otherwise specified, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes that element.
[0102] The above are merely embodiments of this application and are not intended to limit the scope of this application. Various modifications and variations can be made to this application by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the scope of the claims of this application.
Claims
1. A method for identifying the boundaries of volcanic structural facies zones, characterized in that, include: Based on the three-dimensional seismic reflection characteristics and the proportion of pre-defined lithofacies types in drilled wells within a single volcanic eruption period, the drilled well facies zones are divided. Based on the results of facies zone division of drilled wells and the characteristic curves corresponding to various rock types in a single volcanic eruption, pseudo-wells of near-source facies zone and pseudo-wells of far-source facies zone are constructed; the characteristic curves include sonic curves and density curves. Based on the characteristic curves of pseudo-wells in the near-source facies zone and pseudo-wells in the far-source facies zone, the range of values corresponding to the pre-defined type of lithofacies proportion of volcanic structures is determined. Based on the 3D seismic data volume of the volcanic structure, a planar map representing the facies zone is obtained; Based on the preset range of lithofacies proportions corresponding to the volcanic structure, facies zones of the volcanic structure are divided on the planar map representing the facies zones, obtaining the boundaries between the proximal and distal facies zones within the corresponding volcanic eruption period, including: Based on the value A corresponding to the R curve, the volcanic structure facies zone is divided on the planar map representing the facies zone, and the locations of the distal and proximal facies zones of the volcanic structure are obtained. This includes: taking the area on the planar map with a value greater than A as the distal facies zone of the volcanic structure, and taking the area with a value less than A as the proximal facies zone of the volcanic structure. Based on the location of the distant source phagocytic zone and the location of the near source phagocytic zone, the boundary between the near source phagocytic zone and the distant source phagocytic zone is obtained; The determination of the value range corresponding to the preset type lithofacies proportion of volcanic structures based on the characteristic curves of near-source facies pseudo-wells and far-source facies pseudo-wells includes: Based on the acoustic and density curves of pseudo-wells in the near-source facies zone and the far-source facies zone, a two-dimensional seismic forward model was established using the convolution model, ray tracing, and wave equation methods. Based on the two-dimensional seismic forward model, the energy half-life T curve and thickness H curve within the corresponding volcanic eruption period are calculated. Based on the energy half-life T curve and thickness H curve within the corresponding volcanic eruption period, the R curve is calculated using the formula R=(T / H)*100. The value A corresponding to the proportion of volcanic sedimentary phases in the R curve is analyzed, and the range of values for A is obtained.
2. The method for identifying the boundary of a volcanic facies zone according to claim 1, characterized in that, Also includes: Based on the drilling data, the volcanic eruption phases are divided; The lithofacies were classified according to the lithofacies of each volcanic eruption period, and the lithofacies classification results were obtained. The lithofacies classification results include eruptive facies, overflow facies, intrusive facies, volcanic conduit facies, and volcanic sedimentary facies.
3. The method for identifying the boundary of a volcanic facies zone according to claim 2, characterized in that, The characteristic curves include the acoustic wave curve and the density curve; Based on the facies zone division results of drilled wells and the characteristic curves corresponding to various rock types in a single volcanic eruption, pseudo-wells of near-source facies zone and pseudo-wells of far-source facies zone are constructed, including: Based on the sonic and density curves of the drilled wells and the lithofacies classification results, the sonic and density curves of the eruptive facies, the overflow facies, and the volcanic sedimentary facies were determined. Based on the facies zone division results of the drilled wells, pseudo-wells of the near-source facies zone are constructed using the sonic curves and density curves of the eruptive and overflow facies zones of the drilled wells, and pseudo-wells of the distant-source facies zone are constructed using the sonic curves and density curves of the volcanic sedimentary facies zones of the drilled wells.
4. The method for identifying the boundary of a volcanic phagocytic structure according to claim 2, characterized in that, The method of dividing drilled well facies zones based on three-dimensional seismic reflection characteristics and the proportion of pre-defined lithofacies types in drilled wells within a single volcanic eruption period includes: Based on three-dimensional seismic reflection characteristics and the proportion of volcanic sedimentary facies in drilled wells within a single volcanic eruption period, the drilled well facies zones are divided as follows: S = volcanic sedimentary facies ÷ sum of all lithofacies; Where S represents the proportion of volcanic sedimentary facies; S below 50% indicates a near-source facies zone, and S above 50% indicates a far-source facies zone.
5. The method for identifying the boundary of a volcanic facies zone according to claim 1, characterized in that, The three-dimensional seismic data volume based on the volcanic structure yields a planar map characterizing the facies zone, including: Based on the three-dimensional seismic data volume of the volcanic structure, the energy half-decay planar diagram t, the bottom and top surfaces of the volcanic structure for a single volcanic eruption are obtained; Based on the bottom and top surfaces of the volcanic structure of a single eruption, calculate the thickness h planar diagram of a single eruption. According to the formula r = (t / h) * 100, the planar diagram representing the phase band is obtained; where r is the planar diagram representing the parameters corresponding to the phase band.
6. A system for identifying the boundaries of volcanic facies zones, characterized in that, A method for performing the identification of volcanic facies boundaries as described in any one of claims 1 to 5, comprising: The drilled facies zone delineation module is used to delineate drilled facies zones based on three-dimensional seismic reflection characteristics and the proportion of pre-defined lithofacies types in drilled wells within a single volcanic eruption period. The pseudo-well construction module is used to construct near-source facies zone pseudo-wells and far-source facies zone pseudo-wells based on the facies zone division results of drilled wells and the characteristic curves corresponding to various types of rocks in a single volcanic eruption period of drilled wells. The value range determination module is used to determine the value range corresponding to the pre-defined type facies proportion of volcanic structures based on the characteristic curves of near-source facies pseudo-wells and far-source facies pseudo-wells. The planar plot module is used to obtain planar plots representing facies zones from 3D seismic data volumes based on volcanic mechanisms. The volcanic structure facies zone division module is used to divide the volcanic structure facies zones on a planar map representing facies zones according to the value range corresponding to the preset type of rock facies proportion of the volcanic structure, and obtain the boundary between the near-source facies zone and the far-source facies zone within the corresponding volcanic eruption period.
7. A machine-readable storage medium storing instructions thereon, characterized in that, When executed by a processor, the instruction causes the processor to be configured to perform the method for identifying the boundary of a volcanic structure phalanx as described in any one of claims 1 to 5.
8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the method for identifying the boundary of a volcanic phalanx as described in any one of claims 1 to 5.