[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN118818613B - A method for dividing favorable areas based on shale gas exploration prediction - Google Patents

A method for dividing favorable areas based on shale gas exploration prediction Download PDF

Info

Publication number
CN118818613B
CN118818613B CN202411092604.3A CN202411092604A CN118818613B CN 118818613 B CN118818613 B CN 118818613B CN 202411092604 A CN202411092604 A CN 202411092604A CN 118818613 B CN118818613 B CN 118818613B
Authority
CN
China
Prior art keywords
data
class
water saturation
shale gas
speed
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202411092604.3A
Other languages
Chinese (zh)
Other versions
CN118818613A (en
Inventor
周德帅
张兴发
尹思奇
唐铈哲
蒋清山
杨阳
赵骥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Xinghui Kerui Technology Co ltd
Original Assignee
Chengdu Xinghui Kerui Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chengdu Xinghui Kerui Technology Co ltd filed Critical Chengdu Xinghui Kerui Technology Co ltd
Priority to CN202411092604.3A priority Critical patent/CN118818613B/en
Publication of CN118818613A publication Critical patent/CN118818613A/en
Application granted granted Critical
Publication of CN118818613B publication Critical patent/CN118818613B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开一种基于页岩气勘探预测的有利区划分方法,利用地震预测的I类连续储层厚度、含水饱和度、TOC、总含气量、孔隙度、脆性指数进行有利区的划分以及相关参数的选取,I类连续储层厚度及相关参数与产能密切相关,适用效果非常好,以往采用地震预测的优质页岩进行综合评价,从没有采用I类连续储层厚度及相关参数进行评价。另外以往页岩气没有采用地震预测含水饱和度,本发明新增含水饱和度参数,该参数在页岩气藏具有控制因数。总体来说该方法考虑I类连续储层厚度及相关参数,更加合理可靠,另外考虑含水饱和度参数参与评价,使得该发明更具有实用性可靠性。

The present invention discloses a method for dividing favorable areas based on shale gas exploration prediction, which uses the thickness of Class I continuous reservoirs predicted by seismic, water saturation, TOC, total gas content, porosity, and brittleness index to divide favorable areas and select related parameters. The thickness of Class I continuous reservoirs and related parameters are closely related to production capacity, and the applicable effect is very good. In the past, high-quality shales predicted by seismic were used for comprehensive evaluation, and the thickness of Class I continuous reservoirs and related parameters were never used for evaluation. In addition, shale gas did not use seismic prediction of water saturation in the past. The present invention adds a water saturation parameter, which has a control factor in shale gas reservoirs. In general, this method considers the thickness of Class I continuous reservoirs and related parameters, which is more reasonable and reliable. In addition, the water saturation parameter is considered in the evaluation, making the invention more practical and reliable.

Description

Advantageous region division method based on shale gas exploration prediction
Technical Field
The invention belongs to the technical field of exploration, and particularly relates to a beneficial zoning method based on shale gas exploration prediction.
Background
Shale gas is used as an unconventional natural gas which is rich in organic matters and is stored in the shale rich in organic matters, and development of the shale gas has very important significance for global energy structure transformation.
In oil exploration, accurately evaluating and dividing the favorable region has important significance for the development of oil and gas fields. For shale gas, the division of the beneficial zone is more important.
Conventional classification utilizes reservoir characteristic analysis, construction and fault evaluation, resource quantity estimation, resource quality evaluation, cost benefit analysis and the like, but the evaluation on shale gas and dense gas is not specific and is not closely related to productivity.
Disclosure of Invention
In order to solve the problems, the invention provides a beneficial zone dividing method based on shale gas exploration prediction, which evaluates shale gas and productivity in close relation to finish accurate division of beneficial zones and accurately guides production of shale gas in later period.
In order to achieve the aim, the technical scheme adopted by the invention is that the beneficial region dividing method based on shale gas exploration prediction comprises the following steps:
s100, acquiring three-dimensional seismic time domain gather data and well data of a work area and ground elevation data;
S201, performing time migration processing on the three-dimensional seismic time domain gather data to obtain time migration data of a time domain;
S202, converting three-dimensional seismic time domain gather data into angle gather data;
s203, obtaining the correlation of the density, the speed or the wave impedance related to the earthquake in the well data, and determining the optimal correlation parameter;
s300, obtaining an elastic parameter body through inversion by utilizing the horizon data of the target layer and combining with the optimal correlation parameters;
S401, interpreting the plane distribution of faults according to the phase characteristics of the horizon data of the target layer to obtain a fault plane distribution diagram, grading to obtain faults of all levels, converting the horizon data of the target layer into depth domain data by the product of an elastic parameter body and a speed body by utilizing the speed body in the elastic parameter body, and adding the depth domain data with ground elevation data to obtain a buried depth result;
S402, obtaining a class I continuous reservoir thickness result by adopting a sampling point statistical method;
S403, obtaining other parameters of the I-type continuous reservoir by using a fitting mode, wherein the other parameters comprise T0C, total air content, porosity and water saturation;
S404, calculating Young modulus and Poisson 'S ratio by using an elastic parameter body, and then obtaining a brittleness index by taking an average value of the Young modulus and Poisson' S ratio;
s500, dividing the reservoir into favorable areas with different levels by using the burial depth result, the class I continuous reservoir thickness result, the T0C, the total air content, the porosity and the water saturation brittleness index result value and the set reservoir division evaluation standard value.
Further, the horizon data of the target layer is picked up according to the wave group characteristics, the wave group relation and the reflection characteristics on the time offset data, and the horizon data of the target layer is obtained.
Further, obtaining the correlation of the density, the velocity or the wave impedance and the longitudinal and transverse wave velocity ratio related to the earthquake in the well data comprises the steps of selecting the well data of the I-type continuous reservoir section, carrying out histogram analysis on the well data, the density, the velocity or the wave impedance and the longitudinal and transverse wave velocity ratio, and determining which parameter or operation result of the parameter of the earthquake can separate the I-type continuous reservoir from the non-reservoir through the histogram analysis to be used as the optimal correlation parameter.
Further, by utilizing the horizon data of the target layer and combining with the optimal correlation parameters, obtaining an elastic parameter body through inversion, and determining that the elastic parameter body directly related to the well and the earthquake comprises density, longitudinal wave speed and transverse wave speed, wherein the longitudinal wave speed and transverse wave speed ratio is the operation result of the longitudinal wave speed or the transverse wave speed.
Further, inverting to obtain a speed body in the elastic parameter body, obtaining a target layer, adding the obtained horizon of the depth domain with the ground elevation to obtain the plane distribution of the interpretation fault of the burial depth result according to the phase characteristics of the target layer of the offset data, and obtaining a fault plane distribution map.
Further, when the sample point statistics method is adopted to obtain the class I continuous reservoir thickness result, the class I continuous reservoir thickness result is obtained by extracting the thickness of the sample point.
Further, other parameters of the type I continuous reservoir including T0C, total gas content, porosity and water saturation are obtained by using a fitting mode, wherein fitting is performed according to the uphole T0C, total gas content, porosity and water saturation and the uphole density or speed, which correlation degree of the uphole T0C, total gas content, porosity and water saturation and the speed or density is the highest is found, a corresponding relation is found, and the T0C, total gas content, porosity and water saturation are calculated by using the relation.
Further, the evaluation criteria of the favorable region is that the burial depth is less than 500 meters, the thickness of the class I continuous reservoir is more than or equal to 10 meters, the TOC is more than 3 percent, the total air content is more than 3 percent, the porosity is more than 4 percent, the water saturation is less than or equal to 40 percent, and the brittleness index is more than 55 percent;
The evaluation criteria of the second kind of favorable area are that the burial depth is less than 500 meters, the thickness of the class I continuous reservoir is less than 10 meters, TOC is more than 2 percent and less than 3 percent, the total air content is more than 2 percent and less than 3 percent, the porosity is more than 3 percent and less than 4 percent, the water saturation is more than 40 percent, and the brittleness index is more than 45 percent and less than 55 percent.
The beneficial effect of adopting this technical scheme is:
The method utilizes the thickness, the water saturation, the TOC, the total air content, the porosity and the brittleness index of the type I continuous reservoir predicted by the earthquake to divide the favorable region and select related parameters, the type I continuous reservoir thickness and the related parameters are closely related to the production energy, the application effect is very good, the comprehensive evaluation is carried out by adopting the high-quality shale predicted by the earthquake in the past, and the evaluation is carried out without adopting the type I continuous reservoir thickness and the related parameters. In addition, the shale gas does not adopt earthquake to predict the water saturation, and the invention increases the water saturation parameter which has a control factor in the shale gas reservoir. In general, the method considers the thickness of the I-type continuous reservoir and related parameters, is more reasonable and reliable, and additionally considers the water saturation parameters to participate in evaluation, so that the method has more practicability and reliability.
Drawings
FIG. 1 is a flow diagram of an advantageous compartmentalization method based on shale gas exploration prediction in accordance with the present invention;
FIG. 2 is a histogram of longitudinal wave velocity, transverse-to-transverse wave velocity ratio and density in an embodiment of the invention.
Detailed Description
The present invention will be further described with reference to the accompanying drawings, in order to make the objects, technical solutions and advantages of the present invention more apparent.
In this embodiment, referring to fig. 1, the invention provides a beneficial region division method based on shale gas exploration prediction, which includes the steps of:
s100, acquiring three-dimensional seismic time domain gather data and well data of a work area and ground elevation data;
S201, performing time migration processing on the three-dimensional seismic time domain gather data to obtain time migration data of a time domain;
Specifically, the horizon data of the target layer is picked up on the time offset data according to the wave group characteristics, the wave group relation and the reflection characteristics, and the horizon data of the target layer is obtained.
S202, converting the three-dimensional seismic time domain gather data into angle gather data.
S203, obtaining the correlation of the density, the speed or the wave impedance related to the earthquake in the well data, and determining the optimal correlation parameter;
Specifically, the method for acquiring the correlation of the density, the velocity or the wave impedance and the longitudinal and transverse wave velocity ratio related to the earthquake in the well data comprises the steps of selecting the well data of a type I continuous reservoir section, carrying out histogram analysis on the well data, the density, the velocity or the wave impedance and the longitudinal and transverse wave velocity ratio, and determining which parameter or operation result of the parameter of the well and the earthquake can separate the type I continuous reservoir from a non-reservoir through the histogram analysis to be used as the parameter for determining the optimal correlation.
S300, obtaining an elastic parameter body through inversion by utilizing the horizon data of the target layer and combining with the optimal correlation parameters;
And (3) obtaining an elastic parameter body through inversion by utilizing the horizon data of the target layer and combining with the optimal correlation parameter, and determining that the elastic parameter body directly related to the well and the earthquake comprises density, longitudinal wave speed and transverse wave speed, wherein the longitudinal wave speed and transverse wave speed ratio is the longitudinal wave speed or transverse wave speed operation result.
As shown in fig. 2, it is obvious from the graph that the density can better distinguish the type I continuous reservoir from the non-reservoir section, and the longitudinal wave velocity and the transverse wave velocity are more overlapped than the type I continuous reservoir and the non-reservoir section, and cannot be effectively separated.
S401, interpreting the plane distribution of faults according to the phase characteristics of the horizon data of the target layer to obtain a fault plane distribution diagram, grading to obtain faults of all levels, converting the horizon data of the target layer into depth domain data by using a velocity body in the elastic parameter body through the product of the elastic parameter body and the velocity body, and adding the depth domain data with ground elevation data to obtain a buried depth result.
Specifically, inverting to obtain a velocity body in the elastic parameter body, obtaining a target layer, adding the obtained horizon of the depth domain and the ground elevation to obtain the plane distribution of the interpretation fault of the burial depth result according to the phase characteristics of the target layer of the offset data, and obtaining a fault plane distribution map.
S402, obtaining a class I continuous reservoir thickness result by adopting a sampling point statistical method;
specifically, when the sample point statistics method is adopted to obtain the layer thickness result of the class I continuous reservoir, the class I continuous reservoir thickness result is obtained by extracting the thickness of the sample point and the longitudinal wave velocity body.
S403, obtaining other parameters of the I-type continuous reservoir by using a fitting mode, wherein the other parameters comprise T0C, total air content, porosity and water saturation;
Specifically, other parameters of the I-type continuous reservoir are obtained by using a fitting mode, including T0C, total air content, porosity and water saturation, wherein fitting is performed according to the uphole T0C, total air content, porosity and water saturation and the uphole density or speed, which correlation degree of the uphole T0C, total air content, porosity and water saturation and the speed or density is the highest is found, a corresponding relational expression is found, and the T0C, total air content, porosity and water saturation are calculated by using the relational expression.
S404, calculating Young modulus and Poisson 'S ratio by using an elastic parameter body, and then obtaining a brittleness index by taking an average value of the Young modulus and Poisson' S ratio;
specific:
Young's modulus:
Poisson ratio:
s500, dividing the reservoir into favorable areas with different levels by using the burial depth result, the class I continuous reservoir thickness result, the T0C, the total air content, the porosity and the water saturation brittleness index result value and the set reservoir division evaluation standard value.
Specific:
The evaluation standard of the favorable region is that the burial depth is less than 500 meters, the thickness of the I-type continuous reservoir is more than or equal to 10 meters, the TOC is more than 3 percent, the total air content is more than 3 percent, the porosity is more than 4 percent, the water saturation is less than or equal to 40 percent, and the brittleness index is more than 55 percent;
The evaluation criteria of the second kind of favorable area are that the burial depth is less than 500 meters, the thickness of the class I continuous reservoir is less than 10 meters, TOC is more than 2 percent and less than 3 percent, the total air content is more than 2 percent and less than 3 percent, the porosity is more than 3 percent and less than 4 percent, the water saturation is more than 40 percent, and the brittleness index is more than 45 percent and less than 55 percent.
The foregoing has shown and described the basic principles and main features of the present invention and the advantages of the present invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, and that the above embodiments and descriptions are merely illustrative of the principles of the present invention, and various changes and modifications may be made without departing from the spirit and scope of the invention, which is defined in the appended claims. The scope of the invention is defined by the appended claims and equivalents thereof.

Claims (8)

1. An advantageous zoning method based on shale gas exploration prediction is characterized by comprising the following steps:
s100, acquiring three-dimensional seismic time domain gather data and well data of a work area and ground elevation data;
S201, performing time migration processing on the three-dimensional seismic time domain gather data to obtain time migration data of a time domain;
S202, converting three-dimensional seismic time domain gather data into angle gather data;
s203, obtaining the correlation of the density, the speed or the wave impedance related to the earthquake in the well data, and determining the optimal correlation parameter;
s300, obtaining an elastic parameter body through inversion by utilizing the horizon data of the target layer and combining with the optimal correlation parameters;
S401, interpreting the plane distribution of faults according to the phase characteristics of the horizon data of the target layer to obtain a fault plane distribution diagram, grading to obtain faults of all levels, converting the horizon data of the target layer into depth domain data by the product of an elastic parameter body and a speed body by utilizing the speed body in the elastic parameter body, and adding the depth domain data with ground elevation data to obtain a buried depth result;
S402, obtaining a class I continuous reservoir thickness result by adopting a sampling point statistical method;
S403, obtaining other parameters of the I-type continuous reservoir by using a fitting mode, wherein the other parameters comprise T0C, total air content, porosity and water saturation;
S404, calculating Young modulus and Poisson 'S ratio by using an elastic parameter body, and then obtaining a brittleness index by taking an average value of the Young modulus and Poisson' S ratio;
s500, dividing the reservoir into favorable areas with different levels by using the burial depth result, the class I continuous reservoir thickness result, the T0C, the total air content, the porosity and the water saturation brittleness index result value and the set reservoir division evaluation standard value.
2. The method for distinguishing the favorable zone based on shale gas exploration prediction according to claim 1, wherein the horizon data of the target layer is picked up according to the wave group characteristics, the wave group relation and the reflection characteristics on the time offset data, and the horizon data of the target layer is obtained.
3. The method for distinguishing the favorable region based on the shale gas exploration prediction according to claim 2, wherein the method for acquiring the correlation of the density, the speed or the wave impedance and the longitudinal and transverse wave speed ratio related to the earthquake in the well data comprises the steps of selecting the well data of the type I continuous reservoir section, carrying out histogram analysis on the well data and the density, the speed or the wave impedance and the longitudinal and transverse wave speed ratio, and determining which parameter or operation result of the parameter of the well and the earthquake can separate the type I continuous reservoir from the non-reservoir through the histogram analysis as the optimal correlation parameter.
4. The method for distinguishing the favorable zones based on shale gas exploration prediction according to claim 3, wherein the elastic parameter body is obtained through inversion by utilizing the horizon data of the target layer and combining with the optimal correlation parameter, and the elastic parameter body directly related to the well and the earthquake is determined to comprise density, longitudinal wave speed and transverse wave speed, wherein the longitudinal wave speed and transverse wave speed ratio is the operation result of the longitudinal wave speed or the transverse wave speed.
5. The method for dividing the favorable region based on shale gas exploration prediction according to claim 1, wherein the velocity body in the elastic parameter body is obtained through inversion, a target layer is obtained, the obtained horizon of the depth domain is added with the ground elevation to obtain the buried depth achievement, the plane distribution of faults is interpreted according to the phase characteristics of the target layer of the offset data, and a fault plane distribution map is obtained.
6. The method for distinguishing the beneficial areas based on shale gas exploration prediction according to claim 1, wherein the class I continuous reservoir thickness results are obtained by extracting thickness of sampling points from longitudinal wave velocity bodies when the class I continuous reservoir thickness results are obtained by adopting a sampling point statistical method.
7. The method for distinguishing the beneficial zone based on shale gas exploration prediction according to claim 1, wherein the method for obtaining other parameters of the type I continuous reservoir by using a fitting mode comprises the steps of fitting the T0C, the total gas content, the porosity and the water saturation on the well with the density or the speed on the well according to the T0C, the total gas content, the porosity and the water saturation on the well, finding out which correlation degree among the T0C, the total gas content, the porosity and the water saturation on the well is the highest with the speed or the density, finding out a corresponding relation, and calculating the T0C, the total gas content, the porosity and the water saturation by using the relation.
8. The method for partitioning the beneficial zone based on shale gas exploration prediction according to claim 1, wherein the evaluation criteria of the class I beneficial zone are that the burial depth is less than 500 meters, the class I continuous reservoir is greater than or equal to 10 meters, the TOC is greater than 3%, the total air content is greater than 3%, the porosity is greater than 4%, the water saturation is less than or equal to 40%, the brittleness index is greater than 55%, the evaluation criteria of the class II beneficial zone are that the burial depth is less than 500 meters, the class I continuous reservoir is less than 10 meters, the TOC is greater than 2% and less than 3%, the total air content is greater than 2% and less than 3%, the porosity is greater than 3% and less than 4%, the water saturation is greater than 40%, and the brittleness index is greater than 45% and less than 55%.
CN202411092604.3A 2024-08-09 2024-08-09 A method for dividing favorable areas based on shale gas exploration prediction Active CN118818613B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202411092604.3A CN118818613B (en) 2024-08-09 2024-08-09 A method for dividing favorable areas based on shale gas exploration prediction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202411092604.3A CN118818613B (en) 2024-08-09 2024-08-09 A method for dividing favorable areas based on shale gas exploration prediction

Publications (2)

Publication Number Publication Date
CN118818613A CN118818613A (en) 2024-10-22
CN118818613B true CN118818613B (en) 2025-01-21

Family

ID=93076583

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202411092604.3A Active CN118818613B (en) 2024-08-09 2024-08-09 A method for dividing favorable areas based on shale gas exploration prediction

Country Status (1)

Country Link
CN (1) CN118818613B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108222925A (en) * 2017-12-08 2018-06-29 中国石油集团川庆钻探工程有限公司 Shale gas reservoir classification comprehensive evaluation method
CN113534253A (en) * 2020-04-22 2021-10-22 中国石油天然气集团有限公司 Shale gas three-dimensional seismic sweet spot optimization method and device

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2506039A3 (en) * 2011-03-28 2013-08-14 Conocophillips Company Methods for Seismic Fracture Parameter Estimation and Gas Filled Fracture Identification From Vertical Well Log Data
US11519262B2 (en) * 2017-01-17 2022-12-06 Schlumberger Technology Corporation Systematic evaluation of shale plays

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108222925A (en) * 2017-12-08 2018-06-29 中国石油集团川庆钻探工程有限公司 Shale gas reservoir classification comprehensive evaluation method
CN113534253A (en) * 2020-04-22 2021-10-22 中国石油天然气集团有限公司 Shale gas three-dimensional seismic sweet spot optimization method and device

Also Published As

Publication number Publication date
CN118818613A (en) 2024-10-22

Similar Documents

Publication Publication Date Title
CN109102180B (en) Comprehensive parameter evaluation method for double-dessert evaluation of tight sandstone reservoir
CN111461386B (en) Shale gas dessert prediction method based on BP neural network
CN109653725A (en) A layer water flooding degree log interpretation method is stored up based on sedimentary micro and the mixed of rock phase
CN109541685B (en) River channel sand body identification method
CN105528522A (en) Method and device for calculating resource quantity of continuous oil and gas reservoir based on spatial grid
CN111340379A (en) Dissection method for deep sea water channel sand body building structure in thin well network area
CN107829731A (en) A kind of volcanic rock porosity correction method of clay alteration
CN109031436A (en) A kind of hypotonic loose sand reservoir water saturation evaluation method of low-resistance based under Reservoir Classification
Chen et al. Geological risk mapping and prospect evaluation using multivariate and Bayesian statistical methods, western Sverdrup Basin of Canada
CN112051626A (en) Quantitative characterization method of source-sink system
CN113361111B (en) Partitioning method of low-permeability reservoir characteristic model
CN108805158B (en) Compact oil reservoir lithogenous phase dividing method
CN118818613B (en) A method for dividing favorable areas based on shale gas exploration prediction
Iltaf et al. Facies and petrophysical modeling of Triassic Chang 6 tight sandstone reservoir, Heshui oil field, Ordos basin, China
CN112906465A (en) Coal measure stratum acoustic curve reconstruction method and system based on stratum factors
CN110847887B (en) Method for identifying and evaluating cracks of fine-grain sedimentary continental facies shale
Nolan et al. Spatial variability of groundwater recharge and its effect on shallow groundwater quality in southern New Jersey
CN111413737A (en) Method for explaining seismic horizon to deposition unit in close well pattern area
CN113919140A (en) A method for evaluating the static connectivity of thin interbedded oil reservoirs in fan deltas with large well spacing
CN116335648A (en) Method for quantitatively evaluating dessert quality of conglomerate horizontal well based on logging while drilling data
CN110795513B (en) Method for predicting distribution of river facies source storage ectopic type compact oil gas dessert area
CN112983409A (en) Method for predicting mineral components of clastic rock by using machine learning algorithm
CN112415596A (en) Dolomite structure type identification method based on logging information
CN115126479B (en) Method, device, equipment and storage medium for determining oil and gas reservoir exploration position
CN115099578B (en) Oil and gas resource evaluation method for oil-containing gas basin edge region

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant