Mudstone target reservoir prediction method and device, electronic equipment and storage medium

Abstract

This disclosure relates to a method, apparatus, electronic device, and storage medium for predicting target reservoirs in argillaceous limestone. The method includes: determining the range of resistivity and acoustic values ​​of target argillaceous limestone reservoirs within a target work area; constructing a geological model based on seismic data within the target work area and resistivity curves, acoustic curves, and TOC curves from sample wells, respectively, and performing waveform inversion to obtain resistivity inversion data volumes, acoustic inversion data volumes, and TOC inversion data volumes; determining the time domain value of argillaceous limestone thickness based on the resistivity and acoustic inversion data volumes and the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoirs; calculating fusion plane attribute values ​​based on the argillaceous limestone thickness time domain value; extracting the TOC inversion data volume attribute values ​​of the target layer; calculating the planar distribution range of the target argillaceous limestone reservoir based on the fusion plane attribute values ​​and TOC planar distribution attribute values; performing inversion on the target argillaceous limestone layer to achieve semi-quantitative prediction of high-quality argillaceous limestone reservoirs and improve reservoir prediction accuracy.

Description

Technical Field

[0001] This disclosure relates to the field of oil and gas exploration technology, and in particular to a method and apparatus for predicting target reservoirs in argillaceous limestone, electronic equipment and storage medium. Background Technology

[0002] Marine basin reefs and shoals are important oil and gas reservoirs and key exploration areas, with argillaceous limestone traditionally considered source rocks. In recent years, with the deepening of oil and gas exploration, high-yield industrial gas flow shows have been observed in argillaceous limestone formations, indicating significant exploration potential. Criteria for evaluating argillaceous limestone reservoirs include significant limestone thickness, high TOC (Total Organic Carbon), and well-developed fractures.

[0003] In the existing technology, the method for predicting the development area of ​​high-quality argillaceous limestone reservoirs is usually as follows: (1) Calculate the thickness of argillaceous limestone: obtain the logging sensitive parameter curve of the work area; collect sample wells, establish a geological model, and use the obtained seismic data to carry out well-seismic joint inversion to extract the thickness of argillaceous limestone; (2) Calculate the total organic carbon content (TOC): obtain the logging sensitive parameter curve of the work area; collect sample wells, use the linear regression equation to calculate the influence factor of each sensitive parameter curve, and finally use the linear weighting method to calculate the TOC curve, or calculate the TOC curve by the formula (lgRt) / Ac; (3) Calculate the degree of fracture development: obtain seismic data and extract the coherence, curvature or ant body attributes of the seismic data; (4) After calculating and evaluating the three standards of argillaceous limestone reservoirs respectively, the range of the three planar maps is superimposed.

[0004] However, existing technologies have at least the following problems: (1) The calculation of the total organic carbon (TOC) curve of argillaceous limestone based on linear regression of sample wells is greatly affected by the sampling and distribution of sample wells. At the same time, due to the different geological conditions and the number of statistical samples in different regions, the calculation formula varies greatly, making it difficult to depict the TOC distribution range on a plane; the TOC calculated based on the formula (lgRt) / Ac is a TOC prediction method based on continental shale oil, which is not suitable for predicting the TOC of argillaceous limestone; (2) The prediction accuracy is low and the calculation speed is slow when three standard favorable ranges are delineated and superimposed. Therefore, at present, there is no effective method for overall fine prediction of argillaceous limestone reservoir development areas. Summary of the Invention

[0005] To solve the above-mentioned technical problems, or at least partially solve them, embodiments of this disclosure provide a method and apparatus for predicting argillaceous limestone reservoirs, electronic equipment, and storage medium.

[0006] In a first aspect, embodiments of this disclosure provide a method for predicting target reservoirs in argillaceous limestone, comprising:

[0007] Determine the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir within the target work area;

[0008] Geological models were constructed based on seismic data from sample wells in the target work area, along with resistivity curves, acoustic curves, and TOC curves. Waveform inversion was then performed to obtain resistivity inversion data volumes, acoustic inversion data volumes, and TOC inversion data volumes. The resistivity curve and acoustic curve are curves after consistency processing.

[0009] The time domain value of mudstone thickness is determined based on the resistivity inversion data volume, the acoustic inversion data volume, and the range of resistivity and acoustic values ​​of the target mudstone reservoir.

[0010] Based on the time domain value of mudstone thickness, the plane attribute value of the fusion of fracture development and mudstone thickness is calculated according to the relationship between the ant body attribute value and the fracture development threshold value in the seismic data.

[0011] Extract the volumetric attribute values ​​of the TOC inversion data of the target layer to obtain the TOC planar distribution attribute values. Based on the plane attribute values ​​of the fracture development and argillaceous limestone thickness and the TOC planar distribution attribute values, calculate the planar distribution range of the target argillaceous limestone reservoir.

[0012] In one possible implementation, determining the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir within the target work area includes:

[0013] The resistivity values ​​of the resistivity curves and the acoustic values ​​of the acoustic curves of the sample wells in the target work area are intersected by scatter plots to obtain an intersection diagram.

[0014] Based on the distribution of scatter points matching the target argillaceous limestone reservoir in the cross-sectional diagram, the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir is determined.

[0015] In one possible implementation, determining the time domain value of mudstone thickness based on the resistivity inversion data volume, the acoustic inversion data volume, and the range of resistivity and acoustic values ​​of the target mudstone limestone reservoir includes:

[0016] In response to the resistivity and acoustic inversion data of the sample points within the resistivity and acoustic values ​​of the target argillaceous limestone reservoir, the time domain value of the mudstone thickness is calculated based on the time values ​​of the bottom and top interfaces of the target argillaceous limestone layer.

[0017] Since the resistivity and acoustic inversion data of the sample points are not within the range of the resistivity and acoustic values ​​of the target mudstone reservoir, the mudstone thickness time domain value is 0.

[0018] In one possible implementation, the time threshold of mudstone thickness is calculated based on the time values ​​of the bottom and top interfaces of the target mudstone layer using the following expression:

[0019] H = T2 - T1

[0020] Where H is the time domain value of mudstone thickness, T2 is the time value of the bottom interface of the target layer of argillaceous limestone, and T1 is the time value of the top interface of the target layer of argillaceous limestone.

[0021] In one possible implementation, the step of calculating the plane attribute value of fracture development and argillaceous limestone thickness based on the mudstone thickness time domain value and according to the relationship between ant body attribute values ​​and fracture development threshold values ​​in seismic data includes:

[0022] In response to the fact that the ant body attribute value in the seismic data of the sample point is greater than the fracture development threshold value, the plane attribute value of the fusion of fracture development and argillaceous limestone thickness is calculated based on the mudstone thickness time domain value.

[0023] In response to the seismic data of sample points, the ant body attribute value is less than or equal to the crack development threshold value, and the plane attribute value of crack development and argillaceous limestone thickness is the mudstone thickness time domain value.

[0024] In one possible implementation, the calculation of the plane property value of the fusion between fracture development and argillaceous limestone thickness is based on the time domain value of mudstone thickness:

[0025] H k =H+500

[0026] Among them, H k The value represents the plane property of the integration of crack development and mudstone thickness, and H represents the time domain value of mudstone thickness.

[0027] In one possible implementation, calculating the planar distribution range of the target argillaceous limestone reservoir based on the plane attribute value of the fracture development and argillaceous limestone thickness integration and the TOC plane distribution attribute value includes:

[0028] Multiply the fracture development and argillaceous limestone thickness plane attribute value of the sample point with the TOC plane distribution attribute value to obtain the plane attribute value of the target argillaceous limestone reservoir at the sample point;

[0029] Based on the planar attribute values ​​of the target argillaceous limestone reservoir at the sample points, the planar distribution range of the target argillaceous limestone reservoir is calculated.

[0030] Secondly, embodiments of this disclosure provide a target reservoir prediction device for argillaceous limestone, comprising:

[0031] The first determining module is used to determine the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir within the target work area;

[0032] The inversion module is used to construct a geological model based on the seismic data in the target work area and the resistivity curve, sonic curve and TOC curve of the sample well, and to perform waveform indication inversion to obtain resistivity inversion data volume, sonic inversion data volume and TOC inversion data volume, wherein the resistivity curve and sonic curve are curves after consistency processing.

[0033] The second determining module is used to determine the time domain value of mudstone thickness based on the resistivity inversion data volume, the acoustic inversion data volume, and the range of resistivity and acoustic values ​​of the target mudstone reservoir.

[0034] The first calculation module is used to calculate the plane attribute value of the fusion between fracture development and argillaceous limestone thickness based on the time domain value of mudstone thickness and the relationship between the ant body attribute value and the fracture development threshold value in the seismic data.

[0035] The second calculation module is used to extract the volume attribute values ​​of the TOC inversion data of the target layer, obtain the TOC planar distribution attribute values, and calculate the planar distribution range of the target reservoir of argillaceous limestone based on the plane attribute values ​​of the fracture development and the thickness of the argillaceous limestone and the TOC planar distribution attribute values.

[0036] Thirdly, embodiments of this disclosure provide an electronic device, including a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other via the communication bus;

[0037] Memory, used to store computer programs;

[0038] The processor, when executing the program stored in memory, implements the above-mentioned method for predicting argillaceous limestone reservoirs.

[0039] Fourthly, embodiments of this disclosure provide a computer-readable storage medium having a computer program stored thereon, characterized in that the computer program, when executed by a processor, implements the above-described method for predicting argillaceous limestone reservoirs.

[0040] Compared with the prior art, the technical solutions provided in this disclosure have at least some or all of the following advantages:

[0041] The argillaceous limestone reservoir prediction method described in this disclosure determines the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir within the target work area. Based on seismic data within the target work area and the resistivity, acoustic, and TOC curves of sample wells, geological models are constructed, and waveform indication inversion is performed to obtain resistivity inversion data volumes, acoustic inversion data volumes, and TOC inversion data volumes. The resistivity and acoustic curves are curves after consistency processing. Based on the resistivity and acoustic inversion data volumes and the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir... The time domain value of mudstone thickness is determined. Based on the time domain value of mudstone thickness, and according to the relationship between the ant body attribute value and the fracture development threshold value in the seismic data, the plane attribute value of the fusion of fracture development and mudstone thickness is calculated. The volume attribute value of TOC inversion data of the target segment is extracted to obtain the TOC plane distribution attribute value. Based on the plane attribute value of the fusion of fracture development and mudstone thickness and the TOC plane distribution value, the plane distribution range of the target mudstone reservoir is calculated. The target mudstone segment is inverted to achieve semi-quantitative prediction of high-quality mudstone reservoirs and improve the accuracy of reservoir prediction. Attached Figure Description

[0042] The accompanying drawings, which are incorporated in and form a part of this specification, illustrate embodiments consistent with this disclosure and, together with the description, serve to explain the principles of this disclosure.

[0043] To more clearly illustrate the technical solutions in the embodiments of this disclosure or the prior art, the accompanying drawings used in the description of the embodiments or related technologies will be briefly introduced below. Obviously, those skilled in the art can obtain other drawings based on these drawings without creative effort.

[0044] Figure 1 A schematic diagram illustrating the flow of a method for predicting argillaceous limestone reservoirs according to embodiments of the present disclosure is shown.

[0045] Figure 2 A schematic diagram illustrating the intersection of sound waves and lg(RT) according to an embodiment of the present disclosure is shown.

[0046] Figure 3 An inversion profile of argillaceous limestone according to an embodiment of the present disclosure is illustrated schematically;

[0047] Figure 4 A TOC content inversion profile according to an embodiment of the present disclosure is illustrated schematically;

[0048] Figure 5 A schematic diagram illustrating a predicted mudstone thickness according to an embodiment of the present disclosure is shown.

[0049] Figure 6A TOC content prediction planar diagram according to an embodiment of the present disclosure is illustrated schematically;

[0050] Figure 7 A schematic diagram illustrating the prediction of high-quality argillaceous limestone reservoirs according to embodiments of the present disclosure is shown.

[0051] Figure 8 This schematic diagram illustrates a structural block diagram of a target reservoir prediction device for argillaceous limestone according to an embodiment of the present disclosure;

[0052] Figure 9 The diagram illustrates a structural block diagram of an electronic device according to an embodiment of the present disclosure. Detailed Implementation

[0053] To make the objectives, technical solutions, and advantages of the embodiments of this disclosure clearer, the technical solutions of the embodiments of this disclosure will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this disclosure. Based on the embodiments of this disclosure, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this disclosure.

[0054] See Figure 1 The embodiments of this disclosure provide a method for predicting target reservoirs in argillaceous limestone, including:

[0055] S1 determines the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir within the target work area.

[0056] S2. Based on the seismic data in the target work area, a geological model is constructed with the resistivity curve, sonic curve and TOC curve of the sample well, and waveform indicator inversion is performed to obtain resistivity inversion data volume, sonic inversion data volume and TOC inversion data volume. The resistivity curve and sonic curve are curves after consistency processing.

[0057] In some embodiments, the TOC curve is generated using the following expression:

[0058]

[0059] Wherein, TOC is the total organic carbon content, Ac is the acoustic value of the acoustic curve, Rt is the resistivity value of the resistivity curve, and C is the measured TOC coefficient of the region.

[0060] S3. Determine the time domain value of mudstone thickness based on the resistivity inversion data volume, the acoustic inversion data volume, and the range of resistivity and acoustic values ​​of the target mudstone reservoir.

[0061] S4, based on the time domain value of mudstone thickness, calculates the plane attribute value of fracture development and mudstone thickness integration based on the relationship between the ant body attribute value and the fracture development threshold value in the seismic data.

[0062] S5. Extract the volume attribute values ​​of the TOC inversion data of the target layer to obtain the TOC planar distribution attribute values. Based on the plane attribute values ​​of the fracture development and the thickness of the argillaceous limestone and the TOC planar distribution attribute values, calculate the planar distribution range of the target reservoir of argillaceous limestone.

[0063] In this embodiment, step S1, determining the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir within the target work area, includes:

[0064] The resistivity values ​​of the resistivity curves and the acoustic values ​​of the acoustic curves of the sample wells in the target work area are intersected by scatter plots to obtain an intersection diagram.

[0065] Based on the distribution of scatter points matching the target argillaceous limestone reservoir in the cross-sectional diagram, the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir is determined.

[0066] In this embodiment, step S3, determining the time domain value of mudstone thickness based on the resistivity inversion data volume, the acoustic inversion data volume, and the range of resistivity and acoustic values ​​of high-quality mudstone limestone reservoirs, includes:

[0067] In response to the resistivity and acoustic inversion data of the sample points within the resistivity and acoustic values ​​of the target argillaceous limestone reservoir, the time domain value of the mudstone thickness is calculated based on the time values ​​of the bottom and top interfaces of the target argillaceous limestone layer.

[0068] Since the resistivity and acoustic inversion data of the sample points are not within the range of the resistivity and acoustic values ​​of the target mudstone reservoir, the mudstone thickness time domain value is 0.

[0069] In this embodiment, the time threshold of mudstone thickness is calculated based on the time values ​​of the bottom and top interfaces of the target mudstone layer using the following expression:

[0070] H = T2 - T1

[0071] Where H is the time domain value of mudstone thickness, T2 is the time value of the bottom interface of the target layer of argillaceous limestone, and T1 is the time value of the top interface of the target layer of argillaceous limestone.

[0072] In this embodiment, step S4, which involves calculating the plane attribute value of the fusion between fracture development and argillaceous limestone thickness based on the time domain value of mudstone thickness and the relationship between the ant body attribute value and the fracture development threshold value in the seismic data, includes:

[0073] In response to the fact that the ant body attribute value in the seismic data of the sample point is greater than the fracture development threshold value, the plane attribute value of the fusion of fracture development and argillaceous limestone thickness is calculated based on the mudstone thickness time domain value.

[0074] In response to the seismic data of sample points, the ant body attribute value is less than or equal to the crack development threshold value, and the plane attribute value of crack development and argillaceous limestone thickness is the mudstone thickness time domain value.

[0075] In some embodiments, the ant-body attribute value in the seismic data is the ant-body attribute value in the seismic data of the target argillaceous limestone segment. The surface attribute value of the ant-body attribute is extracted to obtain a fracture development distribution map, and the fracture development threshold value is determined based on the fracture development distribution map.

[0076] In this embodiment, the plane attribute value of the fusion between crack development and argillaceous limestone thickness is calculated based on the time domain value of mudstone thickness:

[0077] H k =H+500

[0078] Among them, H k The value represents the plane property of the integration of crack development and mudstone thickness, and H represents the time domain value of mudstone thickness.

[0079] In this embodiment, step S5, calculating the planar distribution range of the target argillaceous limestone reservoir based on the plane attribute value of the fusion of fracture development and argillaceous limestone thickness and the TOC plane distribution attribute value, includes:

[0080] Multiply the fracture development and argillaceous limestone thickness plane attribute value of the sample point with the TOC plane distribution attribute value to obtain the plane attribute value of the target argillaceous limestone reservoir at the sample point;

[0081] Based on the planar attribute values ​​of the target argillaceous limestone reservoir at the sample points, the planar distribution range of the target argillaceous limestone reservoir is calculated.

[0082] In some embodiments, the planar distribution range of the target argillaceous limestone reservoir is calculated based on the planar attribute values ​​of the target reservoir at sample points, including:

[0083] Sample points whose planar attribute values ​​of the target argillaceous limestone reservoir exceed the preset threshold constitute the planar distribution range of the target argillaceous limestone reservoir.

[0084] This disclosure performs waveform indicator inversion on argillaceous limestone strata to predict high-quality argillaceous limestone reservoirs and improves the prediction accuracy of argillaceous limestone reservoir development zones.

[0085] See Figures 2 to 7 The application process of the target reservoir prediction method for argillaceous limestone disclosed herein includes:

[0086] Step 1) Obtain seismic data and well logging curves.

[0087] In this embodiment, both seismic data and well logging curves can be acquired using existing instruments and software.

[0088] Step 2) Perform consistency processing on the well logging curves.

[0089] Well logging data suffers from problems such as a large span of acquisition time, numerous logging instrument models, inconsistent calibration standards, and inconsistent operating methods. To eliminate systematic errors between well logging data measured at different times and with different instruments, consistency processing is generally required. Commonly used consistency processing methods include mean and variance analysis, which can relatively unify well logging data measured at different times and with different instruments.

[0090] Step 3) Perform a scatter plot of the resistivity values ​​from the Rt curve and the acoustic values ​​from the Ac curve. Using this plot, determine the threshold value A for reservoir lgRt and the threshold value B for reservoir Ac. Figure 2 As shown.

[0091] Step 4) Using the consistent resistivity Rt curve and acoustic Ac curve, construct geological models based on the seismic data and the Rt and Ac curves respectively, perform waveform indicator inversion to obtain the Rt inversion data volume and Ac inversion data volume, as shown below. Figure 3 As shown, and execute the statement:

[0092] If lgRt

[0093] H = T2 - T1

[0094] Else H=0

[0095] The time domain value H of the mudstone thickness was obtained, in ms. Where Ac is the acoustic value of the acoustic curve, Rt is the resistivity value of the resistivity curve, T2 is the time value of the bottom interface of the target mudstone layer, and T1 is the time value of the top interface of the target mudstone layer.

[0096] Step 5) Calculate the ant-body attributes of the seismic data for the target argillaceous limestone section, extract its surface attribute values, obtain the fracture development and distribution map, and fuse the ant-body attributes with the thickness H obtained in Step 4). Execute the following statement:

[0097] If Ant>k

[0098] Hk = H + 500

[0099] Else Hk=H

[0100] The plane property value Hk, which combines the fracture development and the thickness of the argillaceous limestone, is obtained, such as Figure 5 ​As shown, Ant is the ant body attribute value of the seismic data, and k is the crack development threshold value;

[0101] Step 6) Through expression Generate the TOC curve. Where TOC is the total organic carbon content, Ac is the acoustic value of the acoustic curve, Rt is the resistivity value of the resistivity curve, and C is the measured TOC coefficient for the region.

[0102] Step 7) Based on the seismic data and TOC curve, construct a geological model, perform waveform indicator inversion, and obtain the TOC inversion data volume, such as... Figure 4 As shown;

[0103] Extract the TOC inversion data volume attribute values ​​of the target layer segment to obtain the TOC plane distribution attribute values, such as... Figure 6 As shown, the TOC planar distribution attribute value is multiplied by the fracture development and argillaceous limestone thickness fusion planar attribute obtained in step 5) to finally obtain the planar distribution range of high-quality argillaceous limestone reservoirs, as shown. Figure 7 As shown.

[0104] See Figure 8 The present disclosure provides an apparatus for predicting target reservoirs in argillaceous limestone, comprising:

[0105] The first determining module 11 is used to determine the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir within the target work area;

[0106] The inversion module 12 is used to construct a geological model based on the seismic data of the sample wells in the target work area and the resistivity curve, acoustic curve and TOC curve respectively, and to perform waveform indication inversion to obtain resistivity inversion data volume, acoustic inversion data volume and TOC inversion data volume, wherein the resistivity curve and acoustic curve are curves after consistency processing.

[0107] The second determining module 13 is used to determine the time domain value of mudstone thickness based on the resistivity inversion data volume, the acoustic inversion data volume, and the range of resistivity and acoustic values ​​of the target mudstone reservoir.

[0108] The first calculation module 14 is used to calculate the plane attribute value of fracture development and argillaceous limestone thickness based on the mudstone thickness time domain value and the relationship between the ant body attribute value and the fracture development threshold value in the seismic data.

[0109] The second calculation module 15 is used to extract the volume attribute values ​​of the TOC inversion data of the target layer, obtain the TOC plane distribution attribute values, and calculate the plane distribution range of the target reservoir of argillaceous limestone based on the plane attribute values ​​of the fracture development and the thickness of the argillaceous limestone and the TOC plane distribution attribute values.

[0110] In some embodiments, the first determining module includes:

[0111] The well logging curve acquisition unit is used to acquire well logging curves;

[0112] Seismic data acquisition unit, used to acquire seismic data;

[0113] A consistency processing unit is used to process the well logging curves for consistency.

[0114] The curve scatter intersection unit is used to calculate the obtained resistivity curve and resistance curve threshold value.

[0115] In some embodiments, the inversion module includes:

[0116] Ant body attribute calculation unit, used to calculate the ant body attribute values ​​for acquiring seismic data;

[0117] The argillaceous limestone TOC content curve calculation unit is used to calculate the argillaceous limestone TOC content curve based on the well logging curves after consistency processing.

[0118] A seismic wavelet estimation unit is used to estimate a seismic wavelet based on the seismic data.

[0119] A geological model construction unit is used to construct a geological model based on the seismic wavelet and reservoir property curves.

[0120] The waveform indication inversion unit is used to perform waveform indication inversion based on the geological model and seismic data constructed from the resistivity Rt curve, acoustic Ac curve and total organic carbon content TOC curve, to obtain the inversion data volume.

[0121] In some embodiments, the second computing module includes:

[0122] The target segment extraction unit is used to extract the target segment from the inverted data volume to obtain the planar distribution pattern of the inverted data volume.

[0123] The data volume extraction unit is used to extract data from the inverted data volume to obtain the mudstone thickness, TOC content, and planar distribution map of the ant bodies.

[0124] The specific implementation process of the functions and roles of each unit in the above device can be found in the implementation process of the corresponding steps in the above method, and will not be repeated here.

[0125] For the device embodiments, since they basically correspond to the method embodiments, the relevant parts can be referred to in the description of the method embodiments. The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of the present invention according to actual needs. Those skilled in the art can understand and implement this without creative effort.

[0126] In the above embodiments, any and more of the first determining module 11, inversion module 12, second determining module 13, first calculation module 14, and second calculation module 15 can be combined into one module, or any one of these modules can be split into multiple modules. Alternatively, at least part of the functionality of one or more of these modules can be combined with at least part of the functionality of other modules and implemented in one module. At least one of the first determining module 11, inversion module 12, second determining module 13, first calculation module 14, and second calculation module 15 can be at least partially implemented as hardware circuitry, such as a field-programmable gate array (FPGA), a programmable logic array (PLA), a system-on-a-chip, a system-on-a-substrate, a system-on-package, an application-specific integrated circuit (ASIC), or implemented in hardware or firmware by any other reasonable means of integrating or packaging the circuitry, or implemented in software, hardware, or firmware, or in any appropriate combination of any of these three implementation methods. Alternatively, at least one of the first determining module 11, the inversion module 12, the second determining module 13, the first calculation module 14, and the second calculation module 15 can be at least partially implemented as a computer program module, which can perform corresponding functions when the computer program module is run.

[0127] See Figure 9 The electronic device provided in the embodiments of this disclosure includes a processor 1110, a communication interface 1120, a memory 1130 and a communication bus 1140, wherein the processor 1110, the communication interface 1120 and the memory 1130 communicate with each other through the communication bus 1140.

[0128] Memory 1130 is used to store computer programs;

[0129] When processor 1110 executes the program stored in memory 1130, it implements the following method for predicting target reservoirs in argillaceous limestone:

[0130] Determine the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir within the target work area;

[0131] Geological models were constructed based on seismic data within the target work area and resistivity curves, acoustic curves, and TOC curves from sample wells. Waveform inversion was then performed to obtain resistivity inversion data volumes, acoustic inversion data volumes, and TOC inversion data volumes. The resistivity curve and acoustic curve are curves after consistency processing.

[0132] The time domain value of mudstone thickness is determined based on the resistivity inversion data volume, the acoustic inversion data volume, and the range of resistivity and acoustic values ​​of the target mudstone reservoir.

[0133] Based on the time domain value of mudstone thickness, the plane attribute value of the fusion of fracture development and mudstone thickness is calculated according to the relationship between the ant body attribute value and the fracture development threshold value in the seismic data.

[0134] Extract the volumetric attribute values ​​of the TOC inversion data of the target layer to obtain the TOC planar distribution attribute values. Based on the plane attribute values ​​of the fracture development and argillaceous limestone thickness and the TOC planar distribution attribute values, calculate the planar distribution range of the target argillaceous limestone reservoir.

[0135] The aforementioned communication bus 1140 can be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus, etc. This communication bus 1140 can be divided into an address bus, a data bus, a control bus, etc. For ease of illustration, only one thick line is used to represent it in the figure, but this does not indicate that there is only one bus or one type of bus.

[0136] The communication interface 1120 is used for communication between the above-mentioned electronic device and other devices.

[0137] The memory 1130 may include random access memory (RAM) or non-volatile memory, such as at least one disk storage device. Optionally, the memory 1130 may also be at least one storage device located remotely from the aforementioned processor 1110.

[0138] The processor 1110 mentioned above can be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it can also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components.

[0139] Embodiments of this disclosure also provide a computer-readable storage medium. The computer-readable storage medium stores a computer program that, when executed by a processor, implements the argillaceous limestone reservoir prediction method as described above.

[0140] The computer-readable storage medium may be included in the device / apparatus described in the above embodiments; or it may exist independently and not assembled into the device / apparatus. The computer-readable storage medium carries one or more programs that, when executed, implement the argillaceous limestone reservoir prediction method according to the embodiments of this disclosure.

[0141] According to embodiments of this disclosure, the computer-readable storage medium can be a non-volatile computer-readable storage medium, such as including, but not limited to: portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof. In this disclosure, the computer-readable storage medium can be any tangible medium that contains or stores a program that can be used by or in conjunction with an instruction execution system, apparatus, or device.

[0142] It should be noted that, in this document, relational terms such as "first" and "second" are used merely to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, 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 a process, method, article, or apparatus. Without further limitations, 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 said element.

[0143] The above description is merely a specific embodiment of this disclosure, enabling those skilled in the art to understand or implement it. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of this disclosure. Therefore, this disclosure is not to be limited to the embodiments shown herein, but is to be accorded the widest scope consistent with the principles and novel features claimed herein.

Claims

1. A method for predicting target reservoirs in argillaceous limestone, characterized in that, The method includes: Determine the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir within the target work area; Geological models were constructed based on seismic data and resistivity, acoustic, and TOC curves of sample wells within the target work area. Waveform inversion was then performed to obtain resistivity inversion data volume, acoustic inversion data volume, and TOC inversion data volume. The resistivity curve and acoustic curve are curves after consistency processing. The time domain value of mudstone thickness is determined based on the resistivity inversion data volume, the acoustic inversion data volume, and the range of resistivity and acoustic values ​​of the target mudstone reservoir. Based on the time domain value of mudstone thickness, the plane attribute value of the fusion of fracture development and mudstone thickness is calculated according to the relationship between the ant body attribute value and the fracture development threshold value in the seismic data. Extract the volumetric attribute values ​​of the TOC inversion data of the target layer to obtain the TOC planar distribution attribute values. Based on the plane attribute values ​​of the fracture development and argillaceous limestone thickness and the TOC planar distribution attribute values, calculate the planar distribution range of the target argillaceous limestone reservoir. The calculation of the plane attribute values ​​for fracture development and argillaceous limestone thickness, based on the time domain value of mudstone thickness and the relationship between ant-body attribute values ​​and fracture development threshold values ​​in seismic data, includes: In response to the fact that the ant body attribute value in the seismic data of the sample point is greater than the fracture development threshold value, the plane attribute value of the fusion of fracture development and argillaceous limestone thickness is calculated based on the mudstone thickness time domain value. In response to the seismic data of sample points, the ant body attribute value is less than or equal to the crack development threshold value, and the plane attribute value of crack development and argillaceous limestone thickness is the mudstone thickness time domain value.

2. The method according to claim 1, characterized in that, The range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir within the target work area includes: The resistivity values ​​of the resistivity curves and the acoustic values ​​of the acoustic curves of the sample wells in the target work area are intersected by scatter plots to obtain an intersection diagram. Based on the distribution of scatter points matching the target argillaceous limestone reservoir in the cross-sectional diagram, the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir is determined.

3. The method according to claim 1, characterized in that, The determination of the mudstone thickness time domain value based on the resistivity inversion data volume, the acoustic inversion data volume, and the range of resistivity and acoustic values ​​of the target mudstone limestone reservoir includes: In response to the resistivity and acoustic inversion data of the sample points within the resistivity and acoustic values ​​of the target argillaceous limestone reservoir, the time domain value of the mudstone thickness is calculated based on the time values ​​of the bottom and top interfaces of the target argillaceous limestone layer. Since the resistivity and acoustic inversion data of the sample points are not within the range of the resistivity and acoustic values ​​of the target mudstone reservoir, the mudstone thickness time domain value is 0.

4. The method according to claim 3, characterized in that, The time threshold of mudstone thickness is calculated based on the time values ​​of the bottom and top interfaces of the target mudstone layer using the following expression: H = T2 - T1 Where H is the time domain value of mudstone thickness, T2 is the time value of the bottom interface of the target layer of argillaceous limestone, and T1 is the time value of the top interface of the target layer of argillaceous limestone.

5. The method according to claim 1, characterized in that, The plane property values ​​of fracture development and argillaceous limestone thickness are calculated based on the time domain value of mudstone thickness. H k =H+500 Among them, H k The value represents the plane property of the integration of crack development and mudstone thickness, and H represents the time domain value of mudstone thickness.

6. The method according to claim 1, characterized in that, The calculation of the planar distribution range of the target argillaceous limestone reservoir based on the plane attribute values ​​of the fracture development and argillaceous limestone thickness, and the TOC plane distribution attribute values, includes: Multiply the fracture development and argillaceous limestone thickness plane attribute value of the sample point with the TOC plane distribution attribute value to obtain the plane attribute value of the target argillaceous limestone reservoir at the sample point; Based on the planar attribute values ​​of the target argillaceous limestone reservoir at the sample points, the planar distribution range of the target argillaceous limestone reservoir is calculated.

7. A device for predicting target reservoirs in argillaceous limestone, characterized in that, include: The first determining module is used to determine the range of resistivity and acoustic values ​​of the target argillaceous limestone reservoir within the target work area; The inversion module is used to construct geological models based on the seismic data of sample wells in the target work area and the resistivity curve, acoustic curve and TOC curve respectively, and to perform waveform indicator inversion to obtain resistivity inversion data volume, acoustic inversion data volume and TOC inversion data volume, wherein the resistivity curve and acoustic curve are curves after consistency processing. The second determining module is used to determine the time domain value of mudstone thickness based on the resistivity inversion data volume, the acoustic inversion data volume, and the range of resistivity and acoustic values ​​of the target mudstone reservoir. The first calculation module is used to calculate the plane attribute value of the fusion between fracture development and argillaceous limestone thickness based on the time domain value of mudstone thickness and the relationship between the ant body attribute value and the fracture development threshold value in the seismic data. The second calculation module is used to extract the volume attribute values ​​of the TOC inversion data of the target layer, obtain the TOC planar distribution attribute values, and calculate the planar distribution range of the target reservoir of argillaceous limestone based on the plane attribute values ​​of the fracture development and the thickness of the argillaceous limestone and the TOC planar distribution attribute values. The calculation of the plane attribute values ​​for fracture development and argillaceous limestone thickness, based on the time domain value of mudstone thickness and the relationship between ant-body attribute values ​​and fracture development threshold values ​​in seismic data, includes: In response to the fact that the ant body attribute value in the seismic data of the sample point is greater than the fracture development threshold value, the plane attribute value of the fusion of fracture development and argillaceous limestone thickness is calculated based on the mudstone thickness time domain value. In response to the seismic data of sample points, the ant body attribute value is less than or equal to the crack development threshold value, and the plane attribute value of crack development and argillaceous limestone thickness is the mudstone thickness time domain value.

8. An electronic device, characterized in that, It includes a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus; Memory, used to store computer programs; A processor, when executing a program stored in memory, implements the method for predicting target reservoirs of argillaceous limestone as described in any one of claims 1-7.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the method for predicting target reservoirs in argillaceous limestone as described in any one of claims 1-7.

Images