CN110133722A - A kind of application velocity of sound qualitative recognition mud distribution type new method - Google Patents

A kind of application velocity of sound qualitative recognition mud distribution type new method Download PDF

Info

Publication number
CN110133722A
CN110133722A CN201910493729.XA CN201910493729A CN110133722A CN 110133722 A CN110133722 A CN 110133722A CN 201910493729 A CN201910493729 A CN 201910493729A CN 110133722 A CN110133722 A CN 110133722A
Authority
CN
China
Prior art keywords
velocity
sandstone
shale
longitudinal wave
sound
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
CN201910493729.XA
Other languages
Chinese (zh)
Other versions
CN110133722B (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.)
Karamay Onco Energy Technology Co Ltd
Original Assignee
Karamay Onco Energy Technology 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 Karamay Onco Energy Technology Co Ltd filed Critical Karamay Onco Energy Technology Co Ltd
Priority to CN201910493729.XA priority Critical patent/CN110133722B/en
Publication of CN110133722A publication Critical patent/CN110133722A/en
Application granted granted Critical
Publication of CN110133722B publication Critical patent/CN110133722B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/303Analysis for determining velocity profiles or travel times
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/622Velocity, density or impedance
    • G01V2210/6222Velocity; travel time
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/624Reservoir parameters

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The present invention provides a kind of application velocity of sound qualitative recognition mud distribution type new methods, it is related to oil and gas prospect technical field, this method is distributed the upper limit as sandstone velocity of longitudinal wave including the use of the velocity of longitudinal wave that CCT model calculates contact argillaceous sandstone, draws the shale content and velocity of longitudinal wave cross plot of contact shale, structure and laminar shale;Dispersed shale sandstone velocity of longitudinal wave is calculated as sandstone velocity of longitudinal wave using HMHS model and is distributed lower limit, draws dispersed shale content and velocity of longitudinal wave cross plot;The distributed areas for marking dispersed shale sandstone, contact argillaceous sandstone, structure sandstone and the laminar shale sandstone velocity of sound obtain argillaceous sandstone velocity of sound distribution plate;The velocity of sound of practical argillaceous sandstone and shale content data are substituted into plate, determine mud distribution form.This method can in qualitative recognition argillaceous sandstone shale distribution pattern, it is representative it is strong, the period is short, can large-scale application in well logging, seismic evaluation etc..

Description

A kind of application velocity of sound qualitative recognition mud distribution type new method
Technical field
The present invention relates to oil and gas prospect technical fields, and in particular to a kind of application velocity of sound qualitative recognition shale point Cloth type new method.
Background technique
Need to consider the influence of mud distribution type and content to the sandstone velocity of sound in the evaluation of actual geophysics.It is logical Often, can using drilling extracting core laboratory under the microscope, grain size analysis determine mud distribution type, but there are samples for this method The technological deficiencies such as this representativeness is limited, experimental period is long, can not large-scale application in well logging, seismic evaluation.
Chinese patent CN109324171A discloses a kind of sedimentary facies qualitative recognition method based on lithology statistics, specific to wrap Include: data collection, data processing, parameter calculate, establish sedimentary facies qualitative recognition plate and the parameter according to each computing unit Calculated result, wherein data collection includes collecting work area prospect pit lithologic log data and known deposition phase data;Data processing Continuous continual stratomere corresponding to a kind of type of sedimentary facies divided using on individual well is rejected as a computing unit Thickness is less than the computing unit of 20m, rejects the invalid rock stratum in each computing unit, remaining rock stratum is divided into sandstone, mud stone two Major class;Parameter calculates than two parameters of sandstone frequency and sandy ground for calculating separately each computing unit;According to parameter calculated result It then crosses in sandstone frequency-sandy ground ratio and carries out cultellation on figure, the type of sedimentary facies and cultellation position further according to computing unit are drawn Make the identification region of different type of sedimentary facies.The working efficiency of sedimentary facies recognition work and accurate is greatly improved in the invention Rate;Have the advantages that speed is fast, at low cost, high-efficient, Subjective Intervention influence it is small.But this method is mainly for continuity stratum Section, and need to reject invalid rock stratum, while there is still a need for subjective judgement is carried out to rock stratum type in advance, application range has office It is sex-limited.
Chinese patent CN106842359B discloses a kind of side using wave impedance qualitative recognition complexity sand-conglomerate body lithology Method, comprising: determine different lithology in array induction, density, interval transit time, middle sub-aperture in conjunction with log data using rock core information Distribution on these borehole log datas of porosity;Carry out the rough segmentation of lithology;Lithology subdivision is carried out, and utilizes wave impedance, into One step determines calcareous sandstone, grey matter conglomerate;By above step, seven class lithology are finely determined.This method has preferable Operability is conducive to promote, and the qualitative recognition for complicated sand-conglomerate body lithology provides practicable new method.But this method needs Want the data such as detailed rock core information, including array induction, density, interval transit time and neutron porosity, need for rock stratum into The parameter of row various aspects measures, complex.
So far, mud distribution kind identification method there are it is complex, multiple parameters need to be relied on and subjectivity is sentenced The problems such as disconnected, it is therefore desirable to find it is a kind of it is representative it is strong, the period is short, subjectivity rely on it is small, can large-scale application in well logging, Shake the recognition methods of evaluation, the distribution pattern to shale in qualitative recognition argillaceous sandstone.
Summary of the invention
In view of the problems of the existing technology the present invention, it is new to provide a kind of application velocity of sound qualitative recognition mud distribution type Method.This method is representative strong, the period is short, can large-scale application in well logging, the recognition methods of seismic evaluation, to qualitative knowledge The distribution pattern of shale in other argillaceous sandstone.
To achieve the above object, The technical solution adopted by the invention is as follows:
A kind of application velocity of sound qualitative recognition mud distribution type new method, comprising the following steps:
(1), it is distributed the upper limit as rock velocity of longitudinal wave using the velocity of longitudinal wave that model calculates contact argillaceous sandstone, draws contact Shale, the shale content of structure and laminar shale and velocity of longitudinal wave cross plot;
(2), dispersed shale sandstone velocity of longitudinal wave is calculated as sandstone velocity of longitudinal wave using model and be distributed lower limit, draw dispersion Shale content and velocity of longitudinal wave cross plot;
(3), dispersed shale sandstone, contact argillaceous sandstone, structure sandstone and the laminar shale sandstone velocity of sound are marked respectively Distributed areas, obtain the argillaceous sandstone velocity of sound distribution plate;
(4), the velocity of sound of practical argillaceous sandstone and shale content data are substituted into plate, and according to the practical argillaceous sandstone velocity of sound Data and boundary relative position determine mud distribution form.
Further, model described in step (1) is CCT model, wherein the application method of CCT model is referring to paper Xue-Hui H,Jun-Xin G,Feng-Bi L,et al.Generalization of the Expression of Cementation Radius in Contact Cement Theory and Its Application[J].Chinese Journal of Geophysics,2015,57(4):439-449.
Further, model described in step (2) is HMHS model, wherein the application method of HMHS model is referring to paper Dvorkin J P,Nur A M.Elasticity of high-porous sandstones:Theory for two North Sea data sets[J].Geophysics,1996,61(5):1363-1370.
Further, the velocity of longitudinal wave of dispersed shale described in step (2) is between Vlow-aVlowBetween;Institute in step (1) The velocity of longitudinal wave of contact shale is stated between cVup-VupBetween;The velocity of longitudinal wave of structure and laminar shale described in step (1) Between dVup-eVupBetween.Wherein, the VupThe upper limit, V are distributed for sandstone velocity of longitudinal wavelowFor under the distribution of sandstone velocity of longitudinal wave Limit.The shale content of the contact shale is less than b;The shale content of the structure and laminar shale is greater than b, wherein b is The interpolation of one shale content is distributed upper limit V in the velocity of soundupOn, it is contact shale that b, which has previously been thought that, and b is considered structure or layer later Shape shale.
Further, described a, b, c, d and e are undetermined coefficient, by empirically determined, in some specific embodiments In, the range of a is 1-1.6, and the range of the b is 0-0.15, and the range of the c is 0.8-1, and the range of the d is The range of 0.8-1, the e are 1-1.3.Preferably, a is 1.5, and the b is 0.1, and the c is 0.9, and the d is 0.9, The e is 1.1.
It is worth noting that the present invention is mixing shale, real velocity of sound distribution by white space surrounded in plate Illustrate mud distribution based on dispersed shale closer to lower limit;Real velocity of sound illustrates mud distribution with structure closer to the upper limit Or based on laminar shale;Real velocity of sound illustrates mud distribution closer to contact shale region to contact based on shale.
It is obtained by the present invention to have the technical effect that
1. method of the invention can in qualitative recognition argillaceous sandstone shale distribution pattern, be suitable for mud distribution type More single situation.
2. method of the invention it is representative it is strong, the period is short, subjectivity rely on it is small, can large-scale application in well logging, earthquake Evaluation provides reference for earth-layer fine explanation, reduces and explain risk.
Detailed description of the invention
Fig. 1 is that the argillaceous sandstone velocity of sound is distributed plate;
Fig. 2 is the distribution of the argillaceous sandstone velocity of sound and argillaceous sandstone acoustic speed plate, wherein scatterplot is Han Dehua (Han, De- Hua the argillaceous sandstone acoustic speed data) published.
Specific embodiment
A kind of application velocity of sound qualitative recognition mud distribution type new method, comprising the following steps:
(1), the upper limit is distributed as sandstone velocity of longitudinal wave according to the velocity of longitudinal wave that CCT model calculates contact argillaceous sandstone, drawn The shale content and velocity of longitudinal wave cross plot of system contact shale, structure and laminar shale;
(2), the dispersed shale sandstone velocity of sound is calculated as sandstone velocity of longitudinal wave according to HMHS model and is distributed lower limit, draw dispersion Shale content and velocity of longitudinal wave cross plot;
(3), point of dispersed shale sandstone, contact argillaceous sandstone, structure sandstone and the laminar shale sandstone velocity of sound is marked Cloth region obtains argillaceous sandstone velocity of sound distribution plate;Wherein, the velocity of longitudinal wave of dispersed shale is between Vlow-aVlowBetween, contact The velocity of longitudinal wave of shale is between cVup-VupBetween, shale content is less than b;The velocity of longitudinal wave of structure or laminar shale between dVup-eVupBetween, shale content is greater than b.Wherein, a, b, c, d and e are undetermined coefficient, and a 1.5, b 0.1, c 0.9, d are 0.9, e 1.1;VlowRepresent sandstone velocity of longitudinal wave distribution lower limit, VupRepresent the sandstone velocity of longitudinal wave distribution upper limit.According to as above ginseng The argillaceous sandstone velocity of sound distribution plate that number is established is as shown in Figure 1;Wherein, the setting of parameters are as follows: quartz body product module amount K= 37GPa;Quartzy shear modulus G=44GPa;Quartzy density p=2.65g/cm3;Shale bulk modulus Kc=25GPa;Shale shearing Modulus Gc=9GPa;Shale density pc=2.55g/cm3;Critical porosityLigancy n=9;Shale contentWherein,For the practical porosity of sandstone.
(4), by paper Han D, Nur A, Morgan D.Effects of porosity and clay content On wave velocities in sandstones [J] .Geophysics, 1986,51 (11): practical shale sand in 2093. The velocity of sound and shale content data of rock substitute into plate (i.e. Fig. 1), obtain Fig. 2, and according to practical argillaceous sandstone sound velocity data and Boundary relative position determines mud distribution form.
As shown in Figure 2, in the present invention, the region that the velocity of sound and shale content data of practical argillaceous sandstone are fallen into largely is Structure or laminar shale, for mud distribution based on structure or stratiform, mud distribution type is more single in argillaceous sandstone sample One, the empirical equation that the relationship and Han Dehua for showing approximately linear sum up the argillaceous sandstone velocity of sound of approximately linear is mutual Verifying, the above results show that the conclusion that method of the invention obtains and the Han Dehua empirical equation put forward and experimental data are mutual It verifies and supplements, the method in the present invention is more accurate.
Finally it should be noted that the above content is merely illustrative of the technical solution of the present invention, rather than the present invention is protected The limitation of range, the simple modification or equivalent replacement that those skilled in the art carry out technical solution of the present invention, All without departing from the spirit and scope of technical solution of the present invention.

Claims (10)

1. a kind of application velocity of sound qualitative recognition mud distribution type new method, it is characterised in that: the following steps are included:
(1), it is distributed the upper limit as sandstone velocity of longitudinal wave using the velocity of longitudinal wave that model calculates contact argillaceous sandstone, draws contact mud Matter, the shale content of structure and laminar shale and velocity of longitudinal wave cross plot;
(2), dispersed shale sandstone velocity of longitudinal wave is calculated as sandstone velocity of longitudinal wave using model and be distributed lower limit, draw dispersed shale Content and velocity of longitudinal wave cross plot;
(3), point of dispersed shale sandstone, contact argillaceous sandstone, structure sandstone and the laminar shale sandstone velocity of sound is marked respectively Cloth region obtains argillaceous sandstone velocity of sound distribution plate;
(4), the velocity of sound of practical argillaceous sandstone and shale content data are substituted into plate, and according to practical argillaceous sandstone sound velocity data Mud distribution form is determined with boundary relative position.
2. according to the method described in claim 1, it is characterized by: model described in step (1) is CCT model.
3. according to the method described in claim 1, it is characterized by: model described in step (2) is HMHS model.
4. according to the method described in claim 1, it is characterized by: the velocity of longitudinal wave of dispersed shale described in step (2) between Vlow-aVlowBetween, the VlowLower limit is distributed for sandstone velocity of longitudinal wave.
5. according to the method described in claim 1, it is characterized by: described in step (1) contact shale velocity of longitudinal wave between cVup-VupBetween, the VupThe upper limit is distributed for sandstone velocity of longitudinal wave.
6. according to the method described in claim 1, it is characterized by: the velocity of longitudinal wave of structure described in step (1) between dVup-eVupBetween, the VupThe upper limit is distributed for sandstone velocity of longitudinal wave.
7. according to the method described in claim 1, it is characterized by: the velocity of longitudinal wave of laminar shale described in step (1) between dVup-eVupBetween, the VupThe upper limit is distributed for sandstone velocity of longitudinal wave.
8. according to the method described in claim 1, it is characterized by: the shale content for contacting shale described in step (1) is less than b。
9. according to the method described in claim 1, it is characterized by: the mud of structure and laminar shale described in step (1) Matter content is greater than b.
10. according to method described in claim 4-9 any one, it is characterised in that: described a, b, c, d and e are undetermined coefficient, The range of a is 1-1.6, and the range of the b is 0-0.15, and the range of the c is 0.8-1, and the range of the d is 0.8-1, The range of the e is 1-1.3.
CN201910493729.XA 2019-06-06 2019-06-06 Method for qualitatively identifying mud distribution type by using sound velocity Active CN110133722B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910493729.XA CN110133722B (en) 2019-06-06 2019-06-06 Method for qualitatively identifying mud distribution type by using sound velocity

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910493729.XA CN110133722B (en) 2019-06-06 2019-06-06 Method for qualitatively identifying mud distribution type by using sound velocity

Publications (2)

Publication Number Publication Date
CN110133722A true CN110133722A (en) 2019-08-16
CN110133722B CN110133722B (en) 2019-12-24

Family

ID=67580471

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910493729.XA Active CN110133722B (en) 2019-06-06 2019-06-06 Method for qualitatively identifying mud distribution type by using sound velocity

Country Status (1)

Country Link
CN (1) CN110133722B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112987096A (en) * 2021-03-15 2021-06-18 中国石油大学(华东) Method for calculating sound velocity of high-argillaceous sandstone

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102426390A (en) * 2011-10-21 2012-04-25 中国石油大学(北京) Method for determining reserve volume of nonhomogeneous sandstone reservoir
CN104570127A (en) * 2013-10-09 2015-04-29 中国石油化工股份有限公司 Method for simultaneously inverting porosity and shaliness by using seismic wave velocity
US20150120197A1 (en) * 2013-10-30 2015-04-30 Chevron U.S.A. Inc. Method of Transforming Reservoir Properties to a Seismic Attribute for Hydrocarbon and Lithology Identification
CN104898161A (en) * 2014-03-05 2015-09-09 中国石油集团东方地球物理勘探有限责任公司 Effective sandstone predicting method based on logging response simulator
CN104950331A (en) * 2015-06-30 2015-09-30 中国石油天然气集团公司 Earthquake prediction method for porosity and shale content of sand shale reservoir
CN108399270A (en) * 2017-02-08 2018-08-14 中国石油化工股份有限公司 A kind of method of anisotropy shale ratio in determining shale formation

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102426390A (en) * 2011-10-21 2012-04-25 中国石油大学(北京) Method for determining reserve volume of nonhomogeneous sandstone reservoir
CN104570127A (en) * 2013-10-09 2015-04-29 中国石油化工股份有限公司 Method for simultaneously inverting porosity and shaliness by using seismic wave velocity
US20150120197A1 (en) * 2013-10-30 2015-04-30 Chevron U.S.A. Inc. Method of Transforming Reservoir Properties to a Seismic Attribute for Hydrocarbon and Lithology Identification
CN104898161A (en) * 2014-03-05 2015-09-09 中国石油集团东方地球物理勘探有限责任公司 Effective sandstone predicting method based on logging response simulator
CN104950331A (en) * 2015-06-30 2015-09-30 中国石油天然气集团公司 Earthquake prediction method for porosity and shale content of sand shale reservoir
CN108399270A (en) * 2017-02-08 2018-08-14 中国石油化工股份有限公司 A kind of method of anisotropy shale ratio in determining shale formation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王飞: "交叉偶极声波测井资料处理及在致密气层评价中的应用研究", 《中国博士学位论文全文数据库 基础科学辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112987096A (en) * 2021-03-15 2021-06-18 中国石油大学(华东) Method for calculating sound velocity of high-argillaceous sandstone
CN112987096B (en) * 2021-03-15 2021-11-26 中国石油大学(华东) Method for calculating sound velocity of high-argillaceous sandstone

Also Published As

Publication number Publication date
CN110133722B (en) 2019-12-24

Similar Documents

Publication Publication Date Title
CN107703560B (en) shale lithofacies fine identification method based on triple information
CN105467464B (en) The effective connected component of Analyze of Nearshore Subaqueous Fans glutenite is divided and control methods
Friend et al. Sedimentation of the Wood Bay Formation (Devonian) of Spitsbergen: regional analysis of a late orogenic basin
CN109613612A (en) A kind of carbonate rock particle beach meticulous depiction and prediction technique
CN105467463B (en) Analyze of Nearshore Subaqueous Fans fluvial deposition cell cube maximum extended distance Quantitative prediction methods
CN106556866B (en) A method of using rock physics template to reservoir quantitative interpretation
CN105701319B (en) A kind of microfacies modeling method under horizontal well rule well pattern
CN103675907A (en) AVO inversion hydrocarbon detection method based on petrographic constraints
CN103529474A (en) Method for realizing accurate description lithofacies by using lithology subdivision
CN104047598A (en) Method for predicating productivity of nonhomogeneity ancient karst carbonate reservoir
CN106707337A (en) Method for sedimentary system interpretation based on fractional-frequency seismic reflection energy stratigraphic slice
CN105467466B (en) A kind of compact reservoir Diagenetic Facies Forecasting Methodology based on multi-scale information
CN112034514B (en) Shale oil-gas system based on seismic sedimentology theory and identification method thereof
CN109298464A (en) Tight sandstone reservoir Diagenetic Facies Logging Identification Method and device
CN106033127B (en) Crustal stress azimuthal seismic Forecasting Methodology based on shear wave velocity rate of change
CN110019594A (en) A kind of geologic datebase method for building up based on numerical simulation
CN110133722A (en) A kind of application velocity of sound qualitative recognition mud distribution type new method
CN102914797A (en) Method and device for acquiring anisotropy coefficient of stratum
CN105468890B (en) Analyze of Nearshore Subaqueous Fans fluvial deposition cell cube width granularity Forecasting Methodology
CN111198406A (en) Lithology recognition method for factor analysis logging of red reservoir
CN115857047B (en) Comprehensive prediction method for earthquake reservoir
Deng et al. Factors controlling reservoir quality of a retreating delta-front in shallow-water lacustrine in the Songliao Basin, Northeast China
Fan et al. Semi-supervised learning–based petrophysical facies division and “Sweet Spot” identification of low-permeability sandstone reservoir
CN113820754B (en) Deep tight sandstone reservoir evaluation method based on artificial intelligence identification of reservoir lithofacies
CN109490988A (en) Establish the method for being suitable for the rock physics new model of hard rock

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant