CN105350959B - The method that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density - Google Patents

The method that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density Download PDF

Info

Publication number
CN105350959B
CN105350959B CN201510753459.3A CN201510753459A CN105350959B CN 105350959 B CN105350959 B CN 105350959B CN 201510753459 A CN201510753459 A CN 201510753459A CN 105350959 B CN105350959 B CN 105350959B
Authority
CN
China
Prior art keywords
well
gas
density
lithology
logging
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
CN201510753459.3A
Other languages
Chinese (zh)
Other versions
CN105350959A (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.)
China Petrochemical Corp
Sinopec Oilfield Service Corp
Sinopec Jianghan Petroleum Engineering Co Ltd
Original Assignee
Sinopec Oilfield Service Corp
Logging Co of Sinopec Jianghan Petroleum Engineering 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 Sinopec Oilfield Service Corp, Logging Co of Sinopec Jianghan Petroleum Engineering Co Ltd filed Critical Sinopec Oilfield Service Corp
Priority to CN201510753459.3A priority Critical patent/CN105350959B/en
Publication of CN105350959A publication Critical patent/CN105350959A/en
Application granted granted Critical
Publication of CN105350959B publication Critical patent/CN105350959B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The present invention relates to a kind of method that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density, obtain the well logging natural gamma of core test analysis gas saturation, rock core bulk density and corresponding depth, the model coefficient of model one, model two is determined using least square method, finally determines best model;Natural gamma, the lithology and density logging of well to be explained are obtained by well-log information, total hydrocarbon, the methane content of well to be explained are obtained by gas measuring logging data, landwaste, the rock core lithology of geological logging to be explained are obtained by geological logging data;Shale gas reservoir is divided by the well logging of acquisition and log data;The best model and well-log information of determination calculate well reservoir gas-bearing saturation degree Sg to be explained;Export result of calculation.Of the invention easy, the scope of application is wider.The present invention uses in certain shale gas field using 120 mouthfuls of wells, and the shale gas-bearing formation gas saturation that the shale gas reservoir gas-bearing saturation degree calculated obtains with core test approaches, and mean error is no more than 10%.

Description

The method that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density
Technical field
The present invention relates to a kind of method that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density.
Background technology
The successful exploitation of shale gas brings new energy revolution to the whole world, and the utilization of shale gas to environmental protection and saves all Have great significance, the searching of shale gas-bearing formation is increasingly valued by people with finding, effective exploitation shale gas turns into ore deposit Produce a kind of new normality of exploitation.Shale gas gas saturation is the one of shale gas reservoir, shale gas reservoir and shale gas selection and appraisal of exploration area The indispensable key parameter of item.
Traditional shale gas reservoir gas-bearing saturation degree determines that method is asked using ECS elements capture well logging method or A Erqifa Take.At present, ECS elements capture well logging can only rely on external Logging Company's service, and cost is very high.A Erqifa asks for shale gas Reservoir gas-bearing saturation degree, mainly using data such as shale reservoir well logging resistivity, porositys.A Erqifa asks for shale gas storage Layer gas saturation precondition is that only stratum water is conductive, and the mineralogical composition for forming formation skeleton is non-conductive.However, shale The conductive minerals that pyrite etc be present, make shale gas reservoir resistivity diminish extremely, the gassiness being calculated more in gas reservoir The saturation degree degree of accuracy is relatively low, and error is larger, and this method has certain limitation.
The content of the invention
The purpose of the present invention is to be directed to the above-mentioned state of the art, it is intended to relies on core experiment to provide one kind by litho-density log Data determines shale gas reservoir gas-bearing saturation computation model, it is possible to increase the shale gas reservoir under the conditions of interlayer containing low-resistivity The method of gas saturation precision.
The implementation of the object of the invention is the method that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density, Concretely comprise the following steps:
1) work area fixed well core sample test analysis gas saturation Sg, rock core bulk density ρ b are obtained and is corresponded to deep The well logging natural gamma GR of degree;The sample of selection must be same work area, identical series of strata, and quantity can not be less than 12 pieces, well number It is unlimited;Then Sg1=a/ ρ b+b of model one, the Sg2=a ρ b+bGR+c of model two model are determined using least square method Coefficient, finally determine best model;
Rock core bulk density ρ b dimensions are g/cm in formula3,
Natural gamma GR dimensions are API,
Gas saturation Sg dimensions are %;
2) natural gamma GR, the lithology and density logging DEN of well to be explained are obtained by well-log information, is obtained by gas measuring logging data Total hydrocarbon, the methane content of well to be explained are taken, landwaste, the rock core lithology of geological logging to be explained are obtained by geological logging data; Lithology and density logging DEN dimensions are g/cm3
3) well logging obtained by step 2) and log data divide shale gas reservoir;
4) step 2) borehole logging tool lithology and density logging (DEN) alternative steps 1 to be explained are used) rock core bulk density (ρ b), by true Fixed best model and well-log information calculate well reservoir gas-bearing saturation degree Sg to be explained;
5) result of calculation is exported.
The invention provides a kind of low-resistance shale reservoir to calculate formation gas bearing saturation degree solution, is determined compared with A Erqifa Shale reservoir gas saturation method is easier, and the scope of application is wider.
The present invention uses in certain shale gas field using 120 mouthfuls of wells, the shale gas reservoir gas-bearing saturation degree and rock calculated The shale gas-bearing formation gas saturation that heart test obtains approaches, and mean error is no more than 10%, can meet that scene is surveyed well logging reservoir and commented Valency and the needs for calculating shale gas gas saturation.
Brief description of the drawings
Fig. 1 is workflow block diagram of the present invention,
Fig. 2 is the inverse saturation degree of model one of the present invention and rock core saturation degree cross plot,
Fig. 3 is the inverse saturation degree of model two of the present invention and rock core saturation degree cross plot,
Fig. 4 is S well exemplary application maps in J shale gas field of the present invention.
Embodiment
The applicant by shale gas reservoir core experimental studies have found that, shale gas reservoir core bulk density (ρ b) and its Obvious quantitative relationship be present in gas saturation (Sg), this relation show as gas saturation and rock core bulk density it is reciprocal into Positive correlation, and gas saturation is equally linear with rock core bulk density and corresponding well logging natural gamma.Survey Well lithology and density logging (DEN) is numerically equal with rock core bulk density, and simply metering system is different.Therefore, WELL LITHOLOGY Equally there is above-mentioned quantitative relationship in density, utilize density of logging well, natural gamma (GR) data with shale gas reservoir gas-bearing saturation degree Shale gas reservoir gas-bearing saturation degree (Sg) can be quantitatively determined.With the calculating of rock core information, using Sg=a1/ ρ b+b1 of model one, The Sg=a2 ρ b+b2GR+c2 of computation model two;Calculating based on well-log information, the Sg=a1/DEN+b1 of model one, calculate mould The Sg=a2DEN+b2GR+c2 of type two.
A1, b1, a2, b2, c2 are equation coefficient in formula.
The present invention is described in detail with reference to the accompanying drawings.
Reference picture 1, of the invention concretely comprises the following steps:
1) work area fixed well core sample test analysis gas saturation Sg, rock core bulk density ρ b are obtained and is corresponded to deep The well logging natural gamma GR of degree;The sample of selection must be same work area, identical series of strata, and quantity can not be less than 12 pieces, well number It is unlimited;Then Sg1=a/ ρ b+b of model one, the Sg2=a ρ b+bGR+c of model two model are determined using least square method Coefficient, then recycle the obtained gas saturation of model inverse to be returned with core test saturation degree, coefficient correlation it is high for most Good model.
2) natural gamma GR, the lithology and density logging DEN of well to be explained are obtained by well-log information, is obtained by gas measuring logging data Total hydrocarbon, the methane content of well to be explained are taken, landwaste, the rock core lithology of geological logging to be explained are obtained by geological logging data;
The natural gamma GR of well to be explained, lithology and density logging DEN data depths at intervals of 0.1m~1.0m, save as wis or Txt forms;Gas detection logging total hydrocarbon and methane content depth interval are 1m or 0.5m, save as wis or txt forms;Geologic debris, The lithology data of rock core saves as excel or txt forms.
3) well logging obtained by step 2) and log data divide shale gas reservoir;
Shale gas reservoir division foundation:Lithology is shale or argiloid;There is obvious exception in total hydrocarbon, methane, without weight Hydrocarbon composition;High natural gamma, low-density.
4) step 2) borehole logging tool lithology and density logging (DEN) alternative steps 1 to be explained are used) rock core bulk density (ρ b), by true Fixed best model and well-log information calculate well reservoir gas-bearing saturation degree Sg to be explained;
5) result of calculation is exported, generally, the gas saturation Sg of shale gas-bearing formation is not less than 0.5.
The present invention is described in detail with J shale gas field S wells specific embodiment below.
1) gas saturation for obtaining 4 mouthfuls of shale wells totally 82 pieces of shale cores is analyzed by J shale gas field core test (Sg), rock core bulk density (ρ b) and the natural gamma (GR) by the corresponding depth of well-log information acquisition, then by 82 pieces of pages Rock rock core Sg and ρ b and GR, least square method is done by the Sg=a1/ ρ b+b1 of model one and the Sg=a2 ρ b+b2GR+c2 of model two Recurrence is handled, and obtains model coefficient a1=1610.600, b1=567.240, a2=-262.896, b2=0.183, c2= 768.147, i.e. Sg=1610.6/ ρ b+567.24 of model one, the Sg=-262.896 ρ b-0.183GR+ of model two 768.147, then with the gas saturation of the inverse of model one do Y-axis, core test gas saturation is done X-axis and is fitted, related Coefficients R is 0.79 (see Fig. 2);The inverse gas saturation of model two does Y-axis, rock core saturation degree is done X-axis and returned, coefficient correlation R is 0.81 (see Fig. 3), contrast model one and model two coefficient correlation, and the coefficient correlation of model two is higher than model one, and model two is Best model, from model two.
2) natural gamma GR, the lithology and density logging DEN of S wells to be explained are obtained by well-log information, passes through gas measuring logging data JingSJing to be explained total hydrocarbon, methane content is obtained, landwaste, the rock core of geological logging to be explained are obtained by geological logging data Lithology;
3) lithology of the natural gamma of the S wells obtained by step 2), lithology and density logging, total hydrocarbon, methane and landwaste rock core is drawn Divide S well shale gas reservoirs, S wells shale gas mainly shows interval in 2170.0~2220.0m;
4) with step 2) using best model Sg=-262.896 ρ b+0.183GR+768.147 and S wells density and Natural gamma material computation S well reservoir gas-bearing saturation degrees Sg;
5) export result of calculation, S 2170.0~2220.0m of well well sections calculate gas saturation substantial scope be 50.0%~ 80.0% (see Fig. 4), contrasted with the well rock core saturation degree test result, mean error is less than 10%, disclosure satisfy that shale gas stores up Layer and shale gas Gas reservoir eraluation need.

Claims (4)

1. the method for shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density, it is characterised in that concretely comprise the following steps:
1) work area fixed well core sample test analysis gas saturation Sg, rock core bulk density ρ b are obtained and corresponds to depth Log well natural gamma GR;The sample of selection must be same work area, identical series of strata, and quantity can not be less than 12 pieces, and well number is unlimited; Then Sg1=a/ ρ b+b of model one, the Sg2=a ρ b+bGR+c of model two model coefficient are determined using least square method, most Best model is determined eventually;
Rock core bulk density ρ b dimensions are g/cm in formula3,
Natural gamma GR dimensions are API,
Gas saturation Sg dimensions are %;
2) natural gamma GR, the lithology and density logging DEN of well to be explained are obtained by well-log information, is obtained and treated by gas measuring logging data Total hydrocarbon, the methane content of well are explained, landwaste, the rock core lithology of geological logging to be explained are obtained by geological logging data;Lithology Density DEN dimensions are g/cm3
3) well logging obtained by step 2) and log data divide shale gas reservoir;
4) step 2) borehole logging tool lithology and density logging DEN alternative steps 1 to be explained are used) rock core bulk density ρ b, pass through the optimal of determination Model and well-log information calculate well reservoir gas-bearing saturation degree Sg to be explained;
5) result of calculation is exported.
2. the method according to claim 1 that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density, its feature Be, returned using the gas saturation that model inverse obtains and core test saturation degree in step 1), coefficient correlation it is high be Best model.
3. the method according to claim 1 that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density, its feature It is, the natural gamma GR of well to be explained, lithology and density logging DEN data depths are at intervals of 0.1m~1.0m in step 2);Gas surveys record Well total hydrocarbon and methane content depth interval are 1m or 0.5m.
4. the method according to claim 1 that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density, its feature It is, step 3) shale gas reservoir division foundation:Lithology is shale or argiloid;Be present obvious exception in total hydrocarbon, methane, be free of Heavy hydrocarbon composition;High natural gamma, low-density.
CN201510753459.3A 2015-11-06 2015-11-06 The method that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density Active CN105350959B (en)

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 CN105350959A (en) 2016-02-24
CN105350959B true 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)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108561126B (en) * 2018-03-13 2021-09-14 中石化石油工程技术服务有限公司 Simple method for determining organic porosity of shale gas reservoir
CN108647417A (en) * 2018-04-28 2018-10-12 中石化石油工程技术服务有限公司 A kind of simple and easy method of determining shale gas reservoir gas-bearing saturation degree

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2385413C1 (en) * 2008-09-30 2010-03-27 Шлюмберже Текнолоджи Б.В. Method of evaluating current gas saturation in bottomhole zone of well in deposit of volatile oil
CN101487390B (en) * 2009-02-23 2012-08-29 大庆油田有限责任公司 Archie mode method for confirming initial oil saturation of oil layer
CN102748004B (en) * 2012-06-06 2015-08-05 中国石油化工股份有限公司 The defining method of mud logging hydrocarbon saturation
CN103912269B (en) * 2014-03-28 2017-04-12 中石化江汉石油工程有限公司测录井公司 Method for determining formation fracture pressure gradient logging of shale gas reservoir
CN103912268B (en) * 2014-03-28 2016-09-28 中石化江汉石油工程有限公司测录井公司 Shale reservoir gas-bearing saturation based on TOC determines method

Also Published As

Publication number Publication date
CN105350959A (en) 2016-02-24

Similar Documents

Publication Publication Date Title
Slatt et al. Merging sequence stratigraphy and geomechanics for unconventional gas shales
CN104453873B (en) Shale oil-gas economical and effective layer section evaluation method
CN101930082B (en) Method for distinguishing reservoir fluid type by adopting resistivity data
CN104514552B (en) A kind of method that coalbed methane reservoir identification is predicted with abundance
Cudjoe et al. A comprehensive approach to sweet-spot mapping for hydraulic fracturing and CO2 huff-n-puff injection in Chattanooga shale formation
CN103513278B (en) Seismic wave groups thickness is utilized to carry out the method for reservoir prediction
Bhattacharya et al. Flow unit modeling and fine-scale predicted permeability validation in Atokan sandstones: Norcan East field, Kansas
Sanderson et al. Quantitative analysis of tin-and tungsten-bearing sheeted vein systems
CN103821505A (en) Sandstone oil and gas carrier bed geophysics-geology-geochemistry detection method and device
CN105275456A (en) Method for identifying high-quality shale by using logging information
CN105350959B (en) The method that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density
CN105488248A (en) Judgment method for deep ore stratum and rock stratum
CN104991286A (en) Sedimentary facies characterization method based on sedimentary modes
Li et al. Fracture permeability evaluation of a coal reservoir using geophysical logging: a case study in the Zhengzhuang area, southern Qinshui Basin
Vinassa et al. A comprehensive approach to sweet-spot mapping for hydraulic fracturing and CO2 huff-n-puff injection in Chattanooga shale formation
Wang et al. Improved methods for determining effective sandstone reservoirs and evaluating hydrocarbon enrichment in petroliferous basins
Li et al. Exploration methods for late Quaternary shallow biogenic gas reservoirs in the Hangzhou Bay area, eastern China
Campo et al. Late Quaternary sequence stratigraphy as a tool for groundwater exploration: Lessons from the Po River Basin (northern Italy)
Chatterjee et al. Relation between coalbed permeability and in-situ stress magnitude for coalbed methane exploration in Jharia and Raniganj coalfields, India
Lemieux et al. Fracture network analysis of the St-Eustache quarry, Quebec, Canada, for groundwater resources management
CN106285623A (en) Determine the method and system of total content of organic carbon
CN103278852B (en) Utilize the method for seismic data volume waveform configuration characteristic model predicting oil/gas
Xiao et al. Log response of ultrasonic imaging and its significance for deep mineral prospecting of scientific drilling borehole-2 in Nanling district, China
Merletti et al. Integration of depositional, petrophysical, and petrographic facies for predicting permeability in tight gas reservoirs
CN107657365A (en) Geological resource extraction value evaluation method and device

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into 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