CN102183788A - Method for synthesizing transverse-wave well-logging curve under stratum condition - Google Patents

Method for synthesizing transverse-wave well-logging curve under stratum condition Download PDF

Info

Publication number
CN102183788A
CN102183788A CN 201110124614 CN201110124614A CN102183788A CN 102183788 A CN102183788 A CN 102183788A CN 201110124614 CN201110124614 CN 201110124614 CN 201110124614 A CN201110124614 A CN 201110124614A CN 102183788 A CN102183788 A CN 102183788A
Authority
CN
China
Prior art keywords
mineral
swal
well
stratum
transverse
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
Application number
CN 201110124614
Other languages
Chinese (zh)
Other versions
CN102183788B (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.)
SHANGHAI PETROLEUM CO Ltd
Original Assignee
SHANGHAI PETROLEUM 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 SHANGHAI PETROLEUM CO Ltd filed Critical SHANGHAI PETROLEUM CO Ltd
Priority to CN2011101246147A priority Critical patent/CN102183788B/en
Publication of CN102183788A publication Critical patent/CN102183788A/en
Application granted granted Critical
Publication of CN102183788B publication Critical patent/CN102183788B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses a method for synthesizing a transverse-wave well-logging curve under a stratum condition. The method comprises the following steps of: standardizing the well-logging curve; selecting a key well with transverse-wave well-logging data in a research area; selecting a marker bed in the key well; computing mineral constituents in the stratum; according to a geological layering condition in the area, on the basis of main mineral constituents in a rock core analysis, establishing a multi-mineral model of each bed; computing transverse-wave well-logging response of each mineral constituent of the stratum rock; re-establishing the transverse-wave well-logging curve; and fitting according to the transverse-wave well-logging responses of the mineral components and the computed volumes of the mineral constituents in each bed to obtain the transverse-wave well-logging curve of the whole stratum. By using the conventional well-logging curve and the fine multi-mineral model and according to the transverse-wave well-logging responses of the constituents of different minerals, pores, fluids and the like, a transverse-wave curve representing the features of the stratum is fitted out and the fitting precision and accuracy are very high.

Description

The synthetic method of SWAL curve under the formation condition
Technical field
The present invention relates to a kind of synthetic method of logging trace, particularly relate to the synthetic method of SWAL curve under a kind of formation condition.
Background technology
The SWAL information material to earthquake instead shake, the resident fluid detection and Identification have important effect.The SWAL information material does not all have measured because cost is higher in most prospecting boreholes in early days and the present exploitation well.For the oil field of exploitation middle and later periods, a new round to tap the latent power, it is essential that seismic inversion and resident fluid detect.Therefore, it is extremely important to carry out shear wave reconstruct (synthesizing) for the well that does not have the SWAL data.The software of shear wave reconstruct aspect is many at present, mainly contains slog, Xu-White theoretical model etc., and is expensive but software is bought, and implements also more complicated; Also have simple volume-based model and p-and s-wave velocity ratio method etc., but the shear wave ratio of precision of reconstruct is relatively poor.Naive model method commonly used just is divided into the stratum rock skeleton and hole two parts, and to carry out shear wave synthetic, and the value of the SWAL of these two parts response and operating personnel's experience have much relations, cause error easily.
Summary of the invention
Technical matters to be solved by this invention provides the synthetic method of SWAL curve under a kind of formation condition, it utilizes the conventional logging curve, adopt meticulous multimineral model, SWAL response according to components such as different minerals, hole, fluids, simulate the shear wave curve of a reflection stratum characteristic, accuracy and precision are all very high.
The present invention solves above-mentioned technical matters by following technical proposals: the synthetic method of SWAL curve under a kind of formation condition is characterized in that it may further comprise the steps:
Step 1, logging trace standardization: choose the key well that the interior a bite of survey region has the SWAL data, in key well, choose reference lamina, well-log information in the zone is carried out standardization according to the log response feature of key well reference lamina;
The calculating of step 2, stratum mineral constituent: according to the geology layering situation in the zone,, set up the multimineral model of each layer, calculate the mineral constituent volume on stratum based on the essential mineral component of core analysis;
The calculating of step 3, each mineral constituent SWAL response of formation rock: according to response of key well SWAL and stratum mineral constituent content, Inversion Calculation obtains each mineral constituent SWAL response of formation rock;
The reconstruct of step 4, SWAL curve: utilize the SWAL response of these mineral constituents and the mineral constituent volume of each layer of having calculated just can match obtains the SWAL curve on whole stratum.
Positive progressive effect of the present invention is: the SWAL curve of institute of the present invention reconstruct and the shear wave curve of actual measurement are overlapping substantially, and accuracy and precision are all very high.
Description of drawings
Fig. 1 is the process flow diagram of the synthetic method of SWAL curve under the formation condition of the present invention.
Fig. 2 is the stratigraphic model synoptic diagram of key well among the present invention.
Fig. 3 is for needing the stratigraphic model synoptic diagram of synthetic shear wave well among the present invention.
Embodiment
Provide preferred embodiment of the present invention below in conjunction with accompanying drawing, to describe technical scheme of the present invention in detail.
As shown in Figure 1, the synthetic method of SWAL curve may further comprise the steps under the formation condition of the present invention:
Step 1, logging trace standardization: choose the key well that the interior a bite of survey region has the SWAL data, the well-log information of key well is wanted complete and quality better.Choose reference lamina in key well, reference lamina generally will have certain thickness, distributional stability in the whole district, and the log response feature is obvious.According to the log response feature of type formation other conventional logging data of the well that do not have the shear wave data in the zone is carried out standardization, promptly well-log information in the zone is carried out standardization according to the log response feature of key well reference lamina.
The calculating of step 2, stratum mineral constituent: according to the geology layering situation in the zone,, set up the multimineral model of each layer, calculate the mineral constituent volume on stratum based on the essential mineral component of core analysis.The multimineral model method is all to regard various mineralogical compositions in the stratum and pore fluid as different mineral volume, and all mineral volume content sums are 1 in the stratum.The stratum can be represented by volume-based model for the response equation of logging instrumentation.Be provided with N-1 bar logging trace, calculate M mineral content (comprising volume of voids), and N＞=M, add that balance equation can list N equation.So just formed linear overdetermined equation group.The uranium (U), thorium (TH), the potassium logging traces such as (K) that for example adopt density (DEN), natural gamma (GR), sound wave (AC), neutron (CNL) and natural gamma spectra to obtain can obtain following log response equation group:
DEN＝ρ 1V 1+ρ 2V 2+…+ρ iV i+ρ φφ
GR＝gr 1V 1+gr 2V 2+…+gr iV i+gr φφ
AC＝ac 1V 1+ac 2V 2+…ac iV i+ac φφ
CNL＝cnl 1V 1+cnl 2V 2+…cnl iV i+cnl φφ
U＝n 1V 1+n 2V 2+…+n iV i+n φφ
ρ in the formula (1) i, gr i, ac i, cnl i, u iLog response Deng the expression mineral; V iExpression different minerals component volume content; ρ φ, gr φ, ac φ, cnl φ, u φLog response (or equivalent log response) Deng the expression pore fluid; φ is the factor of porosity on stratum.M mineral volume is volume of voids V fN equation is material balance equation.By finding the solution above-mentioned equation, just can obtain the volume of each mineral constituent of formation rock and the factor of porosity on stratum.
The calculating of step 3, each mineral constituent SWAL response of formation rock: the volume and the formation porosity that obtain each mineral constituent of formation rock by step 2.Set up equation according to response of key well SWAL and stratum mineral constituent content, Inversion Calculation obtains the SWAL response of each mineral constituent of stratum and the equivalent SWAL response of formation pore fluid, and concrete formula is as follows:
DTSM 1＝dtsm 1V 11+dtsm 2V 22+...+dtsm iV n+dtsm φφ j
DTSM 2＝dtsm 1V 12+dtsm 2V 22+...+dtsm iV n+dtsm φφ j
.
.
.
.
DTSM j＝dtsm 1V 1j++dtsm 2V 2j+…+dtsm iV n+dtsm φφ j
1＝V 1j+V 2j+...V n+φ j
…………(2)
In the formula (2), j represents the sampled point on each degree of depth, DTSM jRepresent the SWAL response of j depth-sampling point; I represents different mineral numberings, dtsm iRepresent the SWAL response of i kind mineral; Dtsm φRepresent equivalent pore fluid SWAL response, φ jThe factor of porosity of representing j depth point.
In some geology layerings, just can form a linear overdetermined equation group like this according to many mineral volume-based model.Utilize least square method to find the solution this system of equations, just can obtain the SWAL response dtsm of stratum different minerals 1, dtsm 2Dtsm i, dtsm φ
The reconstruct of step 4, SWAL curve: utilize the SWAL response of these mineral constituents and the mineral constituent volume of each layer of having calculated just can match obtains the SWAL curve on whole stratum.He Cheng shear wave curve precision height not only like this, and do not need too many loaded down with trivial details parameter.In the well that does not have the SWAL curve, in a certain geology layering, adopt the method for step 2 to calculate the mineral constituent and the factor of porosity V on stratum 1, V 2... V i, φ; The SWAL response of the different minerals component in the zone of calculating according to step 3 in the same geology layering and the equivalent SWAL response of pore fluid, according to many mineral volume-based model principle, the shear wave response on stratum equals the summation of each mineral constituent shear wave response, promptly following formula:
dtsm 1V 1+dtsm 2V 2+…+dtsm kV k+dtsm φφ＝DTSMc?......(3)
In the formula (3), V 1, V 2... V i, φ represents a sets of curves respectively, DTSMc is exactly this interval shear wave slowness of reconstruct.The SWAL curve of institute of the present invention reconstruct and the shear wave curve of actual measurement are overlapping substantially, and accuracy and precision are all very high.
As shown in Figures 2 and 3, lift a concrete example, wherein 1., 2., 3., 4. numbering represents four kinds of different mineral, the volume on this shared stratum of mineral of Regional Representative at each numbering place, and the synthetic method of SWAL curve realizes as follows under the formation condition of the present invention:
Step 1, the essential mineral component of regional rock can be according to the core analysis gained.Through after the standardization, utilize the system of linear equations of conventional logging curve foundation, just can ask for 1., 2., 3., 4. wait the percentage of the shared stratum volume of different mineral.
Step 2 is calculated the mineral volume content of each well geology layering ground floor.First mineral volume content 1. that for example calculates key well is 79.2%, second mineral volume content 2. is that the 15.4%, the 3rd mineral volume content 3. is that the 3.2%, the 4th mineral volume content 4. is 2.2% (factor of porosity, factor of porosity is represented in general last a kind of mineral constituent).
Step 3 is 163.5us/ft (us/ft is the record unit of logging trace) according to known key well at the whole district corresponding DTSM of geological layer ground floor SWAL; By the least square method inverting just can calculate 1., 2., 3., 4. the SWAL individual features of four kinds of mineral is respectively 178.5us/ft, 73us/ft, 91us/ft, 360us/ft (equivalent SWAL is corresponding).
Step 4, calculating synthetic geological layer first mineral of the shear wave well volume content 1. of needs according to step 1, step 2 is that 66.8%, second mineral volume content 2. is that the 22.7%, the 3rd mineral volume content 3. is that the 10%, the 4th mineral volume content 4. is 0.5% (factor of porosity), and the SWAL response of each mineral constituent of step 3 calculating ground floor, the SWAL that just can synthesize ground floor is 178.5*66.8%+73*22.7%+91*10%+360*0.5%=146.7us/ft.
In like manner, just can calculate the SWAL curve of whole well sections such as needing the synthetic shear wave well second layer, the 3rd layer.Then, can synthesize the SWAL curve that obtains other well.
Though more than described the specific embodiment of the present invention, it will be understood by those of skill in the art that these only illustrate, under the prerequisite that does not deviate from principle of the present invention and essence, can make numerous variations or modification to these embodiments.Therefore, protection scope of the present invention is limited by appended claims.

Claims (1)

1. the synthetic method of SWAL curve under the formation condition is characterized in that it may further comprise the steps:
Step 1, logging trace standardization: choose the key well that the interior a bite of survey region has the SWAL data, in key well, choose reference lamina, well-log information in the zone is carried out standardization according to the log response feature of key well reference lamina;
The calculating of step 2, stratum mineral constituent: according to the geology layering situation in the zone,, set up the multimineral model of each layer, calculate the mineral constituent volume on stratum based on the essential mineral component of core analysis;
The calculating of step 3, each mineral constituent SWAL response of formation rock: according to response of key well SWAL and stratum mineral constituent content, Inversion Calculation obtains each mineral constituent SWAL response of formation rock;
The reconstruct of step 4, SWAL curve: utilize the SWAL response of these mineral constituents and the mineral constituent volume of each layer of having calculated just can match obtains the SWAL curve on whole stratum.
CN2011101246147A 2011-05-13 2011-05-13 Method for synthesizing transverse-wave well-logging curve under stratum condition Active CN102183788B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011101246147A CN102183788B (en) 2011-05-13 2011-05-13 Method for synthesizing transverse-wave well-logging curve under stratum condition

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011101246147A CN102183788B (en) 2011-05-13 2011-05-13 Method for synthesizing transverse-wave well-logging curve under stratum condition

Publications (2)

Publication Number Publication Date
CN102183788A true CN102183788A (en) 2011-09-14
CN102183788B CN102183788B (en) 2012-11-21

Family

ID=44569986

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011101246147A Active CN102183788B (en) 2011-05-13 2011-05-13 Method for synthesizing transverse-wave well-logging curve under stratum condition

Country Status (1)

Country Link
CN (1) CN102183788B (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103278850A (en) * 2013-05-24 2013-09-04 西安石油大学 Transverse wave time difference curve construction method based on coal rock industrial component physical volume model
CN103993871A (en) * 2014-05-08 2014-08-20 中国石油化工股份有限公司 Method and device for processing well logging information of thin interbed stratums in standardization mode
CN104047597A (en) * 2014-01-26 2014-09-17 中国石油大学(华东) Fat gas mud shale stratum well log standardizing method
CN104278991A (en) * 2013-07-09 2015-01-14 中国石油化工股份有限公司 Multivariate well logging computing method for total organic carbon and hydrocarbon generating potential of hydrocarbon source rocks in salt-lake facies
CN107290800A (en) * 2016-03-30 2017-10-24 中国石油化工股份有限公司 Log Forecasting Methodology before practical brill
CN107991705A (en) * 2017-10-13 2018-05-04 中国石油天然气股份有限公司 Log curve correction method and apparatus based on Two-dimensional Statistical feature
CN110361780A (en) * 2019-07-26 2019-10-22 西南交通大学 A kind of seismic wave selection method based on conditional mean input energy spectra
CN110656934A (en) * 2019-10-08 2020-01-07 中国石油天然气股份有限公司 Compacted-stratum-removing comparison method for compact sandstone reservoir
CN110656934B (en) * 2019-10-08 2024-04-30 中国石油天然气股份有限公司 Compaction-removing stratum comparison method for tight sandstone reservoir

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1797032A (en) * 2004-12-29 2006-07-05 中国石油天然气集团公司 Method for determining distribution of lithologic character and liquid by using inversion technique of wave impedance
WO2008033770A1 (en) * 2006-09-12 2008-03-20 Services Petroliers Schlumberger Discriminating natural fracture- and stress-induced sonic anisotropy using a combination of image and sonic logs
CN101354444A (en) * 2007-07-25 2009-01-28 中国石油天然气集团公司 Method for determining formation lithologic character and pore fluid
CN101487898A (en) * 2009-02-27 2009-07-22 中国石油集团川庆钻探工程有限公司 Method for oil gas water recognition by employing longitudinal wave seismic exploration post-stack data

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1797032A (en) * 2004-12-29 2006-07-05 中国石油天然气集团公司 Method for determining distribution of lithologic character and liquid by using inversion technique of wave impedance
WO2008033770A1 (en) * 2006-09-12 2008-03-20 Services Petroliers Schlumberger Discriminating natural fracture- and stress-induced sonic anisotropy using a combination of image and sonic logs
CN101354444A (en) * 2007-07-25 2009-01-28 中国石油天然气集团公司 Method for determining formation lithologic character and pore fluid
CN101487898A (en) * 2009-02-27 2009-07-22 中国石油集团川庆钻探工程有限公司 Method for oil gas water recognition by employing longitudinal wave seismic exploration post-stack data

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《地球科学-中国地质大学学报》 20090731 邵才瑞等 利用常规测井资料基于岩石物理和多矿物分析反演横波速度 第699-706页 1 第34卷, 第4期 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103278850B (en) * 2013-05-24 2014-08-20 西安石油大学 Transverse wave time difference curve construction method based on coal rock industrial component physical volume model
CN103278850A (en) * 2013-05-24 2013-09-04 西安石油大学 Transverse wave time difference curve construction method based on coal rock industrial component physical volume model
CN104278991A (en) * 2013-07-09 2015-01-14 中国石油化工股份有限公司 Multivariate well logging computing method for total organic carbon and hydrocarbon generating potential of hydrocarbon source rocks in salt-lake facies
CN104278991B (en) * 2013-07-09 2016-12-28 中国石油化工股份有限公司 Saline Lake Facies hydrocarbon source rock organic carbon and the polynary well logging computational methods of hydrocarbon potential
CN104047597A (en) * 2014-01-26 2014-09-17 中国石油大学(华东) Fat gas mud shale stratum well log standardizing method
CN104047597B (en) * 2014-01-26 2016-06-08 中国石油大学(华东) Oily mud shale stratum log standardized method
CN103993871A (en) * 2014-05-08 2014-08-20 中国石油化工股份有限公司 Method and device for processing well logging information of thin interbed stratums in standardization mode
CN103993871B (en) * 2014-05-08 2017-02-15 中国石油化工股份有限公司 Method and device for processing well logging information of thin interbed stratums in standardization mode
CN107290800A (en) * 2016-03-30 2017-10-24 中国石油化工股份有限公司 Log Forecasting Methodology before practical brill
CN107991705A (en) * 2017-10-13 2018-05-04 中国石油天然气股份有限公司 Log curve correction method and apparatus based on Two-dimensional Statistical feature
CN107991705B (en) * 2017-10-13 2019-09-06 中国石油天然气股份有限公司 Log curve correction method and apparatus based on Two-dimensional Statistical feature
CN110361780A (en) * 2019-07-26 2019-10-22 西南交通大学 A kind of seismic wave selection method based on conditional mean input energy spectra
CN110656934A (en) * 2019-10-08 2020-01-07 中国石油天然气股份有限公司 Compacted-stratum-removing comparison method for compact sandstone reservoir
CN110656934B (en) * 2019-10-08 2024-04-30 中国石油天然气股份有限公司 Compaction-removing stratum comparison method for tight sandstone reservoir

Also Published As

Publication number Publication date
CN102183788B (en) 2012-11-21

Similar Documents

Publication Publication Date Title
CN102183788B (en) Method for synthesizing transverse-wave well-logging curve under stratum condition
Kaven et al. Surface monitoring of microseismicity at the Decatur, Illinois, CO2 sequestration demonstration site
Barnes et al. Geological sequestration of carbon dioxide in the Cambrian Mount Simon Sandstone: Regional storage capacity, site characterization, and large-scale injection feasibility, Michigan Basin
CN104698492B (en) A kind of method for calculating abnormal formation pressure
CN101787884A (en) Method for judging fluid type of reservoir through acoustic porosity-neutron porosity differential
Block et al. Induced seismicity constraints on subsurface geological structure, Paradox Valley, Colorado
CN107327294A (en) A kind of porosity means of interpretation become based on fine and close oily reservoir under the conditions of matrix parameter
CN109458176A (en) The prediction technique and its application of carbonate reservoir pressure
CN106842326B (en) Sand-mud interbed CO2 geological sequestration time-lapse seismic forward modeling method during non-transverse-wave velocity logging
CN106285642A (en) A kind of formation pore pressure Forecasting Methodology based on seismic data
CN103343687B (en) A kind of well-log information obtains the method for sandstone brine layer equivalence sodium chloride salinity
CN103345002A (en) Method for obtaining sandstone brine layer potassium ion content by means of logging information
CN112443322B (en) Hydrocarbon source rock logging evaluation method based on equivalent saturation
Somma et al. Review of recent drilling projects in unconventional geothermal resources at Campi Flegrei Caldera, Cornubian Batholith, and Williston Sedimentary Basin
Estrada et al. Pilot phase of the Aguada Pichana Este block, gas window
CN104698493B (en) A kind of method for calculating abnormal formation pressure
Gunnarsson 3D modeling in Petrel of geological CO2 storage site
Newhouse et al. Geologic, water-chemistry, and hydrologic data from multiple-well monitoring sites and selected water-supply wells in the Santa Clara Valley, California, 1999–2003
Maldonado et al. Mechanical properties of the Niobrara formation
Mandal et al. Multi-Purpose Utility of Constructing 3D Static Geomechanical Model in the Ichthys Field, Browse Basin
Kudrewicz et al. Subsalt Rotliegend Sediments—A New Challenge for Geothermal Systems in Poland
Mennan Well Log Interpretation and 3D Reservoir Property Modeling of the Maui-B Field, Taranaki Basin, New Zealand
Rolfs Integrated geomechanical, geophysical, and geochemical analysis of the Bakken Formation, Elm Coulee field, Williston Basin, Montana
Papiernik et al. Hydrocarbon resources assessment for selected area of Lower Paleozoic basin in Poland based on spatial modeling
Mo et al. Prediction of shale prospectivity from seismically-derived reservoir and completion qualities: Application to a shale-gas field, Horn River Basin, Canada

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant