A multi-information fusion phase-controlled seismic inversion method and apparatus for carbonate hills and shoals

Abstract

This invention discloses a multi-information fusion method and apparatus for seismic inversion of facies-controlled carbonate hillock bodies. The method includes: determining facies-controlled impedance correction values ​​based on digitized sedimentary facies data of the target layer where the hillock body is developed; establishing facies-controlled impedance channels based on paleogeographic thickness data and facies-controlled impedance correction values; establishing a mid-frequency impedance inversion channel based on seismic reflection coefficient data; obtaining an absolute impedance inversion channel by integrating the facies-controlled impedance channel and the mid-frequency impedance inversion channel; and constructing an absolute impedance inversion data volume for identifying hillock bodies developed in the target layer. This method fully considers the influence of paleogeography and sedimentary facies on the distribution of hillock body reservoirs in carbonate strata, and can effectively improve the quantitative prediction accuracy of the distribution pattern, external contour, and internal structure of carbonate hillock body reservoirs.

Description

Technical Field

[0001] This invention relates to the field of oil and gas geophysical exploration technology, and in particular to a multi-information fusion phase-controlled seismic inversion method and apparatus for carbonate hill and shoal bodies. Background Technology

[0002] Seismic inversion technology is an important method in the current quantitative prediction of reservoirs. It combines seismic data, geological laws and drilling data to solve the spatial structure and physical properties of underground reservoirs, and then predicts the spatial distribution law of reservoirs. It is one of the key technologies to solve the reservoir distribution problem in the field of oil and gas exploration.

[0003] There are many seismic inversion methods available today. Among them, the commonly used inversion methods in the oil exploration field include sparse pulse inversion and geostatistical inversion. The sparse pulse inversion method has a low resolution and mainly reflects seismic phenomena. Geostatistical inversion has a high resolution, but requires a large number of wells that are evenly distributed. The inversion results match the well characteristics, but it lacks seismic constraints. Summary of the Invention

[0004] The aforementioned seismic inversion methods are all based on geophysical theories and statistical methods, without comprehensively considering various geological factors and sedimentary concepts. In particular, in the process of predicting carbonate hill and shoal reservoirs, the understanding of paleogeography and sedimentary facies has a clear guiding significance for the distribution of reservoirs. The above methods cannot effectively integrate relevant information, and the reservoir prediction results are difficult to meet the relevant technical requirements in the field of oil and gas exploration.

[0005] In recent years, with the continuous discovery of carbonate oil and gas reservoirs in China, carbonate reservoir resources have received increasing attention from experts and scholars. Developing multi-information fusion seismic inversion methods for carbonate rocks, especially high-precision multi-information fusion seismic inversion methods for carbonate hills and shoals, has become an urgent technical problem to be solved.

[0006] To enrich the process routes and increase the selection space, this invention provides a multi-information fusion phase-controlled carbonate hill and shoal seismic inversion method and apparatus. By combining digitized sedimentary facies maps and paleogeographic maps, a multi-information fusion phase-controlled model is constructed. Then, by combining the mid-frequency impedance of seismic inversion, the distribution of carbonate hill and shoal reservoirs is comprehensively characterized.

[0007] In a first aspect, embodiments of the present invention provide a multi-information fusion phase-controlled carbonate hill and beach seismic inversion method, comprising:

[0008] Based on the digital data of the sedimentary facies of the target layer for the development of hills and shoals, determine the facies control impedance correction values;

[0009] Based on the paleogeographic thickness data and the phase control impedance correction value data, a phase control impedance channel is established;

[0010] Based on seismic reflection coefficient data, a mid-frequency impedance inversion trace is established through inversion.

[0011] By combining the phased impedance channel and the intermediate frequency impedance inversion channel, the absolute impedance inversion channel is obtained;

[0012] An absolute impedance inversion data volume is constructed to identify hills and shoals developed in the target layer.

[0013] Optionally, determining the facies correction value data based on the digitized sedimentary facies data of the target layer for developing hills and shoals includes:

[0014] Based on the digital data of the sedimentary facies of the target layer for the development of hills and shoals, determine the facies control impedance correction coefficient data;

[0015] The phase control impedance correction value is determined by multiplying the phase control impedance correction coefficient by the set initial impedance value.

[0016] Optionally, determining the facies control impedance correction coefficient data based on the digital data of the sedimentary facies of the target layer for developing hills and shoals includes:

[0017] Based on the digital data of the sedimentary facies of the target layer for the development of mound-shoal bodies, the first difference between the numerical values ​​representing the mound-shoal facies and the numerical values ​​representing the background facies is determined.

[0018] For each data point in the digitized sedimentary facies data, a second difference is determined between the value representing the sedimentary facies at that point and the value representing the background facies. The ratio of the second difference to the first difference is determined as the phase control impedance correction coefficient for that point, thus obtaining the phase control impedance correction coefficient data.

[0019] Optionally, establishing a phase-controlled impedance channel based on paleogeographic thickness data and the phase-controlled impedance correction value data includes:

[0020] Based on paleogeographic thickness data, paleogeographic thickness weight data is obtained. The paleogeographic thickness weight is the ratio of paleogeographic thickness to the maximum paleogeographic thickness.

[0021] Based on the paleogeographic thickness weight data and the phase-controlled impedance correction value data, a phase-controlled impedance channel is established by multiplying the paleogeographic thickness weight and the phase-controlled impedance correction value and then adding it to the set initial impedance channel.

[0022] Optionally, the initial impedance at each time point of the initial impedance path is the average impedance of the target layer.

[0023] Optionally, the step of establishing the intermediate frequency impedance inversion channel through inversion includes:

[0024] The intermediate frequency impedance inversion channel is established using the following formula:

[0025] ln(zx,y (t))=ln(Z0)+2∫R x,y (t)dt

[0026] Where x is the horizontal axis variable, y is the vertical axis variable, and t is the time variable; Z x,y (t) represents the inverted intermediate frequency impedance value at time t, along with the x-axis and y-axis; Z0 is the initial intermediate frequency impedance value; R x,y (t) represents the seismic reflection coefficient at the x-axis, y-axis, and time t.

[0027] Optionally, the initial intermediate frequency impedance value is set to the average seismic amplitude of the target layer.

[0028] Optionally, the synthesis of the phased impedance channel and the intermediate frequency impedance inversion channel to obtain the absolute impedance inversion channel includes:

[0029] The absolute impedance inversion channel is obtained by combining the phased impedance channel and the intermediate frequency impedance inversion channel through weighted summation.

[0030] Secondly, embodiments of the present invention provide a multi-information fusion phase-controlled carbonate hill and beach seismic inversion device, comprising:

[0031] The phase control impedance correction value determination module is used to determine the phase control impedance correction value data based on the digital data of the sedimentary facies of the target layer for developing hills and shoals.

[0032] The phase-controlled impedance channel establishment module is used to establish a phase-controlled impedance channel based on the paleogeographic thickness data and the phase-controlled impedance correction value data;

[0033] The mid-frequency impedance inversion module is used to establish a mid-frequency impedance inversion trace based on seismic reflection coefficient data.

[0034] The absolute impedance inversion module is used to integrate the phased impedance channel and the intermediate frequency impedance inversion channel to obtain the absolute impedance inversion channel.

[0035] An absolute impedance inversion data volume construction module is used to construct an absolute impedance inversion data volume for identifying hills and beaches developed in the target layer.

[0036] Thirdly, embodiments of the present invention provide a computer storage medium storing computer-executable instructions, which, when executed by a processor, implement any of the above-mentioned multi-information fusion phase-controlled carbonate hill and beach seismic inversion methods.

[0037] Fourthly, embodiments of this disclosure provide a server, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement any of the above-described multi-information fusion phase-controlled carbonate hill and shoal seismic inversion methods.

[0038] The beneficial effects of the above-described technical solutions provided in the embodiments of the present invention include at least the following:

[0039] The multi-information fusion seismic inversion method for carbonate hillock bodies provided in this invention constructs initial facies impedance data by comprehensively using sedimentary facies data and paleogeomorphic data, fusing multiple information from sedimentary facies and paleogeomorphism. Mid-frequency impedance data is then calculated using seismic data, and absolute impedance is obtained from the facies impedance and mid-frequency impedance. Finally, the absolute impedance from multi-information fusion is integrated to obtain the facies-controlled hillock body inversion results. This method fully considers the influence of paleogeomorphism and sedimentary facies on the distribution of reservoirs in carbonate hillock bodies, effectively improving the quantitative prediction accuracy of reservoir distribution patterns, external contours, and internal structures in carbonate hillock bodies. It has important guiding significance for promoting the exploration and development of oil and gas in deep ancient carbonate rocks and provides key technical guidance for determining the distribution of target geological bodies.

[0040] Other features and advantages of the invention will be set forth in the description which follows, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures particularly pointed out in the written description, claims, and drawings.

[0041] The technical solution of the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. Attached Figure Description

[0042] The accompanying drawings are provided to further illustrate the invention and form part of the specification. They are used in conjunction with embodiments of the invention to explain the invention and do not constitute a limitation thereof. In the drawings:

[0043] Figure 1 This is a flowchart of the multi-information fusion phase-controlled carbonate hill and beach seismic inversion method in Embodiment 1 of the present invention;

[0044] Figure 2 This is a flowchart illustrating the specific implementation of the multi-information fusion phase-controlled carbonate rock hill and beach seismic inversion method in Embodiment 2 of the present invention.

[0045] Figure 3 This is a schematic diagram of the paleogeography of the target layer in Embodiment 2 of the present invention;

[0046] Figure 4 This is a schematic diagram of the target layer deposition phase in Embodiment 2 of the present invention;

[0047] Figure 5 This is a schematic diagram of the initial impedance channel and phase-controlled impedance channel at point B in Embodiment 2 of the present invention;

[0048] Figure 6This is a schematic diagram of the seismic trace and mid-frequency impedance inversion trace at point A in Embodiment 2 of the present invention;

[0049] Figure 7 This is a schematic diagram of the absolute impedance inversion path at point A in Embodiment 2 of the present invention;

[0050] Figure 8 This is a schematic diagram of the absolute impedance inversion path at point B in Embodiment 2 of the present invention;

[0051] Figure 9 This is a schematic diagram of the inversion profile results of the target layered hill and beach body in Embodiment 2 of the present invention;

[0052] Figure 10 This is a schematic diagram of the structure of the multi-information fusion phase-controlled carbonate rock hill and beach seismic inversion device in an embodiment of the present invention. Detailed Implementation

[0053] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.

[0054] It should be understood that the terminology used in this invention is merely for describing particular embodiments and is not intended to limit the invention. Furthermore, with respect to numerical ranges in this invention, it should be understood that each intermediate value between the upper and lower limits of the range is also specifically disclosed. Every smaller range between any stated value or intermediate value within a stated range, and any other stated value or intermediate value within said range, is also included in this invention. The upper and lower limits of these smaller ranges may be independently included or excluded from the range.

[0055] Unless otherwise stated, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains. While only preferred methods and materials have been described herein, any methods and materials similar or equivalent to those described herein may be used in the implementation or testing of this invention. All references to this specification are incorporated by way of citation to disclose and describe methods and / or materials associated with those references. In the event of any conflict with any incorporated reference, the content of this specification shall prevail.

[0056] Example 1

[0057] Embodiment 1 of this invention provides a multi-information fusion phase-controlled seismic inversion method for carbonate hills and beaches, the process of which is as follows: Figure 1 As shown, it includes the following steps:

[0058] Step S11: Determine the facies correction values ​​of the target layer sedimentary facies based on the digital data of the development of the hill and shoal bodies.

[0059] Specifically, the phase control impedance correction value data can be determined through the following steps:

[0060] (1) Determine the phase control impedance correction coefficient data based on the digital data of the sedimentary facies of the target layer for the development of hills and shoals.

[0061] This can be achieved by: using digitized sedimentary facies data of the target layer for developing mound-shoal bodies, determining a first difference between the values ​​representing the mound-shoal facies and those representing the background facies; for each data point in the digitized sedimentary facies data, determining a second difference between the values ​​representing the sedimentary facies at that point and those representing the background facies; and using the ratio of the second difference to the first difference as the facies control impedance correction coefficient for that point, thus obtaining the facies control impedance correction coefficient data.

[0062]

[0063] Where x is the horizontal axis variable and y is the vertical axis variable; C x,y F represents the phase control impedance correction coefficient at the x-axis and y-axis. 丘滩 Let F be the numerical value representing the hill-shoal facies, be a constant, and further, be a positive integer; 背景 Let be the numerical value representing the background phase, be a constant, and further, be less than F. 丘滩 An integer value, for example, can be set to 0.

[0064] Optionally, since the sedimentary facies developed in the target layer are simplified into two types, background facies and hill-shoal facies, it is also possible that for each data point in the digitized sedimentary facies data, if the sedimentary facies value of that point represents the hill-shoal facies, its facies control impedance correction coefficient is set to a positive integer; if the sedimentary facies value of that point represents the background facies, its facies control impedance correction coefficient is set to 0.

[0065] (2) The phase control impedance correction value data is determined by multiplying the phase control impedance correction coefficient by the set initial impedance value.

[0066] If the phase control impedance correction coefficient is determined using the above formula (1), the phase control impedance correction value data can be obtained through the following formula (2):

[0067]

[0068] Wherein, ΔP x,y P0 represents the phase-controlled impedance correction value at the x-axis and y-axis; P0 is the initial impedance value, which can be the average impedance value of the target layer.

[0069] Step S12: Based on the paleogeographic thickness data and phase control impedance correction value data, establish a phase control impedance channel.

[0070] Specifically, this can include: obtaining paleogeographic thickness weight data based on paleogeographic thickness data, where the paleogeographic thickness weight is the ratio of the paleogeographic thickness to the maximum paleogeographic thickness; and establishing a phased impedance channel by multiplying the paleogeographic thickness weight by the phased impedance correction value and then adding it to a set initial impedance channel.

[0071]

[0072] Where x is the horizontal axis variable, y is the vertical axis variable, and t is the time variable; M x,y (t) represents the phase control impedance values ​​at time t, along with the x-axis and y-axis; H x,y H represents the paleogeographic thickness at the x-axis and y-axis. max M represents the maximum thickness of the paleogeography. 初始 (t) represents the value of the initial impedance channel at time t. All values ​​of the initial impedance channel are the same constant, and the average value of the target layer impedance can be taken.

[0073] Step S13: Based on seismic reflection coefficient data, establish the mid-frequency impedance inversion trace through inversion.

[0074] The intermediate frequency impedance inversion channel is established using the following equation (4):

[0075] ln(Z x,y (t))=ln(Z0)+2∫R x,y (t)dt (4)

[0076] Among them, Z x,y (t) represents the mid-frequency impedance inversion value at time t, along with the x-axis and y-axis; Z0 is the initial mid-frequency impedance value, which can be the average seismic amplitude of the target layer; R x,y (t) represents the seismic reflection coefficient at the x-axis, y-axis, and time t.

[0077] Step S14: Combine the phased impedance channel and the intermediate frequency impedance inversion channel to obtain the absolute impedance inversion channel.

[0078] By combining the phased impedance channel and the intermediate frequency impedance inversion channel, the absolute impedance inversion channel is obtained through weighted summation:

[0079] I x,y (t)=C m *M x,y (t)+C z *(Z x,y (t)) (5)

[0080] Among them, I x,y(t) represents the absolute impedance inversion value at the x-axis, y-axis, and time t; C m C is the phase-controlled impedance weighting coefficient, which is a constant; z is the mid-frequency impedance weighting coefficient, which is a constant.

[0081] Step S15: Construct an absolute impedance inversion data volume to identify hills and shoals developed in the target layer.

[0082] The multi-information fusion seismic inversion method for carbonate hillock bodies provided in this invention constructs initial facies impedance data by comprehensively using sedimentary facies data and paleogeomorphic data, fusing multiple information from sedimentary facies and paleogeomorphism. Mid-frequency impedance data is then calculated using seismic data, and absolute impedance is obtained from the facies impedance and mid-frequency impedance. Finally, the absolute impedance from multi-information fusion is integrated to obtain the facies-controlled hillock body inversion results. This method fully considers the influence of paleogeomorphism and sedimentary facies on the distribution of reservoirs in carbonate hillock bodies, effectively improving the quantitative prediction accuracy of reservoir distribution patterns, external contours, and internal structures in carbonate hillock bodies. It has important guiding significance for promoting the exploration and development of oil and gas in deep ancient carbonate rocks and provides key technical guidance for determining the distribution of target geological bodies.

[0083] Example 2

[0084] Embodiment 2 of this invention provides a specific application of the multi-information fusion phase-controlled carbonate hill and beach seismic inversion method. Taking the seismic inversion of carbonate hill and beach in a certain area as an example, the process is as follows: Figure 2 As shown, it includes the following steps:

[0085] Step S21: Read the maximum paleogeographic thickness value from the paleogeographic thickness data of the target layer.

[0086] H max =MAX(H x,y (6)

[0087] Where x is the horizontal axis variable; y is the vertical axis variable; H x,y The x-coordinate and y-coordinate represent the paleogeographic thickness values; MAX is the maximum value function; H max This represents the maximum paleogeographic thickness value in the paleogeographic thickness data.

[0088] See Figure 3 The image shows a schematic diagram of the paleogeography of the target layer. The maximum paleogeographic thickness read is 295.6m. Specifically, the maximum paleogeographic thickness at point A is 200m, and at point B it is 70m.

[0089] Step S22: Determine the phase control impedance correction value based on the digitized sedimentation phase data.

[0090] See Figure 4The diagram shows the sedimentary facies of the target layer. Based on the digitized sedimentary facies data, the phase control impedance correction value data is determined using equation (2) from Example 1.

[0091] by Figure 3 and Figure 4 Taking points A and B as examples, the phase control impedance correction values ​​are shown in Table 1:

[0092] Table 1. Calculation of Phase Control Impedance Correction Values ​​for Points A and B

[0093]

[0094] Step S23: Based on the paleogeographic thickness data and phase control impedance correction value data, establish a phase control impedance channel.

[0095] The phase-controlled impedance channel is established using equation (3) in Example 1.

[0096] by Figure 3 and Figure 4 Taking points A and B as examples, the calculated phase-controlled impedance values ​​are shown in Table 2:

[0097] Table 2. Calculation table of phase control impedance values ​​at points A and B.

[0098]

[0099] See Figure 5 The diagram shows the initial impedance channel and phase-controlled impedance channel at point B.

[0100] Step S24: Based on seismic reflection coefficient data, establish the mid-frequency impedance inversion trace through inversion.

[0101] The intermediate frequency impedance inversion channel is established using equation (4) from Example 1. See [link to example]. Figure 6 The diagram shown is a schematic of the seismic trace and mid-frequency impedance inversion trace at point A.

[0102] Step S25: Combine the phased impedance channel and the intermediate frequency impedance inversion channel to obtain the absolute impedance inversion channel.

[0103] The absolute impedance inversion channel is established using equation (5) in Example 1.

[0104] With phase-controlled impedance weighting coefficient C m Take 0.9, and set the intermediate frequency impedance weighting coefficient C. z Taking 0.8 as an example, Figure 7 This is a schematic diagram of the absolute impedance inversion trace at point A. Figure 8 This is a schematic diagram of the absolute impedance inversion path at point B.

[0105] Step S26: Construct an absolute impedance inversion data volume to identify hills and shoals developed in the target layer.

[0106] Figure 9 A schematic diagram of the inversion profile results of the target hilly beach body.

[0107] Based on the inventive concept of this invention, embodiments of this invention also provide a multi-information fusion phase-controlled carbonate hill and beach seismic inversion device, the structure of which is as follows: Figure 10 As shown, it includes:

[0108] The phase control impedance correction value determination module 101 is used to determine the phase control impedance correction value data based on the digital data of the sedimentary facies of the target layer for developing hills and shoals.

[0109] The phase-controlled impedance channel establishment module 102 is used to establish a phase-controlled impedance channel based on the paleogeographic thickness data and the phase-controlled impedance correction value data.

[0110] The intermediate frequency impedance inversion module 103 is used to establish an intermediate frequency impedance inversion trace based on seismic reflection coefficient data.

[0111] Absolute impedance inversion module 104 is used to integrate the phased impedance channel and the intermediate frequency impedance inversion channel to obtain the absolute impedance inversion channel;

[0112] Absolute impedance inversion data volume construction module 105 is used to construct absolute impedance inversion data volume for identifying hills and beaches developed in the target layer.

[0113] Regarding the apparatus in the above embodiments, the specific manner in which each module performs its operation has been described in detail in the embodiments related to the method, and will not be elaborated upon here.

[0114] Based on the inventive concept of the present invention, embodiments of the present invention also provide a computer storage medium storing computer-executable instructions, which, when executed by a processor, implement the above-mentioned multi-information fusion phase-controlled carbonate hill and beach seismic inversion method.

[0115] Based on the inventive concept of this invention, this embodiment of the invention also provides a server, including: a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, it implements the above-mentioned multi-information fusion phase-controlled carbonate hill and beach seismic inversion method.

[0116] Unless otherwise specifically stated, terms such as processing, calculation, operation, determination, display, etc., may refer to the actions and / or processes of one or more processing or computing systems or similar devices that represent the manipulation and conversion of data representing physical (e.g., electronic) quantities within the registers or memory of the processing system into other data similarly representing physical quantities within the memory, registers, or other such information storage, transmission, or display devices of the processing system. Information and signals can be represented using any of a variety of different techniques and methods. For example, data, instructions, commands, information, signals, bits, symbols, and chips mentioned throughout the above description can be represented by voltage, current, electromagnetic waves, magnetic fields or particles, light fields or particles, or any combination thereof.

[0117] It should be understood that the specific order or hierarchy of steps in the disclosed process is an example of an exemplary method. Based on design preferences, it should be understood that the specific order or hierarchy of steps in the process may be rearranged without departing from the scope of this disclosure. The appended method claims provide elements of various steps in an exemplary order and are not intended to limit the scope to the specific order or hierarchy described.

[0118] In the detailed description above, various features are combined together in a single embodiment to simplify this disclosure. This approach to disclosure should not be construed as reflecting an intention that embodiments of the claimed subject matter require more features than those stated in each claim. Rather, as reflected in the appended claims, the invention is presented with fewer features than all of the features in a single disclosed embodiment. Therefore, the appended claims are hereby clearly incorporated into the detailed description, wherein each claim stands alone as a preferred embodiment of the invention.

[0119] Those skilled in the art will also understand that the various illustrative logic blocks, modules, circuits, and algorithm steps described in conjunction with the embodiments herein can be implemented as electronic hardware, computer software, or a combination thereof. To clearly illustrate the interchangeability between hardware and software, the various illustrative components, blocks, modules, circuits, and steps described above are generally described in terms of their functionality. Whether such functionality is implemented as hardware or software depends on the specific application and the design constraints imposed on the overall system. Those skilled in the art can implement the described functionality in alternative ways for each specific application; however, such implementation decisions should not be construed as departing from the scope of this disclosure.

[0120] The steps of the methods or algorithms described in conjunction with the embodiments herein can be directly embodied in hardware, software modules executed by a processor, or a combination thereof. The software modules can reside in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, hard disks, removable disks, CD-ROMs, or any other form of storage medium well known in the art. An exemplary storage medium is connected to the processor, enabling the processor to read information from and write information to the storage medium. Of course, the storage medium can also be a component of the processor. The processor and storage medium can reside in an ASIC. The ASIC can reside in a user terminal. Alternatively, the processor and storage medium can exist as discrete components in the user terminal.

[0121] For software implementation, the techniques described in this application can be implemented using modules (e.g., procedures, functions, etc.) that perform the functions described in this application. This software code can be stored in memory units and executed by a processor. The memory units can be implemented within the processor or outside the processor; in the latter case, they are communicatively coupled to the processor via various means, as is well known in the art.

[0122] The foregoing description includes examples of one or more embodiments. It is certainly impossible to describe all possible combinations of components or methods in order to describe the above embodiments, but those skilled in the art will recognize that further combinations and arrangements of the various embodiments are possible. Therefore, the embodiments described herein are intended to cover all such changes, modifications, and variations that fall within the scope of the appended claims. Furthermore, the term “comprising” as used in the specification or claims is interpreted in a manner similar to the term “including,” as it is understood when used as a conjunction in the claims. Additionally, the use of any term “or” in the specification of the claims is intended to mean “non-exclusive or.” The terms “first” and “second” are used for descriptive purposes only and should not be construed as indicating or implying relative importance.

Claims

1. A multi-information fusion phase-controlled seismic inversion method for carbonate hills and shoals, characterized in that, include: Based on the digital data of the sedimentary facies of the target layer for the development of mound-shoal bodies, the first difference between the numerical values ​​representing the mound-shoal facies and the numerical values ​​representing the background facies is determined. For each data point in the digitized sedimentary facies data, a second difference is determined between the value representing the sedimentary facies at that point and the value representing the background facies. The ratio of the second difference to the first difference is determined as the phase control impedance correction coefficient for that point, thus obtaining the phase control impedance correction coefficient data. The phase control impedance correction value data is determined by multiplying the phase control impedance correction coefficient by a set initial impedance value. Based on paleogeographic thickness data, paleogeographic thickness weight data is obtained. The paleogeographic thickness weight is the ratio of paleogeographic thickness to the maximum paleogeographic thickness. Based on the paleogeographic thickness weight data and the phase-controlled impedance correction value data, a phase-controlled impedance channel is established by multiplying the paleogeographic thickness weight and the phase-controlled impedance correction value and then adding it to the set initial impedance channel. Based on seismic reflection coefficient data, a mid-frequency impedance inversion trace is established through inversion. By combining the phased impedance channel and the intermediate frequency impedance inversion channel, the absolute impedance inversion channel is obtained; An absolute impedance inversion data volume is constructed to identify hills and shoals developed in the target layer.

2. The method as described in claim 1, characterized in that, The initial impedance at each time point of the set initial impedance channel is the average impedance of the target layer.

3. The method as described in claim 1, characterized in that, The process of establishing the intermediate frequency impedance inversion channel through inversion includes: The intermediate frequency impedance inversion channel is established using the following formula: ； Where x is the horizontal axis variable and y is the vertical axis variable. t is a time variable; x-axis y-axis and the inverted value of the intermediate frequency impedance at time t; To set the initial intermediate frequency impedance value; x-axis y-axis And the seismic reflection coefficient at time t.

4. The method as described in claim 1, characterized in that, The combined phased impedance channel and intermediate frequency impedance inversion channel yield an absolute impedance inversion channel, including: The absolute impedance inversion channel is obtained by combining the phased impedance channel and the intermediate frequency impedance inversion channel through weighted summation.

5. A multi-information fusion phase-controlled carbonate hill and beach seismic inversion device, characterized in that, include: The phase-controlled impedance correction value determination module is used to determine the first difference between the value representing the hill-shoal facies and the value representing the background facies based on the digital data of the sedimentary facies of the target layer for developing hill-shoal bodies. For each data point in the digitized sedimentary facies data, a second difference is determined between the value representing the sedimentary facies at that point and the value representing the background facies. The ratio of the second difference to the first difference is determined as the phase control impedance correction coefficient for that point, thus obtaining the phase control impedance correction coefficient data. The phase control impedance correction value data is determined by multiplying the phase control impedance correction coefficient by a set initial impedance value. The phased impedance channel establishment module is used to obtain paleo-geomorphic thickness weight data based on paleo-geomorphic thickness data. The paleo-geomorphic thickness weight is the ratio of the paleo-geomorphic thickness to the maximum paleo-geomorphic thickness. Based on the paleo-geomorphic thickness weight data and the phased impedance correction value data, the phased impedance channel is established by multiplying the paleo-geomorphic thickness weight and the phased impedance correction value and then adding it to the set initial impedance channel. The mid-frequency impedance inversion module is used to establish a mid-frequency impedance inversion trace based on seismic reflection coefficient data. The absolute impedance inversion module is used to integrate the phased impedance channel and the intermediate frequency impedance inversion channel to obtain the absolute impedance inversion channel. An absolute impedance inversion data volume construction module is used to construct an absolute impedance inversion data volume for identifying hills and beaches developed in the target layer.

6. A computer storage medium, characterized in that, The computer storage medium stores computer-executable instructions, which, when executed by a processor, implement the multi-information fusion phase-controlled carbonate hill and beach seismic inversion method as described in any one of claims 1 to 4.

7. A server, characterized in that, include: The method includes a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the program, implements the multi-information fusion phase-controlled carbonate hill and shoal seismic inversion method according to any one of claims 1 to 4.

Images