CN105370270B - The method that shale gas reservoir gas-bearing saturation degree is determined by the dipole sonic P-wave And S time difference - Google Patents

The method that shale gas reservoir gas-bearing saturation degree is determined by the dipole sonic P-wave And S time difference Download PDF

Info

Publication number
CN105370270B
CN105370270B CN201510752586.1A CN201510752586A CN105370270B CN 105370270 B CN105370270 B CN 105370270B CN 201510752586 A CN201510752586 A CN 201510752586A CN 105370270 B CN105370270 B CN 105370270B
Authority
CN
China
Prior art keywords
well
gas
time difference
shale gas
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
CN201510752586.1A
Other languages
Chinese (zh)
Other versions
CN105370270A (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 CN201510752586.1A priority Critical patent/CN105370270B/en
Publication of CN105370270A publication Critical patent/CN105370270A/en
Application granted granted Critical
Publication of CN105370270B publication Critical patent/CN105370270B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The present invention relates to a kind of method for determining shale gas reservoir gas-bearing saturation degree by the dipole sonic P-wave And S time difference, analyzed by core test and obtain shale gas saturation Sg, the dipole sonic that corresponding depth is obtained by same well dipole well-log information indulges shear wave slowness, calculates transverse and longitudinal time difference ratio;It is fitted with the well logging natural gamma of corresponding depth, transverse and longitudinal time difference ratio and shale gas reservoir core test gas saturation, establishes model, determine model coefficient;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 gas detection logging to be explained are obtained by gas measuring logging data, landwaste, the rock core lithology of well to be explained are obtained by geological logging data;Divide shale gas reservoir, the gas saturation of the shale gas reservoir of computation partition;Export result of calculation.It is of the invention easy, it is applied widely.The present invention applies 80 mouthfuls of wells in certain shale gas field, and the mean error for calculating determination is no more than 10%, can meet needs.

Description

The method that shale gas reservoir gas-bearing saturation degree is determined by the dipole sonic P-wave And S time difference
Technical field
The present invention relates to a kind of method for determining shale gas reservoir gas-bearing saturation degree by the dipole sonic P-wave And S time difference.
Background technology
Shale gas gas saturation is the indispensable important parameter that shale gas reservoir surveys well logging evaluation.Traditional Dipole sonic shale gas reservoir gas-bearing saturation degree determines that method is asked for based on elastic modulus method, is only applicable to sandstone and carbonate Reservoir, more desired zone formation parameter, method complexity, narrow application range are calculated, limitation is larger in shale gas reservoir, substantially It is inapplicable.
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 desirable to provide a kind of easy, applied widely, identified page The shale gas-bearing formation gas saturation that rock gas reservoir gas saturation obtains with core test approaches, and can meet that well logging reservoir is surveyed at scene Evaluation and the side that shale gas reservoir gas-bearing saturation degree is determined by the dipole sonic P-wave And S time difference for calculating free Gas content needs Method.
The implementation of the object of the invention is to determine shale gas reservoir gas-bearing saturation degree by the dipole sonic P-wave And S time difference Method, concretely comprise the following steps:
1) shale gas saturation Sg is obtained by work area fixed well core sample test analysis, passes through dipole well-log information The dipole sonic shear wave slowness DTS and dipole sonic compressional wave time difference DTC of corresponding depth are obtained, and calculate the transverse and longitudinal time difference to compare RMSC =DTS/DTC;
In formula:Dipole sonic shear wave slowness DTS dimensions are μ s/ft,
Dipole sonic compressional wave time difference DTC dimensions are μ s/ft,
The transverse and longitudinal time difference is than RMSC dimensionless;
2) the shale gas saturation Sg that core test analysis is obtained by Conventional Logs corresponds to the well logging of depth certainly Right gamma GR;
3) shale gas reservoir gas-bearing saturation degree model Sg=aRMSC+bGR+c is used, carries out the page of core test acquisition Rock gas-bearing formation gas saturation Sg is fitted with the well logging transverse and longitudinal time difference than RMSC, natural gamma GR, and determines model coefficient, natural in formula Gamma GR dimensions are API;
4) natural gamma, the lithology and density logging 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 gas detection logging are explained, landwaste, the rock core lithology of well to be explained are obtained by geological logging data;
5) well logging obtained by step 4) and log data divide shale gas reservoir;
6) shale gas reservoir gas-bearing is calculated according to the gas saturation model of step 3), the shale gas reservoir of step 5) division Saturation degree (Sg);
7) result of calculation is exported.
The present invention solves conventional dipole acoustic wave methodogy and is calculating shale gas reservoir gas-bearing saturation degree desired zone stratum ginseng More, the narrow application ranges of number, the problem of limitation is big, shale reservoir gas saturation are determined compared with the dipole sonic based on modulus of elasticity Method is easier, and the scope of application is wider;Dipole acoustic log logging data application is further expanded.
The present invention applies 80 mouthfuls of wells in certain shale gas field, and the shale gas reservoir gas-bearing saturation degree and rock core for calculating determination are surveyed Examination obtain shale gas-bearing formation gas saturation approach, mean error be no more than 10%, can meet scene survey well logging evaluating reservoir with Calculate the needs of free Gas content.
Brief description of the drawings
Fig. 1 is workflow block diagram of the present invention,
Fig. 2 is A4 well exemplary application maps in R work areas of the present invention,
Fig. 3 is A1-3HF well exemplary application maps in R work areas of the present invention.
Embodiment
The applicant by shale gas reservoir core experimental studies have found that, in the case of shale lithofacies assemblages are single, shale Gas reservoir gas saturation (Sg) and the quantitative relationship of the horizontal compressional wave time difference ratio (RMSC) of dipole acoustic log are Sg=ARMSCB, A, B is coefficient.Shale gas Reservoir Lithofacies such as certain shale gas field R work area are combined as the A=of the single petrofacies well section of siliceous shale 1950.3, B=-6.606;Shale gas Reservoir Lithofacies are combined as A=162.9, B=- of the single petrofacies well section of culm 1.398.When shale reservoir lithology is two kinds and two or more lithofacies assemblages, by shale gas reservoir gas-bearing saturation degree model Sg= A·RMSCBThe shale gas reservoir gas-bearing saturation degree Sg and horizontal compressional wave time difference ratio RMSC of dipole sonic is fitted, coefficient correlation is relatively low, and And the evaluation of the shale gas reservoir gas-bearing saturation degree of actual field is the evaluation based on whole section of interval, in the absence of only single petrofacies Situation, therefore requirement of engineering can not be met, be suitable only for scientific research needs.Shale gas reservoir gas-bearing saturation degree (Sg) and dipole sonic Log well horizontal compressional wave time difference ratio (RMSC), natural gamma (GR) by model Sg=aRMSC+bGR+c be fitted, coefficient correlation compared with Height, the Sg being calculated can meet that well logging shale gas evaluating reservoir needs are surveyed at scene.
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) shale gas saturation Sg is obtained by work area fixed well core sample test analysis, passes through dipole well-log information The dipole sonic shear wave slowness DTS and dipole sonic compressional wave time difference DTC of corresponding depth are obtained, and calculate the transverse and longitudinal time difference to compare RMSC =DTS/DTC;
In formula:Dipole sonic shear wave slowness DTS dimensions are μ s/ft,
Dipole sonic compressional wave time difference DTC dimensions are μ s/ft,
The transverse and longitudinal time difference is than RMSC dimensionless;
2) the shale gas saturation Sg that core test analysis is obtained by Conventional Logs corresponds to the well logging of depth certainly Right gamma GR;
3) shale gas reservoir gas-bearing saturation degree model Sg=aRMSC+bGR+c is used, carries out the page of core test acquisition Rock gas-bearing formation gas saturation Sg, than RMSC, natural gamma GR least square fittings, determines model coefficient with the well logging transverse and longitudinal time difference A, b, c, natural gamma GR dimensions are API in formula;
4) natural gamma, the lithology and density logging 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 gas detection logging are explained, landwaste, the rock core lithology of well to be explained are obtained by geological logging data;
The natural gamma of well to be explained, lithology and density logging data depth save as wis or txt lattice at intervals of 0.1m~1.0m Formula;Gas detection logging total hydrocarbon and methane content depth interval are 1m or 0.5m, save as wis or txt forms;Geologic debris, rock core Lithology data save as excel or txt forms.
5) well logging obtained by step 4) and log data divide shale gas reservoir.
Reservoir division foundation:Lithology is argiloid or shale;There is obvious exception in total hydrocarbon, methane, without heavy hydrocarbon composition; High natural gamma, low-density.
6) shale gas reservoir gas-bearing is calculated according to the gas saturation model of step 3), the shale gas reservoir of step 5) division Saturation degree (Sg);
7) 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 specific embodiment below.
Example one:Certain shale gas field R work areas A4 well
1) gassiness of 4 mouthfuls of shale gas wells totally 93 pieces of shale cores is obtained by the analysis of certain shale gas field R work area core test Saturation degree (Sg) and the shear wave slowness (DTS) and compressional wave time difference (DTC) that corresponding depth is obtained by dipole well-log information, and The calculating transverse and longitudinal time difference compares RMSC=DTS/DTC;
In formula:Dipole sonic shear wave slowness DTS dimensions are μ s/ft,
Dipole sonic compressional wave time difference DTC dimensions are μ s/ft,
The transverse and longitudinal time difference is than RMSC dimensionless;
2) natural gamma (GR) of depth is corresponded to by Conventional Logs 93 pieces of shale cores of acquisition;
3) gas saturation (Sg) of 93 blocks of rock cores and well logging the transverse and longitudinal time difference ratio RMSC and natural gamma GR are pressed into model Sg =aRMSC+bGR+c least square methods calculate, and it is a=-159.123, b=-0.037, c=to determine model coefficient 335.716;
4) natural gamma, the lithology and density logging data of A4 wells are obtained by well-log information, A4 is obtained by gas measuring logging data Total hydrocarbon, the methane content of well, landwaste, the rock core lithology of A4 wells are obtained by geological logging data;
The natural gamma of well to be explained, lithology and density logging data depth save as wis or txt lattice at intervals of 0.1m~1.0m Formula;Gas detection logging total hydrocarbon and methane content depth interval are 1m or 0.5m, save as wis or txt forms;Geologic debris, rock core Lithology data save as excel or txt forms.
5) by obtaining the natural gammas of A4 wells, lithology and density logging, total hydrocarbon, methane and the lithology classification A4 well pages of landwaste rock core Rock gas reservoir, A4 wells shale gas mainly show interval in 2337.0~2398.0m;
6) according to the gas saturation model Sg=-159.123RMSC-0.037GR+335.716 of step 3), step 5) the shale gas reservoir of division calculates A4 well shale gas reservoir gas-bearing saturation degrees (Sg);
7) result of calculation is exported, it is 50.0% that A4 2337.0~2398.0m of well well sections, which calculate gas saturation substantial scope, ~80.0% (see Fig. 2), contrasted with the well rock core saturation degree test result, error is less than 10%, disclosure satisfy that shale gas reservoir Needed with shale gas Gas reservoir eraluation.
Example two:Certain shale gas field R work areas horizontal well A1-3HF well
1) the A4 wells of A1-3HF wells and example one are all the wells in R work areas, therefore model coefficient is identical, can continue to use reality The model coefficient of example one, model coefficient a=-159.123, b=-0.037, c=335.716;
2) natural gamma, the lithology and density logging of A1-3HF wells to be explained are obtained by well-log information, passes through gas measuring logging data Total hydrocarbon, the methane content of A1-3HF wells to be explained are obtained, landwaste, the rock of A1-3HF wells to be explained are obtained by geological logging data Heart lithology;
The natural gamma of well to be explained, lithology and density logging data depth save as wis or txt lattice at intervals of 0.1m~1.0m Formula;Gas detection logging total hydrocarbon and methane content depth interval are 1m or 0.5m, save as wis or txt forms;Geologic debris, rock core Lithology data save as excel or txt forms.
3) well logging obtained by step 2) and log data divide A1-3HF well shale gas reservoirs, A1-3HF well shale gas Main display interval is in 2710.0~2775.0m;
4) according to the gas saturation model Sg=-159.123RMSC-0.037GR+335.716 of step 1), step 3) the shale gas reservoir of division calculates A1-3HF well shale gas reservoir gas-bearing saturation degrees (Sg);
5) export result of calculation, A1-3HF 2710.0~2775.0m of well well sections calculate gas saturation be 55.0%~ 70.0% (see Fig. 3).It disclosure satisfy that shale gas reservoir and shale gas Gas reservoir eraluation needs.

Claims (4)

1. the method that shale gas reservoir gas-bearing saturation degree is determined by the dipole sonic P-wave And S time difference, it is characterised in that:Specific steps For:
1) shale gas reservoir gas-bearing saturation degree Sg is obtained by work area fixed well core sample test analysis, is logged well and provided by dipole Material obtains the dipole sonic shear wave slowness DTS and dipole sonic compressional wave time difference DTC of corresponding depth, and calculates transverse and longitudinal time difference ratio RMSC=DTS/DTC;
In formula:Dipole sonic shear wave slowness DTS dimensions are μ s/ft,
Dipole sonic compressional wave time difference DTC dimensions are μ s/ft,
The transverse and longitudinal time difference is than RMSC dimensionless;
2) the shale gas reservoir gas-bearing saturation degree Sg analyzed by Conventional Logs acquisition core test corresponds to the well logging of depth Natural gamma GR;
3) shale gas reservoir gas-bearing saturation degree model Sg=aRMSC+bGR+c is used, carries out the shale gas of core test acquisition Reservoir gas-bearing saturation degree Sg is fitted with the well logging transverse and longitudinal time difference than RMSC, natural gamma GR, and determines model coefficient, natural gal in formula Horse GR dimensions are API;
4) natural gamma, the lithology and density logging of well to be explained are obtained by well-log information, is obtained by gas measuring logging data to be explained Total hydrocarbon, the methane content of well gas detection logging, landwaste, the rock core lithology of well to be explained are obtained by geological logging data;
5) well logging obtained by step 4) and log data divide shale gas reservoir;
6) shale gas reservoir gas-bearing saturation is calculated according to the gas saturation model of step 3), the shale gas reservoir of step 5) division Spend Sg;
7) result of calculation is exported.
2. the method according to claim 1 that shale gas reservoir gas-bearing saturation degree is determined by the dipole sonic P-wave And S time difference, It is characterized in that:Step 3) carries out least square fitting with model Sg=aRMSC+bGR+c, determine model coefficient a, b, c。
3. the method according to claim 1 that shale gas reservoir gas-bearing saturation degree is determined by the dipole sonic P-wave And S time difference, It is characterized in that:The natural gamma of step 4) well to be explained, lithology and density logging data depth are at intervals of 0.1m~1.0m;Gas detection logging 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 the dipole sonic P-wave And S time difference, It is characterized in that:Step 5) reservoir division foundation:Lithology is argiloid or shale;Be present obvious exception in total hydrocarbon, methane, be free of Heavy hydrocarbon composition;High natural gamma, low-density.
CN201510752586.1A 2015-11-06 2015-11-06 The method that shale gas reservoir gas-bearing saturation degree is determined by the dipole sonic P-wave And S time difference Active CN105370270B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510752586.1A CN105370270B (en) 2015-11-06 2015-11-06 The method that shale gas reservoir gas-bearing saturation degree is determined by the dipole sonic P-wave And S time difference

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510752586.1A CN105370270B (en) 2015-11-06 2015-11-06 The method that shale gas reservoir gas-bearing saturation degree is determined by the dipole sonic P-wave And S time difference

Publications (2)

Publication Number Publication Date
CN105370270A CN105370270A (en) 2016-03-02
CN105370270B true CN105370270B (en) 2018-04-10

Family

ID=55372828

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510752586.1A Active CN105370270B (en) 2015-11-06 2015-11-06 The method that shale gas reservoir gas-bearing saturation degree is determined by the dipole sonic P-wave And S time difference

Country Status (1)

Country Link
CN (1) CN105370270B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106842303B (en) * 2016-12-30 2019-07-05 中国石油天然气股份有限公司 Method and device for determining gas saturation
CN107422384B (en) * 2017-03-20 2019-03-01 中石化石油工程技术服务有限公司 A method of vertical seismic profile data is imitated using well-log information
CN108733856B (en) * 2017-04-21 2022-06-21 中国石油化工股份有限公司 Shale gas reservoir free gas saturation determination method and computer readable storage medium
CN114060015B (en) * 2020-07-31 2024-05-03 中国石油化工股份有限公司 Method and device for evaluating gas content of compact sandstone
CN112065361B (en) * 2020-09-16 2021-03-12 成都理工大学 Method for determining gas saturation of tight reservoir based on sound wave attenuation

Family Cites Families (6)

* 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
CN102454401B (en) * 2010-10-29 2015-05-20 中国石油化工股份有限公司 Method for obtaining logging saturation of low porosity permeability reservoir
CN102748004B (en) * 2012-06-06 2015-08-05 中国石油化工股份有限公司 The defining method of mud logging hydrocarbon saturation
CN102913238A (en) * 2012-10-29 2013-02-06 陕西联盟物流有限公司 Method for dissolving initial oil saturation degree of oil-water layer
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
CN105370270A (en) 2016-03-02

Similar Documents

Publication Publication Date Title
CN102640018B (en) Workflow for petrophysical and geophysical formation evaluation of wireline and LWD log data
CN105370270B (en) The method that shale gas reservoir gas-bearing saturation degree is determined by the dipole sonic P-wave And S time difference
Perez et al. Calibration of brittleness to elastic rock properties via mineralogy logs in unconventional reservoirs
Bodziak et al. The role of seismic attributes in understanding the hydraulically fracturable limits and reservoir performance in shale reservoirs: An example from the Eagle Ford Shale, south Texas
US9229910B2 (en) Predicting three dimensional distribution of reservoir production capacity
Aderibigbe et al. Integrated rock classification in the Wolfcamp Shale based on reservoir quality and anisotropic stress profile estimated from well logs
Hamza et al. Characterization of anisotropic elastic moduli and stress for unconventional reservoirs using laboratory static and dynamic geomechanical data
Yasin et al. Study on brittleness templates for shale gas reservoirs-A case study of Longmaxi shale in Sichuan Basin, southern China
Manriquez et al. A novel approach to quantify reservoir pressure along the horizontal section and to optimize multistage treatments and spacing between hydraulic fractures
Wojtowicz et al. Reconstructing the mechanical parameters of a transversely-isotropic rock based on log and incomplete core data integration
Gallardo-Giozza et al. Effects of elastic anisotropic models on the prediction of horizontal stresses and hydraulic fracture geometry
Zorn et al. Geomechanics of the Microseismic Response in Devonian Organic Shales at the Marcellus Shale Energy and Environment Laboratory (MSEEL) Site, West Virginia
He et al. Rock skeleton models and seismic porosity inversion
Bratton Stress induced changes in elastic wave attributes in the Wattenberg Field, Colorado, USA
Rasouli Geomechanics of gas shales
RU2585296C1 (en) Method of determining drained hydraulic fracturing crack width and degree of sedimentation of proppant therein
Yu et al. Rock physics diagnostics and modeling for shale gas formation characterization in China
Altowairqi et al. Measuring ultrasonic characterisation to determine the impact of TOC and the stress field on shale gas anisotropy
Tiancheng et al. Source mechanism studies of acoustic emission in large-scale hydraulic fracturing experiment
Tao et al. Application of multipole array sonic logging to acid hydralic fracturing
Sharma et al. Replacing conventional brittleness indices determination with new attributes using true hydrofracturing mechanism
Medina et al. 3D Geomechanical Modeling for Field Development of a Colombian Shale Play
Soltanzadeh et al. Application of an integrated approach for the characterization of mechanical rock properties in the Duvernay formation
Porter et al. A statistical method for determination of water saturation from logs
Masri Amplitude Versus Offset analysis–A possible useful tool for geothermal exploration

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
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.