A method and system for reconstructing wave impedance curves to eliminate the influence of argillaceous material in carbonate rocks.

Abstract

This invention discloses a method and system for reconstructing acoustic impedance curves to eliminate the influence of argillaceous content in carbonate rocks. The method involves statistically analyzing the median of the normal distribution of acoustic impedance in pure limestone sections and the average well logging GR value; determining the median of the normal distribution of acoustic impedance within different well logging GR value intervals; determining the correction amount of the acoustic impedance of limestone with different well logging GR values ​​relative to the acoustic impedance of pure limestone based on the median of the normal distribution of acoustic impedance within different well logging GR value intervals; fitting the relationship between the acoustic impedance correction amount and the well logging GR value; and using the fitted relationship to correct the acoustic impedance of the argillaceous limestone section to obtain the reconstructed acoustic impedance curve. This method eliminates reservoir misjudgment caused by the decrease in acoustic impedance due to increased argillaceous content, enabling better description and sculpting of carbonate reservoirs and effectively guiding carbonate oil and gas exploration and development.

Description

Technical Field

[0001] This invention belongs to the field of petroleum exploration and relates to a method and system for reconstructing wave impedance curves to eliminate the influence of carbonate mudstone. Background Technology

[0002] Carbonate reservoirs hold a significant position in global oil and gas distribution, accounting for approximately 50% of the world's total reserves and over 60% of its total production. Oil and gas fields formed by carbonate reservoirs often boast large reserves and high single-well production, readily forming large-scale oil and gas fields. Of the nine high-yield wells worldwide that have ever produced over 10,000 tons per day, eight are located in carbonate reservoirs. Many important oil and gas producing regions worldwide are dominated by carbonate reservoirs; in my country, carbonate reservoirs are also widely distributed. The reservoir spaces in carbonate rocks are generally classified into three types: primary pores, caverns, and fractures. Compared to sandstone reservoirs, carbonate reservoirs exhibit greater complexity and diversity due to their more diverse reservoir space types and secondary variations.

[0003] Carbonate reservoirs are predominantly fracture-vuggy, exhibiting low impedance characteristics. However, as the clay content in carbonate rocks increases, their acoustic impedance decreases to varying degrees, directly impacting the use of acoustic impedance to determine reservoir characteristics and the selection of reservoir impedance threshold values.

[0004] Currently, reservoir prediction for argillaceous carbonate rocks often involves performing impedance inversion on a separate section and selecting an impedance threshold value for that section individually. This method struggles to eliminate impedance differences caused by varying argillaceous content; increasing argillaceous content in carbonate rocks can also lead to decreased impedance, easily resulting in misjudgments of reservoirs. Furthermore, it causes significant difficulties and workload for reservoir prediction in vertically continuous argillaceous carbonate formations and pure carbonate formations. Summary of the Invention

[0005] The purpose of this invention is to overcome the shortcomings of the prior art and provide a method and system for reconstructing wave impedance curves to eliminate the influence of argillaceous content in carbonate rocks. This eliminates reservoir misjudgment caused by the decrease in wave impedance due to increased argillaceous content, thereby better describing and sculpting carbonate rock reservoirs and effectively guiding carbonate rock oil and gas exploration and development.

[0006] To achieve the above objectives, the present invention employs the following technical solution:

[0007] A method for reconstructing wave impedance curves to eliminate the influence of argillaceous material in carbonate rocks, comprising the following steps;

[0008] Step 1: Calculate the median of the normal distribution of wave impedance and the average value of well logging GR for the pure limestone section;

[0009] Step 2: Based on the results of Step 1, determine the median of the normal distribution of acoustic impedance within different logging GR value ranges;

[0010] Step 3: Based on the median of the normal distribution of wave impedance within different logging GR value ranges, determine the correction amount of the wave impedance of limestone with different logging GR values ​​relative to the wave impedance of pure limestone.

[0011] Step 4: Fit the relationship between the wave impedance correction and the logging GR value;

[0012] Step 5: Correct the wave impedance of the marl-bearing section using the fitted relationship to obtain the reconstructed wave impedance curve.

[0013] Preferably, in step one, a normal distribution map of the wave impedance of the pure limestone strata is generated using well logging data, and the median of the normal distribution of wave impedance of the pure limestone in the study area is statistically calculated.

[0014] Preferably, in step two, the step size of the increase in the logging GR value is determined by statistical analysis of the logging data, a normal distribution map of the wave impedance within different logging GR value intervals is generated, and the median of the normal distribution of wave impedance within different GR value intervals is determined.

[0015] Preferably, the step size for increasing the logging GR value is tested according to the principles of sufficient data samples, concentrated data distribution, and obvious normal distribution characteristics. The step size is determined by comparing and analyzing multiple test values.

[0016] Preferably, in step three, the median wave impedance of pure limestone is subtracted from the median wave impedance of different well logging GR intervals to obtain the correction amount of the wave impedance of marlstone relative to the wave impedance of pure limestone in different well logging GR intervals.

[0017] Preferably, in step four, the wave impedance correction value and the logging GR value are subjected to cross-interaction analysis to fit the relationship between the two and establish the correspondence between the logging GR value and the wave impedance correction value.

[0018] Furthermore, before the cross-analysis, the wave impedance correction within the logging GR value range is used as the single-point wave impedance correction of the average value of the logging GR value range.

[0019] A wave impedance curve reconstruction system for eliminating the influence of carbonate argillaceous material includes:

[0020] The statistics module is used to calculate the median of the normal distribution of wave impedance and the average value of well logging GR in pure limestone sections;

[0021] The median determination module is used to determine the median of the normal distribution of acoustic impedance within different logging GR value intervals based on the results of the statistical module.

[0022] The correction amount determination module is used to determine the correction amount of the wave impedance of limestone with different well logging GR values ​​relative to the wave impedance of pure limestone, based on the median of the normal distribution of wave impedance within different well logging GR value ranges.

[0023] The fitting module is used to fit the relationship between the wave impedance correction and the logging GR value;

[0024] The wave impedance curve reconstruction module is used to correct the wave impedance of the argillaceous section using the fitted relationship to obtain the reconstructed wave impedance curve.

[0025] A computer device includes a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, performs the step of reconstructing the wave impedance curve to eliminate the influence of carbonate mudstone as described in any of the preceding claims.

[0026] A computer-readable storage medium storing a computer program that, when executed by a processor, performs the steps of reconstructing the wave impedance curve to eliminate the influence of carbonate mudstone as described in any of the preceding claims.

[0027] Compared with the prior art, the present invention has the following beneficial effects:

[0028] This scheme corrects the acoustic impedance curve by fitting the effect of different clay contents on the reduction of acoustic impedance in carbonate rocks. This makes the value range of the corrected acoustic impedance curve consistent with that of the acoustic impedance curve of carbonate rocks without clay, thus eliminating the influence of different clay contents on acoustic impedance from the source. The corrected acoustic impedance curve can meet the requirements of reservoir acoustic impedance inversion, thereby reducing reservoir misjudgment caused by the reduction of acoustic impedance due to increased clay content. This allows for better description and sculpting of carbonate reservoirs, and effectively guides the exploration and development of carbonate oil and gas. Attached Figure Description

[0029] Figure 1 This is the frequency distribution diagram of the pure limestone wave impedance according to the present invention;

[0030] Figure 2 This is a waveform impedance frequency distribution diagram of marl-bearing rocks in different GR intervals according to the present invention;

[0031] Figure 3 This is a cross-plot of the logging GR value and the wave impedance correction ΔP of the present invention;

[0032] Figure 4 This is a comparison of the acoustic impedance curves of the target formation in the first well of this invention before and after clay removal correction;

[0033] Figure 5 This is a comparison of the acoustic impedance curves of the target formation in the second well of the present invention before and after the removal of clay.

[0034] Figure 6 This is a comparison diagram of the low-frequency inversion model before and after the wave impedance curve correction of the present invention. Detailed Implementation

[0035] To enable those skilled in the art to better understand the present invention, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of the present invention.

[0036] It should be noted that the terms “front,” “back,” “left,” “right,” “up,” and “down” used in the following description refer to the directions shown in the attached diagram, while the terms “inside” and “outside” refer to the directions toward or away from the geometric center of a specific component, respectively.

[0037] It should be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings of this invention are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of the invention described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover a non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0038] Unless otherwise defined, 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. The terminology used herein in the specification of this invention is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. The term "and / or" as used herein includes any and all combinations of one or more of the associated listed items.

[0039] The wave impedance curve reconstruction method for eliminating the influence of carbonate rock argillaceous material as described in this invention includes the following steps:

[0040] Step 1: Calculate the median of the normal distribution of wave impedance and the average value of well logging GR in the pure limestone section using well logging data.

[0041] First, a normal distribution map of the wave impedance of the pure limestone section was created using well logging data, and the median (A0) of the normal distribution of wave impedance of the pure limestone in the study area was calculated.

[0042] Secondly, the average GR (gamma) value (G0) of the pure limestone section is determined by statistical analysis of well logging data. The well logging GR value of marl is generally greater than G0. Therefore, the statistical analysis of the well logging GR value and its wave impedance distribution law of marl starts from G0.

[0043] Step 2: Statistically determine the step size (B) of the increase in GR value using well logging data, create a normal distribution map of wave impedance within different GR value intervals, and determine the median of the normal distribution of wave impedance within different GR value intervals.

[0044] The step size (B) of increasing the GR value is tested according to the principle of "sufficient data samples, concentrated data distribution, and obvious normal distribution characteristics". The step size (B) is determined by comparing and analyzing multiple test values.

[0045] Using well logging data, normal distribution maps of wave impedance were constructed for different GR value intervals (G0~G0+B, G0+B~G0+2B, G0+2B~G0+3B, ..., G0+(n-1)B~G0+nB). This involved statistically determining the median (A1, A2, A3, ..., A...) of the normal distribution of wave impedance within different GR value intervals in the study area. n );

[0046] Step 3: Determine the correction amount of the wave impedance of limestone with different mud content (GR value) relative to the wave impedance of pure limestone.

[0047] Subtract the median acoustic impedance of pure limestone (A) from the median acoustic impedance of different GR intervals n -A0), to obtain the correction amount of the wave impedance of the marl-bearing rock relative to the wave impedance of pure limestone in different GR intervals, that is, the correction amount of the wave impedance of G0~G0+B, G0+B~G0+2B, G0+2B~G0+3B, ..., G0+(n-1)B~G0+nB in different GR value intervals are (A1-A0, A2-A0, A3-A0, ..., A n -A0).

[0048] Step 4: Fit the relationship between the wave impedance correction and the logging GR value.

[0049] The wave impedance correction values ​​(A1-A0, A2-A0, A3-A0, ..., A1-A0) within the GR value range (G0~G0+B, G0+B~G0+2B, G0+2B~G0+3B, ..., G0+(n-1)B~G0+nB) are calculated. n -A0) is used as the single-point wave impedance correction amount for the average value of the GR value interval (G0+1 / 2B, G0+3 / 2B, G0+5 / 2B, ..., G0+(2n-1) / 2B).

[0050] The wave impedance correction values ​​(A1-A0, A2-A0, A3-A0, ..., A) are applied to the wave impedance correction values. n Cross-plot analysis was performed between -A0) and the logging GR values ​​(G0+1 / 2B, G0+3 / 2B, G0+5 / 2B, ..., G0+(2n-1) / 2B) to fit the relationship between the two and establish the correspondence between the logging GR values ​​and the wave impedance correction.

[0051] Step 5: Correct the wave impedance of the marl-bearing section (the section with GR value greater than G0) using the fitted relationship to obtain the reconstructed wave impedance curve.

[0052] Since the increase in mud content in limestone will cause a decrease in wave impedance, the corrected wave impedance is obtained by adding a correction value to the original wave impedance value, thus realizing the reconstruction of the wave impedance curve.

[0053] The implementation process of this invention will be described in further detail below with reference to specific data:

[0054] First, the logging curves in the study area need to be environmentally corrected and normalized to obtain consistent logging data, which is crucial for the accuracy and scientific validity of subsequent research.

[0055] Step 1: Calculate the median of the normal distribution of wave impedance and the average value of well logging GR in the pure limestone section using well logging data.

[0056] Using well logging data, a normal distribution map of the acoustic impedance of the pure limestone section was created to determine the median of the acoustic impedance distribution in the pure limestone. For example... Figure 1 As shown, the median of the normal distribution of impedance of pure limestone is 16234, and the average logging GR of pure limestone can be statistically calculated to be 15GAPI.

[0057] Step 2: Statistically determine the step size (B) of the increase in GR value using well logging data, create a normal distribution map of wave impedance within different GR value intervals, and determine the median of the normal distribution of wave impedance within different GR value intervals.

[0058] Since the average logging GR value for pure limestone is 15 GAPI, the statistical analysis of the logging GR values ​​and impedance distribution patterns of marlstone begins at 15 GAPI. Following the principles of "sufficient data samples, concentrated data distribution, and obvious normal distribution characteristics," several step sizes (3, 5, 7, 9) for increasing GR values ​​were tested. After comparative analysis, a step size of 5 was determined to be the most reasonable. This step size ensures both the number of statistical intervals and the completeness of the statistical sample size and patterns within each interval. Based on this, the statistical intervals were determined to start from 15, namely 15-20, 20-25, 25-30, ..., 110-115.

[0059] By using well logging data, the normal distribution of wave impedance of marl in multiple well logging GR intervals was statistically analyzed, and the median wave impedance distribution of each interval was obtained as 16150, 16041, 15931, 15682, ..., 11972.

[0060] Step 3: Determine the correction amount of the wave impedance of limestone with different mud content (GR value) relative to the wave impedance of pure limestone.

[0061] By subtracting the median wave impedance of pure limestone from the median wave impedance of different GR intervals, the correction amount of the wave impedance of marl-bearing rocks relative to the wave impedance of pure limestone in different GR intervals is obtained. The wave impedance correction amounts for GR values ​​in the ranges of 15-20, 20-25, 25-30, ..., 110-115 are 84, 193, 303, ..., 4262, respectively.

[0062] Step 4: Fit the relationship between the wave impedance correction and the logging GR value.

[0063] First, the wave impedance correction values ​​(84, 193, 303, ..., 4262) for different GR value ranges (15-20, 20-25, 25-30, 30-35, ..., 110-115) are used as wave impedance correction values ​​for the mean values ​​of the GR value ranges (17.5, 22.5, 27.5, 32.5, ..., 112.5). This establishes a one-to-one correspondence between the two sets of values, laying the foundation for fitting the relationship between them.

[0064] Cross-plot analysis was performed on the mean GR value interval and the corresponding acoustic impedance correction value to fit the relationship between the acoustic impedance correction value and the logging GR value, such as... Figure 3 As shown, the fitting result is y = 49.16x - 705.89 (x: logging GR value; y: wave impedance correction). The wave impedance correction and the logging GR value show a good linear positive correlation. The correlation coefficient squared reaches 0.985, indicating that the fitting effect is good.

[0065] Step 5: Correct the wave impedance of the argillaceous section (GR value greater than 15GAPI) using the fitted relationship to obtain the reconstructed wave impedance curve.

[0066] Since an increase in the mud content in limestone leads to a decrease in wave impedance, adding a correction value to the original wave impedance value yields the corrected wave impedance, thus reconstructing the wave impedance curve. Figure 4 and Figure 5 As shown, Figure 4 and Figure 5 These are comparison charts of the acoustic impedance curves of two wells before and after correction. The left side of the two charts shows the GR curve, and the right side shows the acoustic impedance curves before and after correction, laying a solid foundation for reservoir inversion using acoustic impedance in the future.

[0067] The application in the study area yielded excellent results. It unified the wave impedance range between marl-bearing rocks and pure limestone, such as... Figure 6 As shown, the corrected acoustic impedance curves have been applied to the acoustic impedance inversion of the block, ensuring that the acoustic impedance values ​​of marl-bearing rocks and pure limestone are within the same range. This reduces reservoir misjudgment caused by increased clay content and effectively solves the problem of selecting acoustic impedance threshold values. The acoustic impedance inversion using the reconstructed curves provides valuable support for the study of fractured-vuggy carbonate reservoirs and the deployment of development wells in the study area.

[0068] This invention presents a technical method for reconstructing acoustic impedance curves based on statistical correction of well logging data with good correlation. It is applicable to the prediction and description of carbonate reservoirs in similar carbonate exploration and development areas (such as the central and northern Tarim Basin). This method eliminates the influence of increased clay content on acoustic impedance and provides acoustic impedance curves with a unified value range. It provides practical guidance for establishing initial models for carbonate reservoir acoustic impedance inversion, quantitative research on carbonate reservoirs, and the deployment of development wells.

[0069] The following are embodiments of the apparatus of the present invention, which can be used to execute embodiments of the method of the present invention. For details not omitted in the apparatus embodiments, please refer to the embodiments of the method of the present invention.

[0070] In another embodiment of the present invention, a wave impedance curve reconstruction system for eliminating the influence of carbonate rock argillaceous material is provided. This wave impedance curve reconstruction system for eliminating the influence of carbonate rock argillaceous material can be used to implement the above-mentioned wave impedance curve reconstruction method for eliminating the influence of carbonate rock argillaceous material. Specifically, the wave impedance curve reconstruction system for eliminating the influence of carbonate rock argillaceous material includes a statistical module, a median determination module, a correction amount determination module, a fitting module, and a wave impedance curve reconstruction module.

[0071] The statistics module is used to calculate the median of the normal distribution of wave impedance and the average value of well logging GR in the pure limestone section.

[0072] The median determination module is used to determine the median of the normal distribution of wave impedance within different logging GR value intervals based on the results of the statistical module.

[0073] The correction amount determination module is used to determine the correction amount of the wave impedance of limestone with different well logging GR values ​​relative to the wave impedance of pure limestone, based on the median of the normal distribution of wave impedance within different well logging GR value ranges.

[0074] The fitting module is used to fit the relationship between the wave impedance correction and the logging GR value.

[0075] The wave impedance curve reconstruction module is used to correct the wave impedance of the argillaceous section using the fitted relationship to obtain the reconstructed wave impedance curve.

[0076] In another embodiment of the present invention, a terminal device is provided, comprising a processor and a memory. The memory stores a computer program, the computer program including program instructions, and the processor executes the program instructions stored in the computer storage medium. The processor may be a Central Processing Unit (CPU), or it may be other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), or field-programmable gate arrays (FPGAs). Gate Array (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc., are the computing and control core of the terminal. They are suitable for implementing one or more instructions, specifically suitable for loading and executing one or more instructions to realize the corresponding method flow or corresponding function. The processor described in this embodiment of the invention can be used to operate the wave impedance curve reconstruction method to eliminate the influence of argillaceous carbonate rocks, including step one, statistically analyzing the median of the normal distribution of wave impedance in pure limestone sections and the average value of well logging GR; step two, based on the results of step one, determining the median of the normal distribution of wave impedance within different well logging GR value intervals; step three, based on the median of the normal distribution of wave impedance within different well logging GR value intervals, determining the correction amount of the wave impedance of limestone with different well logging GR values ​​relative to the wave impedance of pure limestone; step four, fitting the relationship between the wave impedance correction amount and the well logging GR value; step five, using the fitted relationship to correct the wave impedance of the argillaceous limestone section to obtain the reconstructed wave impedance curve.

[0077] In another embodiment, the present invention also provides a computer-readable storage medium (Memory), which is a memory device in a terminal device for storing programs and data. It is understood that the computer-readable storage medium here may include both the built-in storage medium in the terminal device and extended storage media supported by the terminal device. The computer-readable storage medium provides storage space that stores the terminal's operating system. Furthermore, the storage space also stores one or more instructions suitable for loading and execution by a processor, which may be one or more computer programs (including program code). It should be noted that the computer-readable storage medium here may be high-speed RAM or non-volatile memory, such as at least one disk storage device.

[0078] One or more instructions stored in a computer-readable storage medium can be loaded and executed by a processor to implement the corresponding steps of the wave impedance curve reconstruction method for eliminating the influence of argillaceous carbonate rocks in the above embodiments. One or more instructions in the computer-readable storage medium are loaded by the processor and executed as follows: Step 1, statistically analyzing the median of the normal distribution of wave impedance and the average value of well logging GR in the pure limestone section; Step 2, determining the median of the normal distribution of wave impedance within different well logging GR value intervals based on the results of Step 1; Step 3, determining the correction amount of the wave impedance of limestone with different well logging GR values ​​relative to the wave impedance of pure limestone based on the median of the normal distribution of wave impedance within different well logging GR value intervals; Step 4, fitting the relationship between the wave impedance correction amount and the well logging GR value; Step 5, using the fitted relationship to correct the wave impedance of the argillaceous limestone section to obtain the reconstructed wave impedance curve.

[0079] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. 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 product embodied 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.

[0080] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program 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 instructions. These computer program instructions can be provided to a processor 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 of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0081] These computer program instructions may also be stored in a computer-readable storage medium 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 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.

[0082] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment 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.

[0083] It should be understood that the above description is for illustrative purposes and not for limitation. Many embodiments and applications beyond the provided examples will be apparent to those skilled in the art upon reading the above description. Therefore, the scope of this teaching should not be determined by reference to the above description, but rather by reference to the foregoing claims and the full scope of their equivalents. For purposes of completeness, all articles and references, including patent applications and publications, are incorporated herein by reference. The omission of any aspect of the subject matter disclosed herein in the foregoing claims is not intended as a waiver of that subject matter, nor should it be construed as an indication that the applicant has not considered that subject matter as part of the disclosed inventive subject matter.

Claims

1. A method for reconstructing wave impedance curves to eliminate the influence of argillaceous material in carbonate rocks, characterized in that, Includes the following steps; Step 1: Calculate the median of the normal distribution of wave impedance and the average value of well logging GR for the pure limestone section; Step 2: Based on the results of Step 1, determine the median of the normal distribution of acoustic impedance within different logging GR value ranges; Step 3: Based on the median of the normal distribution of wave impedance within different logging GR value ranges, determine the correction amount of the wave impedance of limestone with different logging GR values ​​relative to the wave impedance of pure limestone. Step 4: Fit the relationship between the wave impedance correction and the logging GR value; Step 5: Correct the wave impedance of the marl-bearing section using the fitted relationship to obtain the reconstructed wave impedance curve.

2. The method for reconstructing wave impedance curves to eliminate the influence of carbonate rock argillaceous material according to claim 1, characterized in that, In step one, a normal distribution map of the wave impedance of the pure limestone strata is generated using well logging data, and the median of the normal distribution of wave impedance of the pure limestone in the study area is calculated.

3. The method for reconstructing wave impedance curves to eliminate the influence of carbonate rock argillaceous material according to claim 1, characterized in that, In step two, the step size of the increase in logging GR value is determined by statistical analysis of logging data, and a normal distribution map of acoustic impedance within different logging GR value intervals is created to determine the median of the normal distribution of acoustic impedance within different GR value intervals.

4. The method for reconstructing wave impedance curves to eliminate the influence of carbonate rock argillaceous material according to claim 1, characterized in that, The step size for increasing the logging GR value is determined by testing the data sample size, the data distribution size, and the normal distribution characteristics, based on the principles of sufficient data sample size, concentrated data distribution, and obvious normal distribution characteristics. The step size is determined by comparing and analyzing multiple test values.

5. The method for reconstructing wave impedance curves to eliminate the influence of carbonate rock argillaceous material according to claim 1, characterized in that, In step three, the median wave impedance of pure limestone is subtracted from the median wave impedance of GR intervals in different well logging to obtain the correction amount of the wave impedance of marlstone relative to the wave impedance of pure limestone in different well logging GR intervals.

6. The method for reconstructing wave impedance curves to eliminate the influence of carbonate rock argillaceous material according to claim 1, characterized in that, In step four, the wave impedance correction value and the logging GR value are cross-analyzed to fit the relationship between the two and establish the correspondence between the logging GR value and the wave impedance correction value.

7. The method for reconstructing wave impedance curves to eliminate the influence of carbonate rock argillaceous material according to claim 6, characterized in that, Before cross-analysis, the wave impedance correction within the well logging GR value range is used as the single-point wave impedance correction of the well logging GR value range mean.

8. A wave impedance curve reconstruction system for eliminating the influence of argillaceous material in carbonate rocks, characterized in that, include: The statistics module is used to calculate the median of the normal distribution of wave impedance and the average value of well logging GR in pure limestone sections; The median determination module is used to determine the median of the normal distribution of acoustic impedance within different logging GR value intervals based on the results of the statistical module. The correction amount determination module is used to determine the correction amount of the wave impedance of limestone with different well logging GR values ​​relative to the wave impedance of pure limestone, based on the median of the normal distribution of wave impedance within different well logging GR value ranges. The fitting module is used to fit the relationship between the wave impedance correction and the logging GR value; The wave impedance curve reconstruction module is used to correct the wave impedance of the argillaceous section using the fitted relationship to obtain the reconstructed wave impedance curve.

9. A computer 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 step of reconstructing the wave impedance curve to eliminate the influence of carbonate mudstone as described in any one of claims 1 to 7.

10. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the step of reconstructing the wave impedance curve to eliminate the influence of carbonate mudstone as described in any one of claims 1 to 7.

Images