CN105350959A - Method for determining gas saturation of shale gas reservoir through well-logging lithologic density - Google Patents
Method for determining gas saturation of shale gas reservoir through well-logging lithologic density Download PDFInfo
- Publication number
- CN105350959A CN105350959A CN201510753459.3A CN201510753459A CN105350959A CN 105350959 A CN105350959 A CN 105350959A CN 201510753459 A CN201510753459 A CN 201510753459A CN 105350959 A CN105350959 A CN 105350959A
- Authority
- CN
- China
- Prior art keywords
- well
- logging
- gas
- density
- lithology
- 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a method for determining gas saturation of shale gas reservoir through well-logging lithologic density. The method comprises following steps: obtaining gas saturation of core test analyses, core volume and density and well-logging natural gamma at the corresponding depth and utilizing a least square method to determine model coefficients of a first model and a second model in order to ultimately determine an optimal model; acquiring natural gamma and lithologic density of a well to be interpreted, acquiring total hydrocarbon and methane content of the well to be interpreted through gas-logging resources and acquiring lithic fragments, cores, lithology of geological logging to be explained through geological-logging resources; dividing the shale gas reservoir through acquired gas-logging and mud-logging resources; calculating gas saturation Sg of the reservoir of the well to be interpreted through the determined optimal model and well-logging resources; and outputting a calculation result.The method for determining gas saturation of the shale gas reservoir through well-logging lithologic density has following beneficial effects: the method is easy and convenient and is extensively applied and is ready to use in 120 wells of shale gas; and calculated gas saturation of the shale gas reservoir is approximate to gas saturation of the shale gas reservoir, with the average error less than 10%.
Description
Technical field
The present invention relates to a kind of method by WELL LITHOLOGY density determination shale gas reservoir gas-bearing saturation ratio
Background technology
The success of shale gas is developed and is brought new forms of energy revolution to the whole world; the utilization of shale gas is to environmental protection and energy-conservationly all have great significance; the searching of shale gas-bearing formation is more and more subject to people's attention with discovery, and effective exploitation shale gas becomes the new normality of one of mineral exploration and exploitation.Shale gas gas saturation is an indispensable key parameter of shale gas reservoir, shale gas reservoir and shale gas selection and appraisal of exploration area.
Traditional shale gas reservoir gas-bearing saturation ratio defining method is that employing ECS element captures well logging method or A Erqifa asks for.At present, ECS element captures well logging can only rely on the service of external Logging Company, and cost is very high.A Erqifa asks for shale gas reservoir gas-bearing saturation ratio, mainly uses data such as the well logging of shale reservoir resistivity, degree of porosity etc.A Erqifa asks for shale gas reservoir gas-bearing saturation ratio precondition, only layer water conduction, and the mineralogical composition forming formation skeleton is non-conductive.But there is the conductive mineral of pyrite and so in shale gas reservoir, make shale gas reservoir resistivity extremely diminish, the gas saturation degree of accuracy calculated is lower more, and error is comparatively large, and the method has certain limitation.
Summary of the invention
The object of the invention is for the above-mentioned state of the art, be intended to rely on core experiment to provide a kind of by litho-density log data determination shale gas reservoir gas-bearing saturation computation model, the method containing the shale gas reservoir gas-bearing saturation ratio precision under low-resistivity interlayer condition can be improved.
The implementation of the object of the invention is, by the method for WELL LITHOLOGY density determination shale gas reservoir gas-bearing saturation ratio, concrete steps are:
1) the well logging natural gamma GR of work area fixed well core sample test analysis gas saturation Sg, rock core bulk density ρ b and the corresponding degree of depth is obtained; The sample chosen must be same work area, identical series of strata, and quantity can not be less than 12 pieces, and well number is not limit; Then utilize the model coefficient of least square method Confirming model one Sg1=a/ ρ b+b, model two Sg2=a ρ b+bGR+c, finally determine best model;
In formula, rock core bulk density ρ b dimension is g/cm
3,
Natural gamma GR dimension is API,
Gas saturation Sg dimension is %;
2) obtained natural gamma GR, the lithology and density logging DEN of well to be explained by well-log information, obtained total hydrocarbon, the methane content of well to be explained by gas measuring logging data, obtained landwaste, the rock core lithology of geological logging to be explained by geological logging data; Lithology and density logging DEN dimension is g/cm
3;
3) by step 2) well logging that obtains and log data divide shale gas reservoir;
4) by step 2) borehole logging tool lithology and density logging (DEN) alternative steps 1 to be explained) rock core bulk density (ρ b), calculate well reservoir gas-bearing saturation ratio Sg to be explained by the best model determined and well-log information;
5) result of calculation is exported.
The invention provides a kind of low-resistance shale reservoir and calculate formation gas bearing saturation ratio solution, determine that shale reservoir gas-bearing saturation intent is easier compared with A Erqifa, the scope of application is wider.
The present invention has applied 120 mouthfuls of wells and has used in certain shale gas field, the shale gas-bearing formation gas saturation that the shale gas reservoir gas-bearing saturation ratio calculated and core test obtain is close, average error is no more than 10%, can meet the on-the-spot needs surveyed well logging evaluating reservoir and calculate shale gas gas saturation.
Accompanying drawing explanation
Fig. 1 is workflow block diagram of the present invention,
Fig. 2 is model one inverse saturation ratio of the present invention and rock core saturation ratio cross plot,
Fig. 3 is model two inverse saturation ratio of the present invention and rock core saturation ratio cross plot,
Fig. 4 is J shale gas field S well application example figure of the present invention.
Detailed description of the invention
The applicant be experimental studies have found that by shale gas reservoir core, there is obvious quantitative relationship in shale gas reservoir core bulk density (ρ b) and its gas saturation (Sg), this relation shows as gas saturation and becomes positive correlation with the inverse of rock core bulk density, and gas saturation equally with rock core bulk density and corresponding well logging natural gamma linear.WELL LITHOLOGY density (DEN) and rock core bulk density are numerically equal, and just metering system is different.Therefore, WELL LITHOLOGY density is same exists above-mentioned quantitative relationship with shale gas reservoir gas-bearing saturation ratio, utilizes density of logging well, natural gamma (GR) data quantitatively can determine shale gas reservoir gas-bearing saturation ratio (Sg).With the calculating of rock core information, the Sg=a1/ ρ b+b1 that uses a model, computation model two Sg=a2 ρ b+b2GR+c2; Based on the calculating of well-log information, model one Sg=a1/DEN+b1, computation model two Sg=a2DEN+b2GR+c2.
In formula, a1, b1, a2, b2, c2 are equation coefficient.
With reference to the accompanying drawings in detail the present invention is described in detail.
With reference to Fig. 1, concrete steps of the present invention are:
1) the well logging natural gamma GR of work area fixed well core sample test analysis gas saturation Sg, rock core bulk density ρ b and the corresponding degree of depth is obtained; The sample chosen must be same work area, identical series of strata, and quantity can not be less than 12 pieces, and well number is not limit; Then utilize the model coefficient of least square method Confirming model one Sg1=a/ ρ b+b, model two Sg2=a ρ b+bGR+c, and then the gas saturation utilizing model inverse to obtain and core test saturation ratio return, what index of correlation was high is best model.
2) obtained natural gamma GR, the lithology and density logging DEN of well to be explained by well-log information, obtained total hydrocarbon, the methane content of well to be explained by gas measuring logging data, obtained landwaste, the rock core lithology of geological logging to be explained by geological logging data;
Natural gamma GR, the lithology and density logging DEN data depth of well to be explained are spaced apart 0.1m ~ 1.0m, save as wis or txt form; Gas detection logging total hydrocarbon and methane content depth interval are 1m or 0.5m, save as wis or txt form; The lithology data of geologic debris, rock core saves as excel or txt form.
3) by step 2) well logging that obtains and log data divide shale gas reservoir;
Shale gas reservoir division foundation: lithology is shale or argiloid; Total hydrocarbon, methane exist obviously abnormal, not containing heavy hydrocarbon composition; High natural gamma, low-density.
4) by step 2) borehole logging tool lithology and density logging (DEN) alternative steps 1 to be explained) rock core bulk density (ρ b), calculate well reservoir gas-bearing saturation ratio Sg to be explained by the best model determined and well-log information;
5) export result of calculation, usually, the gas saturation Sg of shale gas-bearing formation is not less than 0.5.
With J shale gas field S well specific embodiment in detail the present invention is described in detail below.
1) gas saturation (Sg) of 4 mouthfuls of shale wells totally 82 pieces of shale core is obtained by the core test analysis of J shale gas field, rock core bulk density (ρ b) and obtained the natural gamma (GR) of the corresponding degree of depth by well-log information, then by 82 pieces of shale core Sg and ρ b and GR, do least square method by model one Sg=a1/ ρ b+b1 and model two Sg=a2 ρ b+b2GR+c2 and return process, obtain model coefficient a1=1610.600, b1=567.240, a2=-262.896, b2=0.183, c2=768.147, i.e. model one Sg=1610.6/ ρ b+567.24, model two Sg=-262.896 ρ b-0.183GR+768.147, Y-axis is done again by the gas saturation of model one inverse, core test gas saturation is done X-axis and is carried out matching, coefficient R is 0.79 (see Fig. 2), model two inverse gas saturation does Y-axis, rock core saturation ratio is done X-axis and returned, coefficient R is 0.81 (see Fig. 3), the index of correlation of contrast model one and model two, and model two index of correlation is higher than model one, model two is best model, selects model two.
2) obtained natural gamma GR, the lithology and density logging DEN of S well to be explained by well-log information, obtained total hydrocarbon, the methane content of JingSJing to be explained by gas measuring logging data, obtained landwaste, the rock core lithology of geological logging to be explained by geological logging data;
3) by step 2) the lithology classification S well shale gas reservoir of the natural gamma of the S well that obtains, lithology and density logging, total hydrocarbon, methane and landwaste rock core, S well shale gas mainly shows interval at 2170.0 ~ 2220.0m;
4) by step 2) utilize density and the natural gamma material computation S well reservoir gas-bearing saturation ratio Sg of best model Sg=-262.896 ρ b+0.183GR+768.147 and S well;
5) result of calculation is exported, it is 50.0% ~ 80.0% (see Fig. 4) that S well 2170.0 ~ 2220.0m well section calculates gas saturation substantial scope, contrast with this well rock core saturation ratio test result, average error is less than 10%, can meet shale gas reservoir and shale gas Gas reservoir eraluation needs.
Claims (4)
1., by the method for WELL LITHOLOGY density determination shale gas reservoir gas-bearing saturation ratio, it is characterized in that concrete steps are:
1) the well logging natural gamma GR of work area fixed well core sample test analysis gas saturation Sg, rock core bulk density ρ b and the corresponding degree of depth is obtained; The sample chosen must be same work area, identical series of strata, and quantity can not be less than 12 pieces, and well number is not limit; Then utilize the model coefficient of least square method Confirming model one Sg1=a/ ρ b+b, model two Sg2=a ρ b+bGR+c, finally determine best model;
In formula, rock core bulk density ρ b dimension is g/cm
3,
Natural gamma GR dimension is API,
Gas saturation Sg dimension is %;
2) obtained natural gamma GR, the lithology and density logging DEN of well to be explained by well-log information, obtained total hydrocarbon, the methane content of well to be explained by gas measuring logging data, obtained landwaste, the rock core lithology of geological logging to be explained by geological logging data; Lithology and density logging DEN dimension is g/cm
3;
3) by step 2) well logging that obtains and log data divide shale gas reservoir;
4) by step 2) borehole logging tool lithology and density logging (DEN) alternative steps 1 to be explained) rock core bulk density (ρ b), calculate well reservoir gas-bearing saturation ratio Sg to be explained by the best model determined and well-log information;
5) result of calculation is exported.
2. the method by WELL LITHOLOGY density determination shale gas reservoir gas-bearing saturation ratio according to claim 1, it is characterized in that step 1) in utilize model inverse to obtain gas saturation and core test saturation ratio return, what index of correlation was high is best model.
3. the method by WELL LITHOLOGY density determination shale gas reservoir gas-bearing saturation ratio according to claim 1, is characterized in that step 2) in the natural gamma GR of well to be explained, lithology and density logging DEN data depth be spaced apart 0.1m ~ 1.0m; Gas detection logging total hydrocarbon and methane content depth interval are 1m or 0.5m.
4. the method by WELL LITHOLOGY density determination shale gas reservoir gas-bearing saturation ratio according to claim 1, is characterized in that step 3) shale gas reservoir division foundation: lithology is shale or argiloid; Total hydrocarbon, methane exist obviously abnormal, not containing heavy hydrocarbon composition; High natural gamma, low-density.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510753459.3A CN105350959B (en) | 2015-11-06 | 2015-11-06 | The method that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510753459.3A CN105350959B (en) | 2015-11-06 | 2015-11-06 | The method that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105350959A true CN105350959A (en) | 2016-02-24 |
| CN105350959B CN105350959B (en) | 2018-03-16 |
Family
ID=55327005
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510753459.3A Active CN105350959B (en) | 2015-11-06 | 2015-11-06 | The method that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105350959B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108561126A (en) * | 2018-03-13 | 2018-09-21 | 中石化石油工程技术服务有限公司 | A kind of simple and easy method of determining shale gas reservoir organic porosity |
| CN108647417A (en) * | 2018-04-28 | 2018-10-12 | 中石化石油工程技术服务有限公司 | A kind of simple and easy method of determining shale gas reservoir gas-bearing saturation degree |
| CN109063296A (en) * | 2018-07-23 | 2018-12-21 | 中国石油集团川庆钻探工程有限公司 | Shale gas content while-drilling calculation method |
| CN112649858A (en) * | 2019-10-11 | 2021-04-13 | 中国石油化工股份有限公司 | Shale brittleness prediction method and system based on core test |
| CN115749760A (en) * | 2022-11-28 | 2023-03-07 | 中海石油(中国)有限公司海南分公司 | Reservoir fluid property evaluation method based on measurement and recording combination |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101487390A (en) * | 2009-02-23 | 2009-07-22 | 大庆油田有限责任公司 | Archie mode method for confirming initial oil saturation of oil layer |
| CN102748004A (en) * | 2012-06-06 | 2012-10-24 | 中国石油化工股份有限公司 | Method for determining hydrocarbon saturation of gas logging |
| US8606522B2 (en) * | 2008-09-30 | 2013-12-10 | Schlumberger Technology Corporation | Method to determine current gas saturation in a near-wellbore zone in a volatile oil formation |
| CN103912269A (en) * | 2014-03-28 | 2014-07-09 | 中石化江汉石油工程有限公司测录井公司 | Method for determining formation fracture pressure gradient logging of shale gas reservoir |
| CN103912268A (en) * | 2014-03-28 | 2014-07-09 | 中石化江汉石油工程有限公司测录井公司 | Shale reservoir gas saturation determining method based on TOC |
-
2015
- 2015-11-06 CN CN201510753459.3A patent/CN105350959B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8606522B2 (en) * | 2008-09-30 | 2013-12-10 | Schlumberger Technology Corporation | Method to determine current gas saturation in a near-wellbore zone in a volatile oil formation |
| CN101487390A (en) * | 2009-02-23 | 2009-07-22 | 大庆油田有限责任公司 | Archie mode method for confirming initial oil saturation of oil layer |
| CN102748004A (en) * | 2012-06-06 | 2012-10-24 | 中国石油化工股份有限公司 | Method for determining hydrocarbon saturation of gas logging |
| CN103912269A (en) * | 2014-03-28 | 2014-07-09 | 中石化江汉石油工程有限公司测录井公司 | Method for determining formation fracture pressure gradient logging of shale gas reservoir |
| CN103912268A (en) * | 2014-03-28 | 2014-07-09 | 中石化江汉石油工程有限公司测录井公司 | Shale reservoir gas saturation determining method based on TOC |
Non-Patent Citations (3)
| Title |
|---|
| 任元等: "基于TOC的页岩储层含气饱和度确定方法及应用", 《江汉石油职工大学学报》 * |
| 王濡岳等: "页岩气储层地球物理测井评价研究现状", 《地球物理学进展》 * |
| 石文睿等: "涪陵页岩气田焦石坝页岩气储层含气量测井评价", 《测井技术》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108561126A (en) * | 2018-03-13 | 2018-09-21 | 中石化石油工程技术服务有限公司 | A kind of simple and easy method of determining shale gas reservoir organic porosity |
| CN108647417A (en) * | 2018-04-28 | 2018-10-12 | 中石化石油工程技术服务有限公司 | A kind of simple and easy method of determining shale gas reservoir gas-bearing saturation degree |
| CN109063296A (en) * | 2018-07-23 | 2018-12-21 | 中国石油集团川庆钻探工程有限公司 | Shale gas content while-drilling calculation method |
| CN112649858A (en) * | 2019-10-11 | 2021-04-13 | 中国石油化工股份有限公司 | Shale brittleness prediction method and system based on core test |
| CN115749760A (en) * | 2022-11-28 | 2023-03-07 | 中海石油(中国)有限公司海南分公司 | Reservoir fluid property evaluation method based on measurement and recording combination |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105350959B (en) | 2018-03-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lai et al. | Brittleness index estimation in a tight shaly sandstone reservoir using well logs | |
| CN104278991B (en) | Saline Lake Facies hydrocarbon source rock organic carbon and the polynary well logging computational methods of hydrocarbon potential | |
| CN102175832B (en) | Method for determining optimal saturation calculation model of typical reservoir | |
| CN106950347B (en) | A method of evaluation mud shale each group partial volume | |
| CN103122762B (en) | Method and device for detecting effective fractured interval of unconventional shale oil and gas reservoir | |
| CN104514552B (en) | A kind of method that coalbed methane reservoir identification is predicted with abundance | |
| CN103470250B (en) | A kind of method and apparatus measuring formation pore structure and fluid behaviour | |
| CN105134195A (en) | Shale gas reservoir quality evaluation method based on logging information | |
| CN107795320B (en) | Calculation method for horizontal well carbonate reservoir parameters | |
| CN103698811A (en) | Carbonate rock structure component well logging quantitative identification method and application thereof | |
| Ali et al. | Hydrocarbon potential assessment of carbonate-bearing sediments in a meyal oil field, Pakistan: Insights from logging data using machine learning and quanti elan modeling | |
| CN104570128A (en) | Unconventional oil and gas reservoir predicting method based on rock brittleness factors | |
| Li et al. | A rock physics model for estimating elastic properties of upper Ordovician-lower Silurian mudrocks in the Sichuan Basin, China | |
| CN105350959A (en) | Method for determining gas saturation of shale gas reservoir through well-logging lithologic density | |
| Sanderson et al. | Quantitative analysis of tin-and tungsten-bearing sheeted vein systems | |
| CN105064986A (en) | Method for building reservoir four-property relationship spectrum by using conventional well detection and logging information | |
| CN105242307A (en) | Complex carbonate stratum earthquake porosity obtaining method and apparatus | |
| Guoqiang | Challenges and countermeasures of log evaluation in unconventional petroleum exploration and development | |
| CN105205296A (en) | Method for calculating porosity of shale gas reservoir | |
| Feng et al. | Accurate determination of water saturation in tight sandstone gas reservoirs based on optimized Gaussian process regression | |
| CN110529106B (en) | Method for determining content of coal seam micro-components by using logging information | |
| CN105370270A (en) | Method for determining gas saturation of shale gas reservoir by longitudinal-transverse wave time difference of dipole acoustic waves | |
| Zappone et al. | SAPHYR: Swiss Atlas of Physical Properties of Rocks: the continental crust in a database | |
| CN104514553A (en) | Recognition method of gas reservoir in lithic sandstone layer and application thereof | |
| Zhang et al. | Estimation of total porosity in shale formations from element capture logging and conventional logging data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200812 Address after: 100728 Beijing, Chaoyangmen, North Street, No. 22, No. Co-patentee after: SINOPEC OILFIELD SERVICE Corp. Patentee after: SINOPEC Group Co-patentee after: SINOPEC JIANGHAN OILFIELD SERVICE Corp. Address before: 100029, Beijing, Chaoyang District Xin Tung Street, a No. 6, Twelfth layers Co-patentee before: Logging company of Sinopec Jianghan Petroleum Engineering Co.,Ltd. Patentee before: SINOPEC OILFIELD SERVICE Corp. |
|
| TR01 | Transfer of patent right |