Microtopographic delineation method and apparatus
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2020-11-09
- Publication Date
- 2026-03-03
Smart Images

Figure CN114460641B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of petroleum exploration, and more particularly to a method and apparatus for micro-topography characterization. Background Technology
[0002] The Tarim Basin is rich in Paleozoic carbonate oil and gas resources, among which ultra-deep fractured-vuggy carbonate oil and gas reservoirs are a type of reservoir. These reservoirs are characterized by deep burial, strong reservoir heterogeneity, and special reflection characteristics. Their burial depth exceeds 5500m, the reservoir connectivity is complex, the reflection characteristics are beaded, and the local oil-water relationship is complex. In order to characterize the distribution characteristics of remaining oil in carbonate oil and gas reservoirs, the existing method is to use structural characterization technology.
[0003] Existing structural characterization techniques directly interpret the structure at the top of the target layer and then correct and map it based on the actual drilling data. The resulting structural map can reflect the overall structural elevation of a block, but it is difficult to reflect the structural elevation of individual reservoirs. As a result, when characterizing the distribution characteristics of ultra-deep fractured-vuggy carbonate oil and gas reservoirs, it is difficult to accurately describe the local structure and thus cannot provide a favorable basis for the selection of development well locations and the sidetracking of failed or inefficient wells.
[0004] Therefore, how to structurally characterize ultra-deep fractured-vuggy carbonate oil and gas reservoirs is an urgent problem to be solved. Summary of the Invention
[0005] This invention provides a micro-topography characterization method to achieve a detailed description of reservoir distribution characteristics, thereby providing a favorable basis for the optimal selection of development well locations and the sidetracking of failed or inefficient wells.
[0006] In a first aspect, the present invention provides a method for micro-topography characterization, comprising:
[0007] Based on the geological information and seismic data of the drilled wells, the drilled wells are calibrated to obtain calibration results;
[0008] Based on the calibration results, the structural interpretation of the oil and gas reservoirs is determined;
[0009] The trend surface analysis method is applied to process the structural interpretation and determine the structural trend;
[0010] Based on the structural interpretation and structural trend, structural anomalies are identified, and these anomalies are used to represent micro-topography.
[0011] In one possible design, based on the geological information and seismic data of the drilled wells, the drilled wells are calibrated to obtain calibration results, including:
[0012] Based on geological information, the drilled wells were sorted into layers to determine the results of the layer sorting.
[0013] Based on seismic data, spectral analysis was performed on the drilled wells to determine the spectral analysis results;
[0014] Based on the results of the stratification analysis and the spectral analysis, the drilled wells were calibrated, and the calibration results were obtained.
[0015] In one possible design, the structural interpretation of the oil and gas reservoir is determined based on the calibration results, including:
[0016] Based on the location of the oil and gas reservoirs, determine the distance between each marker layer and the oil and gas reservoir;
[0017] Based on the calibration results, the marker layer with the smallest distance was selected for structural interpretation.
[0018] In one possible design, trend surface analysis is applied to process the structural interpretation and determine the structural trend, including:
[0019] Define the data scatter values of the flag layer as (x i ,y i ,z i ), where x i Let y be the x-coordinate of the i-th data point. i Let z be the ordinate of the i-th data point. i Let z be the observation value of the i-th data point; i The corresponding trend value is
[0020] The first-order surface equation for the trend value is:
[0021] According to the least squares principle, the parameters b0, b1, and b2 are obtained by differentiation; where b0 is a number greater than 0, b1 is a number greater than 0, and b2 is a number greater than 0.
[0022] Substitute the parameters into the equation of a linear surface to determine the trend value.
[0023] In one possible design, structural anomalies are identified based on structural interpretation and trends. These anomalies are used to represent micro-topography, including:
[0024] According to the formula Identify and construct outliers;
[0025] Based on the geological information from the drilled wells and structural anomalies, the micro-topography is determined.
[0026] Secondly, the present invention also provides a log detection device, comprising:
[0027] The first determination module is used to calibrate the drilled wells based on the geological information and seismic data, obtain calibration results, and determine the structural interpretation of the oil and gas reservoirs based on the calibration results.
[0028] The processing module is used to apply trend surface analysis to process the construction interpretation and determine the construction trend;
[0029] The second determination module is used to determine structural anomalies based on structural interpretation and structural trends. Structural anomalies are used to represent micro-topography.
[0030] In one possible design, the first determined module is specifically used for:
[0031] Based on geological information, the drilled wells were sorted into layers to determine the results of the layer sorting.
[0032] Based on seismic data, spectral analysis was performed on the drilled wells to determine the spectral analysis results;
[0033] Based on the results of the stratification analysis and the spectral analysis, the drilled wells were calibrated, and the calibration results were obtained.
[0034] In one possible design, the first determined module is specifically used for:
[0035] Based on the location of the oil and gas reservoirs, determine the distance between each marker layer and the oil and gas reservoir;
[0036] Based on the calibration results, the marker layer with the smallest distance was selected for structural interpretation.
[0037] In one possible design, the processing module is specifically used for:
[0038] Define the data scatter values of the flag layer as (x i ,y i ,z i ), where x i Let y be the x-coordinate of the i-th data point. i Let z be the ordinate of the i-th data point. i Let z be the observation value of the i-th data point; i The corresponding trend value is
[0039] The first-order surface equation for the trend value is:
[0040] According to the least squares principle, the parameters b0, b1, and b2 are obtained by differentiation; where b0 is a number greater than 0, b1 is a number greater than 0, and b2 is a number greater than 0.
[0041] Substitute the parameters into the equation of a linear surface to determine the trend value.
[0042] In one possible design, the second determining module is specifically used for:
[0043] According to the formula Identify and construct outliers;
[0044] Based on the geological information from the drilled wells and structural anomalies, the micro-topography is determined.
[0045] Thirdly, the present invention also provides a construction interpretation software platform, comprising:
[0046] Processor; and,
[0047] Memory is used to store the processor's executable instructions;
[0048] The processor is configured to execute any of the micro-topography methods in the first aspect by executing executable instructions.
[0049] Fourthly, embodiments of the present invention also provide a storage medium storing a computer program thereon, which, when executed by a processor, implements any of the micro-topography characterization methods in the first aspect.
[0050] This invention provides a method and apparatus for micro-topography characterization. The method involves calibrating drilled wells based on geological information and seismic data to obtain calibration results; determining the structural interpretation of oil and gas reservoirs based on the calibration results; processing the structural interpretation using trend surface analysis to determine structural trends; and identifying structural anomalies based on the structural interpretation and trends. These anomalies are used to represent micro-topography, enabling a detailed description of reservoir distribution characteristics and providing a favorable basis for optimal well location selection and sidetracking of unsuccessful or inefficient wells. Attached Figure Description
[0051] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0052] Figure 1 This is an application scenario diagram illustrating the micro-topography characterization method of the present invention according to an example embodiment;
[0053] Figure 2 This is a schematic flowchart illustrating a micro-topography characterization method according to an exemplary embodiment of the present invention;
[0054] Figure 3 This is a schematic diagram of the overall process of the micro-topography characterization method according to an example embodiment of the present invention;
[0055] Figure 4 This is a schematic diagram of the micro-topography characterization device according to an exemplary embodiment of the present invention;
[0056] Figure 5 This is a schematic diagram illustrating the structure of the construction explanation software platform according to an example embodiment of the present invention. Detailed Implementation
[0057] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.
[0058] The terms “first,” “second,” “third,” “fourth,” etc. (if present) 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 embodiments of the invention described herein can be implemented, for example, 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.
[0059] The technical solutions of the present invention and how they solve the above-mentioned technical problems will be described in detail below with reference to specific embodiments. The embodiments of the present invention will be described below in conjunction with the accompanying drawings.
[0060] Figure 1 This is an application scenario diagram illustrating the micro-topography characterization method of the present invention according to an example embodiment, such as... Figure 1As shown, the geological anomaly 101, micro-topography 102, and micro-topography 103 are included. Geological anomaly 101 is the envelope surface derived by performing spectral analysis on seismic data from the location of drilled wells and analyzing the subsurface strata of the drilled wells. Because oil and gas reservoirs differentiate according to gravity, their density is relatively low, and they tend to migrate towards structural highs. Based on this characteristic of oil and gas migration, the marker layer where geological anomaly 101 is located is structurally interpreted. Furthermore, a smoothing method is selected to smooth the structural interpretation results, eliminating the influence of regional variation factors and obtaining the overall macroscopic structural trend. Finally, the result of the refined interpretation is subtracted from the smoothed trend surface result to highlight local structural anomalies, obtaining the "structural remnants" after the macroscopic structural trend, namely micro-topography 102 and micro-topography 103. This processing method achieves a refined description of the distribution characteristics of oil and gas reservoirs, thus providing a favorable basis for the selection of development well locations and the sidetracking of unsuccessful or inefficient wells.
[0061] Figure 2 This is a schematic flowchart illustrating a micro-topography characterization method according to an exemplary embodiment of the present invention; as shown below. Figure 2 As shown, the micro-topography characterization method provided in this embodiment includes:
[0062] Step 201: Based on the geological information and seismic data of the drilled wells, calibrate the drilled wells and obtain the calibration results;
[0063] Specifically, based on geological information, the drilled wells are classified into layers to determine the classification results; based on seismic data, the drilled wells are subjected to spectral analysis to determine the spectral analysis results; based on the classification results and spectral analysis results, the drilled wells are calibrated to obtain the calibration results.
[0064] Seismic data refers to artificially generated seismic waves, such as those generated by explosives. These waves propagate downwards, are reflected, and then received at the surface from underground geological bodies. Therefore, a seismic wave is a superposition of reflected waves from all underground geological bodies. Since the reflection characteristics of fractured-vuggy carbonate oil and gas reservoirs resemble a beaded pattern, software algorithms are used to perform spectral analysis on the development pattern of this beaded pattern to extract a specific frequency band of seismic waves sensitive to fractured-vuggy carbonate reservoirs. Analysis of this seismic wave reveals its propagation velocity and the physical structure of the strata in which it originates.
[0065] The aforementioned seismic data analyzes seismic waves from drilled well locations. Next, the sonic transit time data and density logging data from drilled wells are analyzed. The sonic transit time data is acquired at the drilling velocity, and the density logging data is acquired at the time of drilling, representing the physical structure of the formation. These two data points are compared with their corresponding counterparts in the seismic data. When discrepancies are found, geological anomalies are identified. Considering that in actual production, well logging data may not cover the entire well section from top to bottom, leading to incomplete correspondence with seismic data during stratigraphic calibration, a small tolerance range is pre-set in practice. Based on this tolerance, the formation is stratified. During layer-by-layer calibration, the aforementioned method is used to identify laterally continuous stratigraphic layers that show discrepancies between the seismic and drilling data. The resulting calibration results for multiple stratigraphic layers serve as a preliminary screening for geological anomalies.
[0066] To more clearly determine the micro-geomorphology of oil and gas reservoirs, it is necessary to further refine the geological anomalies selected in the preliminary screening to narrow down the scope for subsequent processing, i.e., proceed to step 202.
[0067] Step 202: Based on the calibration results, determine the structural interpretation of the oil and gas reservoir;
[0068] Specifically, based on the location of the oil and gas reservoir, the distance between each marker layer and the oil and gas reservoir is determined; based on the calibration results of multiple layers obtained above, the marker layer with the smallest distance is selected as the geological anomaly for structural interpretation.
[0069] This method of finely calibrating drilled wells and then using the calibration results to select the marker layers closest to the reservoir and distributed in the region as geological anomalies for structural interpretation can ensure that the true structural location of oil and gas reservoirs is reflected to the greatest extent.
[0070] After identifying a geological anomaly, further processing is required because the anomaly contains multiple fissures and cavities, which may be rich in oil and gas or contain other impurities such as soil.
[0071] Step 203: Apply trend surface analysis to process the structural interpretation and determine the structural trend;
[0072] Specifically, the data scatter values of the flag layer are defined as (x i ,y i ,z i ), where x i Let y be the x-coordinate of the i-th data point. i Let z be the ordinate of the i-th data point. i Let z be the observation value of the i-th data point, i.e., the absolute depth of the i-th data point; observation value z iThe corresponding trend value is That is, the relative depth of the i-th data point. These scattered data points constitute the construction explanation of the marker layer.
[0073] The tilt of the marker layer is determined based on these scatter points. The tilt is divided into three types: gentle slope and sloping areas, areas with relatively small structural undulations, and areas with relatively large structural undulations.
[0074] When applying trend surface analysis, if the indicator layer has a gentle slope or relatively small structural undulations,
[0075] The first-order surface equation for the trend value is:
[0076]
[0077] To find the undetermined coefficient b i According to the least squares principle
[0078]
[0079] The derivative equation is as follows:
[0080]
[0081] Solving for the parameters b0, b1, and b2 yields the values b0, b1, and b2, where b0 is a number greater than 0, b1 is a number greater than 0, and b2 is a number greater than 0.
[0082] Substitute these parameters into the equation of a first-order surface to determine the trend value.
[0083] When the inclination of the marker layer results in relatively large undulations, the quadratic surface equation is constructed simultaneously with the trend value of the first-order surface equation:
[0084]
[0085] The parameters are obtained by differentiating according to the least squares principle, and these parameters are substituted into the equation of the quadratic surface to determine the trend value.
[0086] Here, the relationship between the number of terms in the equation and the degree of the equation is as follows:
[0087]
[0088] Where q represents the number of terms in the equation; p represents the degree of the equation; the number of terms increases with the degree.
[0089] The above processing method uses trend surface analysis to construct different surface equations for different marker layers, and selects different smoothing methods to smooth the construction interpretation results, thereby eliminating the influence of regional variation factors and obtaining the overall macro trend of the construction.
[0090] Step 204: Based on the structural interpretation and structural trend, determine the structural anomalies. The structural anomalies are used to represent micro-topography.
[0091] Specifically, according to the formula Identify structural anomalies; based on the geological information from drilled wells and structural anomalies, determine the micro-topography.
[0092] That is, after determining the trend value, the difference between the observed value and the trend value is calculated to obtain the residual value:
[0093]
[0094] Based on the residual value Δz i Highlighting local structural anomalies, a planar distribution map of these anomalies is drawn and compared with existing drilling data for verification, thereby accurately depicting the micro-topography of oil and gas reservoirs.
[0095] Figure 3 This is a schematic diagram illustrating the overall process of the micro-topography characterization method according to an exemplary embodiment of the present invention; as shown below. Figure 3 As shown, the overall process combines the processing methods of steps 201-204 above. On the one hand, based on the trend values of the data scatter points obtained by solving, a certain interval is preset to draw a contour map of regional changes, i.e., a trend map, so as to study the regional change patterns. On the other hand, by calculating the residual values, local anomalies are discovered, local change patterns are studied, and a detailed micro-topographical characterization is obtained.
[0096] More specifically, based on a detailed micro-topographical characterization, combined with reservoir scale and oil source fracture distribution, favorable oil and gas accumulation areas can be effectively identified. Based on the characteristics of these identified favorable accumulation areas, analogy methods can be applied to guide the selection of well locations in undrilled areas or the sidetracking of unsuccessful or inefficient well sites. In other words, this micro-topographical characterization method precisely characterizes the local structural relationships of oil and gas reservoirs, thereby providing a direct basis for the development of carbonate reservoirs.
[0097] Figure 4 This is a schematic diagram of the micro-topography characterization device according to an exemplary embodiment of the present invention; as shown below. Figure 4 As shown, the micro-topography characterization device 40 provided in this embodiment includes:
[0098] The first determining module 401 is used to calibrate the drilled wells based on the geological information and seismic data of the drilled wells, obtain calibration results, and determine the structural interpretation of the oil and gas reservoirs based on the calibration results.
[0099] Processing module 402 is used to apply trend surface analysis to process the construction interpretation and determine the construction trend;
[0100] The second determining module 403 is used to determine structural anomalies based on structural interpretation and structural trends. Structural anomalies are used to represent micro-topography.
[0101] In one possible design, the first determining module 401 is specifically used for:
[0102] Based on geological information, the drilled wells were sorted into layers to determine the results of the layer sorting.
[0103] Based on seismic data, spectral analysis was performed on the drilled wells to determine the spectral analysis results;
[0104] Based on the results of the stratification analysis and the spectral analysis, the drilled wells were calibrated, and the calibration results were obtained.
[0105] In one possible design, the first determining module 401 is specifically used for:
[0106] Based on the location of the oil and gas reservoirs, determine the distance between each marker layer and the oil and gas reservoir;
[0107] Based on the calibration results, the marker layer with the smallest distance was selected for structural interpretation.
[0108] In one possible design, processing module 402 is specifically used for:
[0109] Define the data scatter values of the flag layer as (x i ,y i ,z i ), where x i Let y be the x-coordinate of the i-th data point. i Let z be the ordinate of the i-th data point. i Let z be the observation value of the i-th data point; i The corresponding trend value is
[0110] The first-order surface equation for the trend value is:
[0111] According to the least squares principle, the parameters b0, b1, and b2 are obtained by differentiation; where b0 is a number greater than 0, b1 is a number greater than 0, and b2 is a number greater than 0.
[0112] Substitute the parameters into the equation of a linear surface to determine the trend value.
[0113] In one possible design, the second determining module 403 is specifically used for:
[0114] According to the formula Identify and construct outliers;
[0115] Based on the geological information from the drilled wells and structural anomalies, the micro-topography is determined.
[0116] Figure 5 This is a schematic diagram illustrating the structure of the software platform according to an example embodiment of the present invention. Figure 5 As shown, this embodiment provides a software platform 50, including:
[0117] Processor 501; and,
[0118] Memory 502 is used to store executable instructions of the processor, and the memory may also be flash memory;
[0119] The processor 501 is configured to execute the various steps of the above method by executing executable instructions. See the relevant descriptions in the preceding method embodiments for details.
[0120] Alternatively, the memory 502 can be either standalone or integrated with the processor 501.
[0121] When the memory 502 is a device independent of the processor 501, the database 50 may also include:
[0122] Bus 503 is used to connect processor 501 and memory 502.
[0123] Furthermore, embodiments of this application also provide a computer-readable storage medium storing computer-executable instructions, wherein when at least one processor of a user device executes the computer-executable instructions, the user device performs the various possible methods described above.
[0124] Computer-readable media include computer storage media and communication media, wherein communication media include any medium that facilitates the transfer of computer programs from one place to another. Storage media can be any available medium accessible to a general-purpose or special-purpose computer. An exemplary storage medium is coupled to a 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. Alternatively, the ASIC can reside in a user device. Of course, the processor and storage medium can also exist as discrete components in a communication device.
[0125] Those skilled in the art will understand that all or part of the steps of the above-described method embodiments can be implemented by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When executed, the program performs the steps of the above-described method embodiments; and the aforementioned storage medium includes various media capable of storing program code, such as ROM, RAM, magnetic disks, or optical disks.
[0126] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.
Claims
1. A method of microtopographic profiling, characterized by, The method comprises the following steps: stratigraphic analysis is performed on the drilled well according to geological information of the drilled well, and stratigraphic analysis results are determined; spectrum analysis is performed on the drilled well according to seismic data of the drilled well, and a frequency band of seismic waves sensitive to fracture-vug carbonate reservoirs is extracted, and spectrum analysis results are determined; calibration is performed on the drilled well according to the stratigraphic analysis results and the spectrum analysis results, and calibration results are obtained; distances between each marker bed and the oil and gas reservoir are determined according to a position of the oil and gas reservoir; a marker bed with a minimum distance is selected for structural interpretation according to the calibration results; trend surface analysis is applied to process the structural interpretation, and a structural trend is determined; structural anomaly points are determined according to the structural interpretation and the structural trend, and the structural anomaly points are used to represent micro-landscapes.
2. The method of claim 1, wherein, The trend surface analysis is applied to process the structural interpretation, and the structural trend is determined, which comprises the following steps: The data scatter point value of the mark layer is defined as (x i ,y i ,z i ), wherein x i is the horizontal coordinate of the i th data, y i is the vertical coordinate of the i th data, and z i is the observation value of the i th data; the corresponding trend value of the observation value z i is ; The equation of the quadratic surface that constructs the trend value is: ; parameters b0, b1 and b2 are derived according to the least square principle, wherein b0 is a number greater than 0, b1 is a number greater than 0, and b2 is a number greater than 0; the parameters are substituted into the first-order surface equation to determine the trend value.
3. The method of claim 2, wherein, The structural anomaly points are determined according to the structural interpretation and the structural trend, and the structural anomaly points are used to represent micro-landscapes, which comprises the following steps: An abnormal point is determined according to the formula ; the micro-landscapes are determined according to the geological information of the drilled well and the structural anomaly points.
4. A microtopography profiling device, characterized by The method comprises the following steps: a first determination module is configured to calibrate the drilled well according to geological information and seismic data of the drilled well, and calibration results are obtained; structural interpretation of an oil and gas reservoir is determined according to the calibration results; a processing module is configured to apply trend surface analysis to process the structural interpretation, and a structural trend is determined; a second determination module is configured to determine structural anomaly points according to the structural interpretation and the structural trend, and the structural anomaly points are used to represent micro-landscapes; The first determination module is specifically configured to: stratigraphic analysis is performed on the drilled well according to the geological information, and stratigraphic analysis results are determined; spectrum analysis is performed on the drilled well according to the seismic data, and a frequency band of seismic waves sensitive to fracture-vug carbonate reservoirs is extracted, and spectrum analysis results are determined; calibration is performed on the drilled well according to the stratigraphic analysis results and the spectrum analysis results, and calibration results are obtained; distances between each marker bed and the oil and gas reservoir are determined according to a position of the oil and gas reservoir; a marker bed with a minimum distance is selected for structural interpretation according to the calibration results.
5. The apparatus of claim 4, wherein, The processing module is specifically configured to: The data scatter point value of the mark layer is defined as (x i ,y i ,z i ), wherein x i is the horizontal coordinate of the i th data, y i is the vertical coordinate of the i th data, and z i is the observation value of the i th data; the corresponding trend value of the observation value z i is ; The equation of the quadratic surface that constructs the trend value is: ; parameters b0, b1 and b2 are derived according to the least square principle, wherein b0 is a number greater than 0, b1 is a number greater than 0, and b2 is a number greater than 0; the parameters are substituted into the first-order surface equation to determine the trend value.
6. The apparatus of claim 5, wherein, The second determination module is specifically configured to: According to the formula , the abnormal point is determined; the micro-landscapes are determined according to the geological information of the drilled well and the structural anomaly points.
7. A configuration interpretation software platform, characterized by The method comprises the following steps: a processor; and a memory for storing executable instructions of the processor; wherein the processor is configured to execute the micro-landscape delineation method according to any one of claims 1 to 3 by executing the executable instructions.
8. A storage medium having stored thereon a computer program, characterized in that The program, when executed by the processor, implements the microtopographic delineation method of any one of claims 1 to 3.
Citation Information
Patent Citations
Method, device and equipment for analyzing distribution of remaining oil in fractured-vuggy carbonate oil and gas reservoir
CN118057213A