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