CN105588883A - 三维岩石力学参数获取方法和*** - Google Patents
三维岩石力学参数获取方法和*** Download PDFInfo
- Publication number
- CN105588883A CN105588883A CN201410641421.2A CN201410641421A CN105588883A CN 105588883 A CN105588883 A CN 105588883A CN 201410641421 A CN201410641421 A CN 201410641421A CN 105588883 A CN105588883 A CN 105588883A
- Authority
- CN
- China
- Prior art keywords
- wave
- velocity
- time difference
- rock
- density
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 148
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000003068 static effect Effects 0.000 claims abstract description 34
- 230000001915 proofreading effect Effects 0.000 claims description 15
- 238000012937 correction Methods 0.000 claims description 9
- 238000007619 statistical method Methods 0.000 claims description 9
- 239000004575 stone Substances 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 7
- 241000208340 Araliaceae Species 0.000 claims description 5
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 5
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 235000008434 ginseng Nutrition 0.000 claims description 5
- 239000004576 sand Substances 0.000 claims description 5
- 238000000611 regression analysis Methods 0.000 claims description 4
- 239000004568 cement Substances 0.000 claims 1
- 238000004458 analytical method Methods 0.000 description 5
- 241001415846 Procellariidae Species 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 241001076939 Artines Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- JTJMJGYZQZDUJJ-UHFFFAOYSA-N phencyclidine Chemical compound C1CCCCN1C1(C=2C=CC=CC=2)CCCCC1 JTJMJGYZQZDUJJ-UHFFFAOYSA-N 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种三维岩石力学参数获取方法和***,该方法包括:根据已知井的声波测井资料获得横波时差和纵波时差,并确定横纵波时差关系,并据此确定横波波速和纵波波速;根据密度测井数据和纵波波速,确定岩石密度和纵波波速间的密度波速关系;根据纵波波速对地震速度场数据进行校正;根据校正后的地震速度场数据和岩石密度计算三维岩石力学参数,得到动态三维岩石力学参数值;将得到的动态三维岩石力学参数值转化为对应的静态值。通过该方法所得的岩石力学参数综合考虑了测井数据和地震数据,适应于三维地质体的空间非连续性,适用性更强。
Description
技术领域
本发明属于地球物理勘探领域,具体是涉及一种三维岩石力学参数获取方法和***。
背景技术
准确求取岩石力学参数对于制定完善的钻井、完井与油气开发方案和技术措施,降低作业风险具有重要的作用。岩石力学参数有动静态两种值,静态值是通过室内岩芯样品试验直接求得;动态值是根据地球物理测井资料推求计算出来的。室内实验可参考性强,然而由于岩芯样品的获取有限,不能得到整个目的层的相关参数。声波测井具有测量的连续性和较高分辨率的优点,可以精确地确定探测范围内的物性参数和地层参数,但只能得出井点附近岩体的力学参数。
人们通常根据地球物理测井资料推求岩石力学参数,并根据室内试验结果进行校正,用井位附近岩体的参数代替一定区块的参数来进行相关计算。而这种简化忽略了岩体岩性和物理特性的空间非连续性,往往导致随后的相关计算中误差增加,甚至在岩性突变位置的计算结果完全失真。
发明内容
针对现有技术中存在的问题,本发明提供一种三维岩石力学参数获取方法和***,用以克服现有技术中获得的三维岩石力学参数准确性较差的缺陷。
本发明提供了一种三维岩石力学参数获取方法,包括:
根据目标工区已知井的声波测井资料获得横波时差和纵波时差,并确定所述横波时差和所述纵波时差间的横纵波时差关系;
根据所述横纵波时差关系,确定与所述横波时差对应的横波波速,并确定与所述纵波时差对应的纵波波速;
根据所述已知井的密度测井数据和所述纵波波速,确定岩石密度和所述纵波波速间的密度波速关系,并根据所述密度波速关系确定与所述纵波波速分别对应的岩石密度;
根据所述纵波波速对所述目标工区的地震速度场数据进行校正;
根据校正后的地震速度场数据和所述岩石密度计算三维岩石力学参数,得到动态三维岩石力学参数值;
建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系,并将得到的所述动态三维岩石力学参数值转化为对应的静态值。
本发明提供了一种三维岩石力学参数获取***,包括:
第一确定模块,用于根据目标工区已知井的声波测井资料获得横波时差和纵波时差,并确定所述横波时差和所述纵波时差间的横纵波时差关系;
第二确定模块,用于根据所述横纵波时差关系,确定与所述横波时差对应的横波波速,并确定与所述纵波时差对应的纵波波速;
第三确定模块,用于根据所述已知井的密度测井数据和所述纵波波速,确定岩石密度和所述纵波波速间的密度波速关系,并根据所述密度波速关系确定与所述纵波波速分别对应的岩石密度;
校正模块,用于根据所述纵波波速对所述目标工区的地震速度场数据进行校正;
计算模块,用于根据校正后的地震速度场数据和所述岩石密度计算三维岩石力学参数,得到动态三维岩石力学参数值;
转化模块,用于建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系,并将得到的所述动态三维岩石力学参数值转化为对应的静态值。
本发明提供的三维岩石力学参数获取方法和***,根据目标工区已知井的声波测井资料获得横波时差和纵波时差,进而确定与横波时差对应的横波波速,确定与纵波时差对应的纵波波速,进而,根据已知井的密度测井数据和纵波波速,确定岩石密度。为保证结果的准确性,根据纵波波速对目标工区的地震速度场数据进行校正,使得能够根据校正后的地震速度场数据和之前获得的岩石密度计算三维岩石力学参数,得到动态三维岩石力学参数值。为了便于工程分析应用,将得到的所述动态三维岩石力学参数值转化为对应的静态值。通过该方法所得的岩石力学参数综合考虑了测井数据和地震数据,适应于三维地质体的空间非连续性,适用性更强。
附图说明
图1为本发明实施例提供的三维岩石力学参数获取方法的流程图;
图2是本发明实施例中密度与纵波波速回归图;
图3a、图3b是本发明实施例中动静态弹性模量回归直线图和动静态泊松比回归直线图;
图4a、图4b、图4c和图4d是本发明实施例中目标工区4个层面的弹性模量分布云图;
图5a、图5b、图5c和图5d是本发明实施例中目标工区4个层面的泊松比分布云图;
图6为本发明实施例提供的三维岩石力学参数获取***的结构示意图。
具体实施方式
图1为本发明实施例提供的三维岩石力学参数获取方法的流程图,如图1所示,本实施例提供的所述方法具体包括:
步骤101、根据目标工区已知井的声波测井资料获得横波时差和纵波时差,并确定所述横波时差和所述纵波时差间的横纵波时差关系;
具体地,若目标工区有偶极声波测井资料,则直接从偶极声波测井资料中获得横波时差;若目标有声波全波列测井资料,则从声波全波列测井资料中分离出横波时差。而一般测井资料中都会包含纵波时差信息,可以从中获得纵波时差。
以对xx区实地应用为例进行说明。xx区具有进行过纵波和横波测井的典型井的测井资料,通过统计分析,拟合纵横波时差得到公式(1):
Δts=1.777Δtp+0.546(1)
值得说明的是,经过长时间的经验积累,已经获得了针对某些储层的经验公式。比如,若目标工区为中砂岩储层时,根据公式(2)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为:
其中,Δts为横波时差,Δtp为纵波时差,ρb为中砂岩岩石密度;
或者,若目标工区为泥岩储层时,根据公式(3)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为:
Δts/Δtp=B-0.8(ρsh-ρshmin)/(ρshmax-ρshmin)(3)
其中,B为经验常数,ρsh为泥岩岩石密度,当ρsh≤2.2g/cm3时,B=2.5;当ρsh≥2.65g/cm3时,B=1.7。ρshmin常取2.2g/cm3,ρshmax常取2.65g/cm3。
步骤102、根据所述横纵波时差关系,确定与所述横波时差对应的横波波速,并确定与所述纵波时差对应的纵波波速;
仍以上述xx工区为例,在得到横波时差和纵波时差间的对应关系后,根据统计分析结果,根据公式(4)确定与所述横波时差对应的横波波速,并确定与所述纵波时差对应的纵波波速:
Vs=0.3046×106/Δts
Vp=0.3046×106/Δtp(4)
其中,Vs和Vp分别为横波波速和纵波波速。
值得说明的是,也可以根据经验公式:来估算横波波速。
步骤103、根据所述已知井的密度测井数据和所述纵波波速,确定岩石密度和所述纵波波速间的密度波速关系,并根据所述密度波速关系确定与所述纵波波速分别对应的岩石密度;
由于要进行岩石力学参数即泊松比和弹性模量的计算,在计算这些参数时,除了需要用到上述横波波速和纵波波速外,还需要用到岩石密度。本实施例中,可以根据已知井的密度测井数据和纵波波速间的对应关系来估计获得其他未知井的岩石密度,从而获得与目标工区中每个纵波波速对应的岩石密度。
具体地,可以采用指数模型对密度测井数据和纵波波速进行统计分析,确定岩石密度与纵波波速间的密度波速关系。该指数模型的方式与现有技术中相似,不赘述。获得的岩石密度和纵波波速的密度波速关系回归图如图2所示,图2是本发明实施例中密度与纵波波速回归图。
步骤104、根据所述纵波波速对所述目标工区的地震速度场数据进行校正;
为了能够利用根据地震谱资料得到的地震速度场数据计算三维岩石参数,需要解决两个问题:①地震速度场数据与声波测井数据相比,准确性较差并存在概率误差,不能直接用于计算;②横波波速和岩石密度是在纵波波速基础上建立的,不能直接用于地震速度场。因此有必要根据声波测井得到的纵波波速进行地震速度场数据的校正。而地震速度场数据离散性较大,很难进行一般的数据回归分析校正其误差。为了解决这个问题,本实施例提出采用距离加权反比插值法对所述目标工区的地震速度场数据进行校正。
具体地,设空间待插值点即地震速度场中的某速度为P(xp,yp,zp),P点邻域内有已知散乱点即纵波波速Qi(xi,yi,zi),i=1,2,…,n,利用距离加权反比法对P点的属性值Zp进行插值。待插值点P(xp,yp,zp)的属性值是待插值点邻域内已知散乱点属性值的加权平均,权的大小与待插值点与邻域内散乱点之间的距离有关,是距离的k(k>0,一般取2或3)次方的倒数,即:
其中,di为待插值点与其邻域内第i个点之间的距离。
将待插值点与确定的对应的权重值相乘,即得到校正后的速度。
步骤105、根据校正后的地震速度场数据和所述岩石密度计算三维岩石力学参数,得到动态三维岩石力学参数值;
本实施例中的岩石力学参数包括泊松比和弹性模量。根据校正后的速度场数据得到的岩石力学参数为动态值。具体地:
根据公式(5)计算得到动态泊松比和动态弹性模量:
其中,μd为动态泊松比,Ed为动态弹性模量,ρ为岩石密度,V'p和V's分别为校正后的纵波波速和横波波速。
步骤106、建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系,并将得到的所述动态三维岩石力学参数值转化为对应的静态值。
岩石力学参数的静态值和动态值存在着一定的差值,利用声波测井和地震速度场得到的参数是动态值,不能直接用于工程分析中。为此,在三轴应力下对砂、泥岩等岩芯进行岩石力学参数的动静态同步测试,并进行动静态参数的线形回归。通常情况下根据室内岩芯样品三轴试验得到的静态值和相应深度段的计算动态值进行统计分析,获得动静态值的回归关系,如图3a和3b所示。结果表明:岩石的动静态杨氏模量之间有较好的相关性,而动静态泊松比之间的关系不明显。
具体地,上述举例中,对所述动态三维岩石力学参数值和所述静态三维岩石力学参数值进行统计回归分析,得到回归公式(6):
其中,μs为静态泊松比,Es为静态弹性模量;
根据公式(6)将得到的所述动态三维岩石力学参数值转化为对应的静态值。
可选的,在步骤106之后,为了便于工程分析,还包括如下步骤:
三维展布所述静态三维岩石力学参数值,并根据所述目标工区的地质特征信息确定所述静态三维岩石力学参数值是否合理;
若不合理,则调整校正参数,重新根据所述纵波波速对所述目标工区的地震速度场数据进行校正。
本实施例中计算得到的三维岩石力学参数可以采用Surfer、Petrel等绘图软件进行三维展布,方便对比分析岩石力学参数与工程地质力学特性之间的关系。
仍以上述xx工区为例,将计算结果导入petrel地质构造图中,可以绘出每个层面的岩石力学参数。其中,图4a、图4b、图4c和图4d分别给出了该区巴一顶、巴二顶、巴三顶和巴三底4个层面的弹性模量分布云图。由图可知,区域横向上整体表现出南低北高,东西两侧低中间高的趋势。其中北部地区弹性模量较高,达到55GPa以上;在东西两个构造鞍部附进出现了低值区,约为40GPa;井区所在的构造高部位弹性模量约为50GPa。垂向上比较巴一、巴二和巴三层段,随着深度的增加,弹性模量稍有变化,但不很明显。图5a、图5b、图5c和图5d分别给出了该区巴一顶、巴二顶、巴三顶和巴三底4个层面的泊松比分布云图。由图可知,区域横向上整体表现出南低北高,东西两侧低中间高的趋势。最高值达到0.24左右,集中分布于克深2区块北部边界断层附近;最低值约为0.23,零星分布与东西两侧的构造鞍部;井区所在的构造高部位泊松比约在0.235。垂向上比较巴一、巴二和巴三层段,随着深度的增加,泊松比稍有变化,但不很明显。
本实施例中,根据目标工区已知井的声波测井资料获得横波时差和纵波时差,进而确定与横波时差对应的横波波速,确定与纵波时差对应的纵波波速,进而,根据已知井的密度测井数据和纵波波速,确定岩石密度。为保证结果的准确性,根据纵波波速对目标工区的地震速度场数据进行校正,使得能够根据校正后的地震速度场数据和之前获得的岩石密度计算三维岩石力学参数,得到动态三维岩石力学参数值。为了便于工程分析应用,将得到的所述动态三维岩石力学参数值转化为对应的静态值。通过该方法所得的岩石力学参数综合考虑了测井数据和地震数据,适应于三维地质体的空间非连续性,适用性更强。
图6为本发明实施例提供的三维岩石力学参数获取***的结构示意图,如图6所示,该***包括:
第一确定模块11,用于根据目标工区已知井的声波测井资料获得横波时差和纵波时差,并确定所述横波时差和所述纵波时差间的横纵波时差关系;
第二确定模块12,用于根据所述横纵波时差关系,确定与所述横波时差对应的横波波速,并确定与所述纵波时差对应的纵波波速;
第三确定模块13,用于根据所述已知井的密度测井数据和所述纵波波速,确定岩石密度和所述纵波波速间的密度波速关系,并根据所述密度波速关系确定与所述纵波波速分别对应的岩石密度;
校正模块14,用于根据所述纵波波速对所述目标工区的地震速度场数据进行校正;
计算模块15,用于根据校正后的地震速度场数据和所述岩石密度计算三维岩石力学参数,得到动态三维岩石力学参数值;
转化模块16,用于建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系,并将得到的所述动态三维岩石力学参数值转化为对应的静态值。
进一步地,所述***还包括:
展示模块21,用于三维展布所述静态三维岩石力学参数值,并根据所述目标工区的地质特征信息确定所述静态三维岩石力学参数值是否合理;
调整模块22,用于若所述静态三维岩石力学参数值不合理,则调整校正参数,重新根据所述纵波波速对所述目标工区的地震速度场数据进行校正。
进一步地,所述第一确定模块11具体用于:
若所述声波测井资料中有偶极声波测井资料,则直接从所述偶极声波测井资料中获得所述横波时差;
若所述声波测井资料中有声波全波列测井资料,则从所述声波全波列测井资料中分离出所述横波时差;
从所述声波测井资料中获得所述纵波时差。
进一步地,所述第一确定模块11还用于:
对所述横波时差和所述纵波时差进行统计分析,得到所述横纵波时差关系为:
Δts=1.777Δtp+0.546(1)
或者,若所述目标工区为中砂岩储层时,根据公式(2)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为:
其中,Δts为横波时差,Δtp为纵波时差,ρb为中砂岩岩石密度;
或者,若所述目标工区为泥岩储层时,根据公式(3)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为:
Δts/Δtp=B-0.8(ρsh-ρshmin)/(ρshmax-ρshmin)(3)
其中,B为经验常数,ρsh为泥岩岩石密度,当ρsh≤2.2g/cm3时,B=2.5;当ρsh≥2.65g/cm3时,B=1.7。ρshmin常取2.2g/cm3,ρshmax常取2.65g/cm3。
进一步地,所述第二确定模块12具体用于:
根据公式(4)确定与所述横波时差对应的横波波速,并确定与所述纵波时差对应的纵波波速:
Vs=0.3046×106/Δts
Vp=0.3046×106/Δtp(4)
其中,Vs和Vp分别为横波波速和纵波波速。
进一步地,所述第三确定模块13具体用于:
采用指数模型对所述密度测井数据和所述纵波波速进行统计分析,确定岩石密度与所述纵波波速间的密度波速关系。
进一步地,所述校正模块14具体用于:
根据所述纵波波速,采用距离加权反比插值法对所述目标工区的地震速度场数据进行校正。
进一步地,所述岩石力学参数包括泊松比和弹性模量,所述计算模块15具体用于:
根据公式(5)计算得到动态泊松比和动态弹性模量:
其中,μd为动态泊松比,Ed为动态弹性模量,ρ为岩石密度,V'p和V's分别为校正后的纵波波速和横波波速。
进一步地,所述转化模块16具体用于:
对所述动态三维岩石力学参数值和所述静态三维岩石力学参数值进行统计回归分析,得到回归公式(6):
其中,μs为静态泊松比,Es为静态弹性模量;
根据公式(6)将得到的所述动态三维岩石力学参数值转化为对应的静态值。
本实施例的***可以用于执行图1所示方法实施例的技术方案,其实现原理和技术效果类似,此处不再赘述。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (18)
1.一种三维岩石力学参数获取方法,其特征在于,包括:
根据目标工区已知井的声波测井资料获得横波时差和纵波时差,并确定所述横波时差和所述纵波时差间的横纵波时差关系;
根据所述横纵波时差关系,确定与所述横波时差对应的横波波速,并确定与所述纵波时差对应的纵波波速;
根据所述已知井的密度测井数据和所述纵波波速,确定岩石密度和所述纵波波速间的密度波速关系,并根据所述密度波速关系确定与所述纵波波速分别对应的岩石密度;
根据所述纵波波速对所述目标工区的地震速度场数据进行校正;
根据校正后的地震速度场数据和所述岩石密度计算三维岩石力学参数,得到动态三维岩石力学参数值;
建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系,并将得到的所述动态三维岩石力学参数值转化为对应的静态值。
2.根据权利要求1所述的方法,其特征在于,所述建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系,并将得到的所述动态三维岩石力学参数值转化为对应的静态三维岩石力学参数值之后,还包括:
三维展布所述静态三维岩石力学参数值,并根据所述目标工区的地质特征信息确定所述静态三维岩石力学参数值是否合理;
若不合理,则调整校正参数,重新根据所述纵波波速对所述目标工区的地震速度场数据进行校正。
3.根据权利要求1所述的方法,其特征在于,所述根据目标工区已知井的声波测井资料获得横波时差和纵波时差,包括:
若所述声波测井资料中有偶极声波测井资料,则直接从所述偶极声波测井资料中获得所述横波时差;
若所述声波测井资料中有声波全波列测井资料,则从所述声波全波列测井资料中分离出所述横波时差;
从所述声波测井资料中获得所述纵波时差。
4.根据权利要求3所述的方法,其特征在于,所述确定所述横波时差和所述纵波时差间的横纵波时差关系,包括:
对所述横波时差和所述纵波时差进行统计分析,得到所述横纵波时差关系为:
Δts=1.777Δtp+0.546(1)
或者,若所述目标工区为中砂岩储层时,根据公式(2)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为:
其中,Δts为横波时差,Δtp为纵波时差,ρb为中砂岩岩石密度;
或者,若所述目标工区为泥岩储层时,根据公式(3)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为:
Δts/Δtp=B-0.8(ρsh-ρshmin)/(ρshmax-ρshmin)(3)
其中,B为经验常数,ρsh为泥岩岩石密度,当ρsh≤2.2g/cm3时,B=2.5;
当ρsh≥2.65g/cm3时,B=1.7。ρshmin常取2.2g/cm3,ρshmax常取2.65g/cm3。
5.根据权利要求4所述的方法,其特征在于,所述根据所述横纵波时差关系,确定与所述横波时差对应的横波波速,并确定与所述纵波时差对应的纵波波速,包括:
根据公式(4)确定与所述横波时差对应的横波波速,并确定与所述纵波时差对应的纵波波速:
Vs=0.3046×106/Δts
Vp=0.3046×106/Δtp(4)
其中,Vs和Vp分别为横波波速和纵波波速。
6.根据权利要求1至5中任一项所述的方法,其特征在于,所述根据所述已知井的密度测井数据和所述纵波波速,确定岩石密度和所述纵波波速间的密度波速关系,包括:
采用指数模型对所述密度测井数据和所述纵波波速进行统计分析,确定岩石密度与所述纵波波速间的密度波速关系。
7.根据权利要求6所述的方法,其特征在于,所述根据所述纵波波速对所述目标工区的地震速度场数据进行校正,包括:
根据所述纵波波速,采用距离加权反比插值法对所述目标工区的地震速度场数据进行校正。
8.根据权利要求7所述的方法,其特征在于,所述根据校正后的地震速度场数据和所述岩石密度计算三维岩石力学参数,得到动态三维岩石力学参数值,包括:
所述岩石力学参数包括泊松比和弹性模量,根据公式(5)计算得到动态泊松比和动态弹性模量:
其中,μd为动态泊松比,Ed为动态弹性模量,ρ为岩石密度,V'p和V's分别为校正后的纵波波速和横波波速。
9.根据权利要求8所述的方法,其特征在于,所述建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系,并将得到的所述动态三维岩石力学参数值转化为对应的静态值,包括:
对所述动态三维岩石力学参数值和所述静态三维岩石力学参数值进行统计回归分析,得到回归公式(6):
其中,μs为静态泊松比,Es为静态弹性模量;
根据公式(6)将得到的所述动态三维岩石力学参数值转化为对应的静态值。
10.一种三维岩石力学参数获取***,其特征在于,包括:
第一确定模块,用于根据目标工区已知井的声波测井资料获得横波时差和纵波时差,并确定所述横波时差和所述纵波时差间的横纵波时差关系;
第二确定模块,用于根据所述横纵波时差关系,确定与所述横波时差对应的横波波速,并确定与所述纵波时差对应的纵波波速;
第三确定模块,用于根据所述已知井的密度测井数据和所述纵波波速,确定岩石密度和所述纵波波速间的密度波速关系,并根据所述密度波速关系确定与所述纵波波速分别对应的岩石密度;
校正模块,用于根据所述纵波波速对所述目标工区的地震速度场数据进行校正;
计算模块,用于根据校正后的地震速度场数据和所述岩石密度计算三维岩石力学参数,得到动态三维岩石力学参数值;
转化模块,用于建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系,并将得到的所述动态三维岩石力学参数值转化为对应的静态值。
11.根据权利要求10所述的***,其特征在于,还包括:
展示模块,用于三维展布所述静态三维岩石力学参数值,并根据所述目标工区的地质特征信息确定所述静态三维岩石力学参数值是否合理;
调整模块,用于若所述静态三维岩石力学参数值不合理,则调整校正参数,重新根据所述纵波波速对所述目标工区的地震速度场数据进行校正。
12.根据权利要求10所述的***,其特征在于,所述第一确定模块具体用于:
若所述声波测井资料中有偶极声波测井资料,则直接从所述偶极声波测井资料中获得所述横波时差;
若所述声波测井资料中有声波全波列测井资料,则从所述声波全波列测井资料中分离出所述横波时差;
从所述声波测井资料中获得所述纵波时差。
13.根据权利要求12所述的***,其特征在于,所述第一确定模块还用于:
对所述横波时差和所述纵波时差进行统计分析,得到所述横纵波时差关系为:
Δts=1.777Δtp+0.546(1)
或者,若所述目标工区为中砂岩储层时,根据公式(2)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为:
其中,Δts为横波时差,Δtp为纵波时差,ρb为中砂岩岩石密度;
或者,若所述目标工区为泥岩储层时,根据公式(3)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为:
Δts/Δtp=B-0.8(ρsh-ρshmin)/(ρshmax-ρshmin)(3)
其中,B为经验常数,ρsh为泥岩岩石密度,当ρsh≤2.2g/cm3时,B=2.5;
当ρsh≥2.65g/cm3时,B=1.7。ρshmin常取2.2g/cm3,ρshmax常取2.65g/cm3。
14.根据权利要求13所述的***,其特征在于,所述第二确定模块具体用于:
根据公式(4)确定与所述横波时差对应的横波波速,并确定与所述纵波时差对应的纵波波速:
Vs=0.3046×106/Δts
Vp=0.3046×106/Δtp(4)
其中,Vs和Vp分别为横波波速和纵波波速。
15.根据权利要求10至14中任一项所述的***,其特征在于,所述第三确定模块具体用于:
采用指数模型对所述密度测井数据和所述纵波波速进行统计分析,确定岩石密度与所述纵波波速间的密度波速关系。
16.根据权利要求15所述的***,其特征在于,所述校正模块具体用于:
根据所述纵波波速,采用距离加权反比插值法对所述目标工区的地震速度场数据进行校正。
17.根据权利要求16所述的***,其特征在于,所述岩石力学参数包括泊松比和弹性模量,所述计算模块具体用于:
根据公式(5)计算得到动态泊松比和动态弹性模量:
其中,μd为动态泊松比,Ed为动态弹性模量,ρ为岩石密度,V'p和V's分别为校正后的纵波波速和横波波速。
18.根据权利要求17所述的***,其特征在于,所述转化模块具体用于:
对所述动态三维岩石力学参数值和所述静态三维岩石力学参数值进行统计回归分析,得到回归公式(6):
其中,μs为静态泊松比,Es为静态弹性模量;
根据公式(6)将得到的所述动态三维岩石力学参数值转化为对应的静态值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410641421.2A CN105588883B (zh) | 2014-11-13 | 2014-11-13 | 三维岩石力学参数获取方法和*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410641421.2A CN105588883B (zh) | 2014-11-13 | 2014-11-13 | 三维岩石力学参数获取方法和*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105588883A true CN105588883A (zh) | 2016-05-18 |
CN105588883B CN105588883B (zh) | 2018-05-08 |
Family
ID=55928625
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410641421.2A Active CN105588883B (zh) | 2014-11-13 | 2014-11-13 | 三维岩石力学参数获取方法和*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105588883B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110261485A (zh) * | 2019-07-26 | 2019-09-20 | 北京领示科技有限公司 | 一种超声波测量材料内部各处弹性模量及泊松比的方法 |
CN111060416A (zh) * | 2019-12-23 | 2020-04-24 | 中国石油大学(北京) | 一种地表露头与岩心岩石力学层划分方法及*** |
WO2020098203A1 (zh) * | 2018-11-13 | 2020-05-22 | 浙江大学 | 一种用于超重力地质构造物理模拟实验的数据处理方法 |
CN112394397A (zh) * | 2019-08-13 | 2021-02-23 | 中国石油化工股份有限公司 | 一种页岩气藏三维岩石力学参数场建模方法 |
CN113777668A (zh) * | 2020-06-10 | 2021-12-10 | 中国石油化工股份有限公司 | 用于砂泥岩互层致密气藏的地应力计算方法及装置 |
CN113775329A (zh) * | 2020-06-09 | 2021-12-10 | 中国石油化工股份有限公司 | 一种获取深层碳酸盐岩岩石力学参数的方法及*** |
CN114112651A (zh) * | 2020-08-27 | 2022-03-01 | 中国石油化工股份有限公司 | 一种用于人造岩心的岩石动静态力学参数转换方法及*** |
CN115407045A (zh) * | 2022-08-02 | 2022-11-29 | 西南石油大学 | 岩石力学参数评价模型构建方法及岩石力学特性评价方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1719979A1 (ru) * | 1990-01-17 | 1992-03-15 | Ленинградский горный институт им.Г.В.Плеханова | Способ определени физико-механических параметров в плоско-параллельных объектах |
CN101231346A (zh) * | 2008-02-26 | 2008-07-30 | 中国矿业大学(北京) | 应用地震波速度预测煤、岩体物理力学参数方法 |
KR20090077379A (ko) * | 2008-01-11 | 2009-07-15 | 한국지질자원연구원 | 탄성파 속도 측정 홀더 및 탄성파 속도 측정기 |
CN103267678A (zh) * | 2013-04-28 | 2013-08-28 | 北京大学 | 一种动静态岩石力学参数同步测量方法及装置 |
CN103323352A (zh) * | 2013-06-07 | 2013-09-25 | 中国石油天然气股份有限公司 | 天然气水合物沉积物动三轴力学-声学-电学同步测试的实验装置及方法 |
-
2014
- 2014-11-13 CN CN201410641421.2A patent/CN105588883B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1719979A1 (ru) * | 1990-01-17 | 1992-03-15 | Ленинградский горный институт им.Г.В.Плеханова | Способ определени физико-механических параметров в плоско-параллельных объектах |
KR20090077379A (ko) * | 2008-01-11 | 2009-07-15 | 한국지질자원연구원 | 탄성파 속도 측정 홀더 및 탄성파 속도 측정기 |
CN101231346A (zh) * | 2008-02-26 | 2008-07-30 | 中国矿业大学(北京) | 应用地震波速度预测煤、岩体物理力学参数方法 |
CN103267678A (zh) * | 2013-04-28 | 2013-08-28 | 北京大学 | 一种动静态岩石力学参数同步测量方法及装置 |
CN103323352A (zh) * | 2013-06-07 | 2013-09-25 | 中国石油天然气股份有限公司 | 天然气水合物沉积物动三轴力学-声学-电学同步测试的实验装置及方法 |
Non-Patent Citations (2)
Title |
---|
MOHAMMED S.AMEEN ET.AL: "Predicting rock mechanical properties of carbonates from wireline logs (A case study:Arab-D reservoir,Ghawar field,Saudi Arabia)", 《MARINE AND PETROLEUM GEOLOGY》 * |
王贤君 等: "井间岩石力学参数预测", 《大庆石油地质与开发》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020098203A1 (zh) * | 2018-11-13 | 2020-05-22 | 浙江大学 | 一种用于超重力地质构造物理模拟实验的数据处理方法 |
CN110261485A (zh) * | 2019-07-26 | 2019-09-20 | 北京领示科技有限公司 | 一种超声波测量材料内部各处弹性模量及泊松比的方法 |
CN112394397A (zh) * | 2019-08-13 | 2021-02-23 | 中国石油化工股份有限公司 | 一种页岩气藏三维岩石力学参数场建模方法 |
CN111060416A (zh) * | 2019-12-23 | 2020-04-24 | 中国石油大学(北京) | 一种地表露头与岩心岩石力学层划分方法及*** |
CN113775329A (zh) * | 2020-06-09 | 2021-12-10 | 中国石油化工股份有限公司 | 一种获取深层碳酸盐岩岩石力学参数的方法及*** |
CN113777668A (zh) * | 2020-06-10 | 2021-12-10 | 中国石油化工股份有限公司 | 用于砂泥岩互层致密气藏的地应力计算方法及装置 |
CN114112651A (zh) * | 2020-08-27 | 2022-03-01 | 中国石油化工股份有限公司 | 一种用于人造岩心的岩石动静态力学参数转换方法及*** |
CN115407045A (zh) * | 2022-08-02 | 2022-11-29 | 西南石油大学 | 岩石力学参数评价模型构建方法及岩石力学特性评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105588883B (zh) | 2018-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105588883A (zh) | 三维岩石力学参数获取方法和*** | |
CN109490965B (zh) | 一种定量评价地层非均匀性的方法及装置 | |
CN104502997B (zh) | 一种利用裂缝密度曲线预测裂缝密度体的方法 | |
CN108072902B (zh) | 一种碳酸盐岩岩石物理模型的建模方法 | |
CN102096107B (zh) | 一种根据声波时差和密度反演孔隙扁度进行储层渗透性评价的方法 | |
CN105445791A (zh) | 一种基于多种地震属性的地层孔隙压力预测方法 | |
CN109655903B (zh) | 页岩层横波速度预测方法及*** | |
CN104677734A (zh) | 一种多尺度信息融合的致密地层脆性指数测量方法 | |
CN104181585A (zh) | 一种地球物理勘探横波估算方法和*** | |
CN107329171A (zh) | 深度域储层地震反演方法及装置 | |
CN103954999A (zh) | 一种适用于低孔隙度砂泥岩地层的横波速度预测新方法 | |
CN109799540B (zh) | 基于地质信息约束下的火山岩型铀矿床磁化率反演方法 | |
CN105738952A (zh) | 一种水平井区储层岩石相建模方法 | |
CN105549088A (zh) | 裂缝性致密砂岩中气层的识别方法和装置 | |
CN105370270B (zh) | 由偶极声波纵、横波时差确定页岩气储层含气饱和度的方法 | |
RU2485553C1 (ru) | Способ оценки трещинной пористости по данным скважинной сейсморазведки | |
CN103628866B (zh) | 获取地层的地层压力系数的方法 | |
CN107728205A (zh) | 一种地层压力预测方法 | |
CN104749630A (zh) | 构建微地震监测速度模型的方法 | |
US20220050224A1 (en) | Evaluating anisotropic effective permeability in rock formations having natural fracture networks | |
CN107728208B (zh) | 基于地层速度异常预测孔隙压力的方法 | |
CN110927787B (zh) | 基于虚拟岩性的横波速度计算方法及*** | |
CN111308558A (zh) | 页岩气水平井纵波时差校正方法 | |
CN103424774B (zh) | 一种地震油气藏成像方法及装置 | |
CN106353837B (zh) | 基于加权gardner公式的密度预测方法 |
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 |