CN109374743A - 一种致密砂岩杨氏模量修正方法及装置 - Google Patents
一种致密砂岩杨氏模量修正方法及装置 Download PDFInfo
- Publication number
- CN109374743A CN109374743A CN201811186028.3A CN201811186028A CN109374743A CN 109374743 A CN109374743 A CN 109374743A CN 201811186028 A CN201811186028 A CN 201811186028A CN 109374743 A CN109374743 A CN 109374743A
- Authority
- CN
- China
- Prior art keywords
- modulus
- young
- content
- tight sand
- kinetic youngs
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/011—Velocity or travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/024—Mixtures
- G01N2291/0245—Gases in porous solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/02827—Elastic parameters, strength or force
Landscapes
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Acoustics & Sound (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种致密砂岩杨氏模量修正方法及装置,该方法首先根据阵列声波波速,计算动态杨氏模量;然后采集致密砂岩的硅元素含量、铁元素含量、泥质含量和孔隙度;根据硅元素含量、铁元素含量、泥质含量和孔隙度,对计算出来的动态杨氏模量进行修正,得到修正后的杨氏模量。本发明对计算出来的动态杨氏模量加入反映石英含量的硅元素含量、反映岩屑含量的铁元素含量、以及反映塑性的泥质含量和孔隙度作为修正系数,对计算出来的动态杨氏模量进行修正,使得修正后的动态杨氏模量更准确,与静态杨氏模量的相关性较高。
Description
技术领域
本发明属于油气技术领域,具体涉及一种致密砂岩杨氏模量修正方法及装置。
背景技术
致密砂岩气是发育于低孔、特低渗透砂岩中的一种非常规油气资源,油气储量丰富,但由于其储层致密、孔隙结构复杂、非均质强,开采前需要进行压裂施工以获取更多的天然气。压裂施工中岩石的杨氏模量对压裂施工参数的选择具有重要的意义。
常规动态杨氏模量的计算方法为,通过阵列声波、交叉偶极子声波等特殊测井获取纵、横波速度,计算动态杨氏模量,但致密砂岩受陆相沉积多期河道叠置影响,矿物成分变化大,孔隙结构复杂,薄隔层受围岩等因素影响,单一利用声波计算动态杨氏模量误差大。
例如,鄂尔多期盆地北部气田储层岩石随着岩屑含量增加,声波时差逐渐增大,动态杨氏模量变小,而实验室分析岩屑砂岩动态杨氏模量明显变大,通过阵列声波计算得到的动态杨氏模量与实验室分析的杨氏模量差异明显。岩石的杨氏模量随泥质含量先增加再减小,单一利用声波计算得到的动态杨氏模量不能反映泥质含量对杨氏模量的影响。
发明内容
本发明的目的在于提供一种致密砂岩杨氏模量修正方法及装置,用以解决单一利用声波计算动态杨氏模量造成的误差大的问题。
为解决上述技术问题,本发明的技术方案为:
本发明提供了一种致密砂岩杨氏模量修正方法,包括如下步骤:
1)根据阵列声波波速,计算动态杨氏模量;
2)采集致密砂岩的硅元素含量、铁元素含量、泥质含量和孔隙度;根据硅元素含量、铁元素含量、泥质含量和孔隙度,对计算出来的动态杨氏模量进行修正,得到修正后的杨氏模量。
本发明还提供了一种致密砂岩杨氏模量修正装置,包括处理器,所述处理器用于执行指令实现如下方法:
1)根据阵列声波波速,计算动态杨氏模量;
2)采集致密砂岩的硅元素含量、铁元素含量、泥质含量和孔隙度;根据硅元素含量、铁元素含量、泥质含量和孔隙度,对计算出来的动态杨氏模量进行修正,得到修正后的杨氏模量。
本发明的有益效果:
该方法及装置引入反映石英含量的硅元素和反映岩屑含量的铁元素作为修正系数,对计算出来的动态杨氏模量进行修正,同时,在对致密砂岩进行实验过程中,通过分析样本点,发现杨氏模量随泥质含量的增加先增加后减小,孔隙度与杨氏模量成负相关,故根据此特性对计算出来的动态杨氏模量再加入反映塑性的泥质含量和孔隙度作为修正系数,对计算出来的动态杨氏模量进行修正,使得修正后的动态杨氏模量更准确,与静态杨氏模量的相关性较高。
作为方法及装置的进一步改进,为了进一步提高动态杨氏模量的准确性,所述修正后的杨氏模量为:
其中,E2为修正后的杨氏模量,E1为计算出来的动态杨氏模量,k1、k2、k3、k4和k5均为系数,φ为孔隙度,Vsh为泥质含量,Si为硅元素含量,Fe为铁元素含量。
作为方法及装置的进一步改进,为了进一步提高动态杨氏模量的准确性,k1=-0.24,k2=-0.005,k3=0.245,k4=-0.003,k5=0.029。
作为方法及装置的进一步改进,为了得到准确的动态杨氏模量,所述动态杨氏模量为:
其中,E1为动态杨氏模量,ρ为致密砂岩的密度,Vs为横波波速,Vp为纵波波速。
附图说明
图1是E0与采用现有技术方法得到的E1的相关性示意图;
图2是E0与采用本发明方法进行修正得到的E2的相关性示意图。
具体实施方式
方法实施例
首先,根据阵列声波波速,计算动态杨氏模量E1:
其中,E1为动态杨氏模量,ρ为致密砂岩的密度,Vs为横波波速,Vp为纵波波速。
具体如何计算动态杨氏模量的过程为现有技术,这里不再赘述。
然后,采集致密砂岩的泥质含量、孔隙度、硅元素含量(Si)、铁元素含量(Fe),引入测井孔隙度、泥质含量、元素录井Fe、元素录井Si四个校正因素,分析各因素对杨氏模量的影响,根据式(2),对计算出来的动态杨氏模量进行修正,得到修正后的杨氏模量E2:
其中,E2为修正后的杨氏模量,E1为计算出来的动态杨氏模量,k1、k2、k3、k4和k5均为系数,φ为孔隙度,Vsh为泥质含量,Si为硅元素含量,Fe为铁元素含量。
其中,系数k1、k2、k3、k4和k5可选择如下:k1=-0.24,k2=-0.005,k3=0.245,k4=-0.003,k5=0.029,以使修正后的杨氏模量更准确。即:
其中,E2为修正后的杨氏模量,E1为计算出来的动态杨氏模量,φ为孔隙度,Vsh为泥质含量,Si为硅元素含量,Fe为铁元素含量。
当然,还可根据实验结果将系数k1、k2、k3、k4和k5进行调整与修改,以使修正后的杨氏模量E2与静态杨氏模量E0相关性更高。其中,具体获得静态杨氏模量E0的方法,可采用现有技术中的方法,对致密砂岩做岩石力学实验得到。
下面通过具体的实验对上述方法的有效性进行验证。
选取已测阵列声波或者交叉偶极子声波的井,选取系列岩心做岩石力学实验,获得岩石的静态杨氏模量E0,具体内容见表1。
通过岩心归位,确定岩石力学实验样品点的阵列声波波速值,利用阵列声波波速,根据式(1)计算动态杨氏模量E1,与对应深度岩石力学实验获取的静态杨氏模量E0做相关性分析,如图1所示,二者总体上呈线性正相关,拟合后的直线方程为y=0.6467x-0.3199,其中,椭圆形圈住部分为岩屑砂岩点、泥质粉砂岩点。但是由于孔隙、岩屑和泥质点影响,从图1中可以看出,部分样本点偏离较大,相关系数R2=0.604,E0和E1的相关性较差。此时,需要引入相关参数对E1进行校正。
通过分析样本点,发现杨氏模量随泥质含量的增加先增加后减小、孔隙度与杨氏模量呈负相关的特性,加入泥质含量和孔隙度对动态杨氏模量E1进行校正,发现修正后的E1与E0的相关性变好。
但岩屑砂岩和泥质砂岩受矿物成分影响较大,如本样本中的泥质砂岩和岩屑砂岩偏离较大,故引入反映岩石矿物成分的Fe、Si元素,并根据式(2)对计算出来的动态杨氏模量E1进行校正,校正后得到修正后的杨氏模量E2。此时,如图2所示,再对E0和E2做相关性分析,得到的拟合后的直线方程为y=0.6447x+1.4991,相关系数R2=0.8679,相关性明显变好。
表1样本点数据
装置实施例
本发明还提供了一种致密砂岩杨氏模量修正装置,该装置实际上为计算机等具备数据处理能力的设备,该设备包括处理器,该处理器可以是通用处理器,还可以是数字信号处理器、专用集成电路等,该处理器用于执行指令实现本发明的致密砂岩杨氏模量修正方法,具体方法见上述介绍的方法实施例,这里不再赘述。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (8)
1.一种致密砂岩杨氏模量修正方法,其特征在于,包括如下步骤:
1)根据阵列声波波速,计算动态杨氏模量;
2)采集致密砂岩的硅元素含量、铁元素含量、泥质含量和孔隙度;根据硅元素含量、铁元素含量、泥质含量和孔隙度,对计算出来的动态杨氏模量进行修正,得到修正后的杨氏模量。
2.根据权利要求1所述的致密砂岩杨氏模量修正方法,其特征在于,所述修正后的杨氏模量为:
其中,E2为修正后的杨氏模量,E1为计算出来的动态杨氏模量,k1、k2、k3、k4和k5均为系数,φ为孔隙度,Vsh为泥质含量,Si为硅元素含量,Fe为铁元素含量。
3.根据权利要求2所述的致密砂岩杨氏模量修正方法,其特征在于,k1=-0.24,k2=-0.005,k3=0.245,k4=-0.003,k5=0.029。
4.根据权利要求1所述的致密砂岩杨氏模量修正方法,其特征在于,所述动态杨氏模量为:
其中,E1为动态杨氏模量,ρ为致密砂岩的密度,Vs为横波波速,Vp为纵波波速。
5.一种致密砂岩杨氏模量修正装置,其特征在于,包括处理器,所述处理器用于执行指令实现如下方法:
1)根据阵列声波波速,计算动态杨氏模量;
2)采集致密砂岩硅元素含量、铁元素含量、泥质含量和孔隙度;根据硅元素含量、铁元素含量、泥质含量和孔隙度,对计算出来的动态杨氏模量进行修正,得到修正后的杨氏模量。
6.根据权利要求5所述的致密砂岩杨氏模量修正装置,其特征在于,所述修正后的杨氏模量为:
其中,E2为修正后的杨氏模量,E1为计算出来的动态杨氏模量,k1、k2、k3、k4和k5均为系数,φ为孔隙度,Vsh为泥质含量,Si为硅元素含量,Fe为铁元素含量。
7.根据权利要求6所述的致密砂岩杨氏模量修正装置,其特征在于,k1=-0.24,k2=-0.005,k3=0.245,k4=-0.003,k5=0.029。
8.根据权利要求5所述的致密砂岩杨氏模量修正装置,其特征在于,所述动态杨氏模量为:
其中,E1为动态杨氏模量,ρ为致密砂岩的密度,Vs为横波波速,Vp为纵波波速。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811186028.3A CN109374743B (zh) | 2018-10-11 | 2018-10-11 | 一种致密砂岩杨氏模量修正方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811186028.3A CN109374743B (zh) | 2018-10-11 | 2018-10-11 | 一种致密砂岩杨氏模量修正方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109374743A true CN109374743A (zh) | 2019-02-22 |
CN109374743B CN109374743B (zh) | 2021-02-12 |
Family
ID=65403097
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811186028.3A Active CN109374743B (zh) | 2018-10-11 | 2018-10-11 | 一种致密砂岩杨氏模量修正方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109374743B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278614A (zh) * | 2013-04-28 | 2013-09-04 | 北京大学 | 一种动静态岩石力学参数校正方法及装置 |
CN103293563A (zh) * | 2013-05-13 | 2013-09-11 | 中国石油大学(华东) | 一种确定油气储层岩石裂隙发育度和流体性质的方法 |
CN108460204A (zh) * | 2018-02-26 | 2018-08-28 | 浙江中科电声研发中心 | 一种通过扬声器振动部件的受力和位移反推其材料动态力学参数的方法 |
-
2018
- 2018-10-11 CN CN201811186028.3A patent/CN109374743B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278614A (zh) * | 2013-04-28 | 2013-09-04 | 北京大学 | 一种动静态岩石力学参数校正方法及装置 |
CN103293563A (zh) * | 2013-05-13 | 2013-09-11 | 中国石油大学(华东) | 一种确定油气储层岩石裂隙发育度和流体性质的方法 |
CN108460204A (zh) * | 2018-02-26 | 2018-08-28 | 浙江中科电声研发中心 | 一种通过扬声器振动部件的受力和位移反推其材料动态力学参数的方法 |
Non-Patent Citations (4)
Title |
---|
E. A. EISSA 等: "Technical Note Relation Between Static and Dynamic Young’s Moduli of Rocks", 《ROCK MECH。》 * |
P. S. THEOCARIS A N D C. LAZOPOULOS: "AN I M P R O V E D M E T H O D FOR THE EVALUATION OF M O D E I DYNAMIC STRESS INTENSITY FACTORS", 《IMPACT ENONA》 * |
SHIYU XU 等: "A physical model for shear-wave velocity prediction", 《GEOPHYSICAL PROSPECTING》 * |
耿乃光 等: "岩石动态与静态杨氏模量的对比测量研究", 《第二节全国岩石动力学学术会议》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109374743B (zh) | 2021-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Grew et al. | Evolution of structural complexity in boron minerals | |
Pessagno | Upper Cretaceous stratigraphy of the western Gulf Coast area of Mexico, Texas, and Arkansas | |
CN105221141B (zh) | 一种泥页岩脆性指数预测方法 | |
Fildani et al. | U-Pb zircon ages from the southwestern Karoo Basin, South Africa—Implications for the Permian-Triassic boundary | |
Lindskog et al. | Lower–Middle Ordovician carbon and oxygen isotope chemostratigraphy at Hällekis, Sweden: implications for regional to global correlation and palaeoenvironmental development | |
Reches | Development of monoclines: Part I. Structure of the Palisades Creek branch of the East Kaibab monocline, Grand Canyon, Arizona | |
CN103984027A (zh) | 基于椭球体双重孔隙模型的岩石纵波速度预测方法 | |
Shi et al. | Weakening laws of rock uniaxial compressive strength with consideration of water content and rock porosity | |
FU et al. | Transitional source rock and its contribution to hydrocarbon accumulation in superimpose rift-subsidence basin of northern South China Sea: Taking Baiyun Sag of Zhu Ⅱ Depression as an example | |
US20230358918A1 (en) | Forward physical simulation method for seismic response characteristics of marine natural gas hydrate system | |
CN112814669A (zh) | 一种页岩油藏全生命周期采收率预测方法和系统 | |
Bissell | Tectonic control of late Paleozoic and early Mesozoic sedimentation near the hinge line of the Cordilleran miogeosynclinal belt | |
Hupp et al. | Effects of size-dependent sediment mixing on deep-sea records of the Paleocene-Eocene Thermal Maximum | |
Stone et al. | Quantifying compaction, pressure solution and quartz cementation in moderately-and deeply-buried quartzose sandstones from the Greater Green River Basin, Wyoming | |
Kahraman et al. | Estimating the wet-rock P-wave velocity from the dry-rock P-wave velocity for pyroclastic rocks | |
Song et al. | Correlations between rock mass intactness index (K v) and other rock mass classification indices (RMR 89 system and GSI) | |
CN109374743A (zh) | 一种致密砂岩杨氏模量修正方法及装置 | |
Zhang et al. | Chloride transport in undersea concrete tunnel | |
Wu et al. | New biostratigraphic evidence of Late Permian to Late Triassic deposits from Central Tibet and their paleogeographic implications | |
Jian-Zhong et al. | One-dimensional viscous behavior of clay and its constitutive modeling | |
Zhang et al. | Prediction of diagenetic facies via well logs and petrophysical properties in tight sandstone from Zhu-III sag: Pearl River mouth basin, South China Sea | |
CN110161208B (zh) | 页岩非均质性定量表征方法 | |
Li et al. | New zircon U-Pb age and its restriction on the warming time of the interglacial paleoclimate during the Cryogenian in the Yangtze Block | |
Yu et al. | Influences of Tibetan Plateau uplift on provenance evolution of the paleo-Pearl River | |
G. Ghedan et al. | An improved approach for generating saturation functions for simulation models using dynamic rock types |
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 |