CN105588883A - 三维岩石力学参数获取方法和系统 - Google Patents

Abstract

本发明提供一种三维岩石力学参数获取方法和系统，该方法包括：根据已知井的声波测井资料获得横波时差和纵波时差，并确定横纵波时差关系，并据此确定横波波速和纵波波速；根据密度测井数据和纵波波速，确定岩石密度和纵波波速间的密度波速关系；根据纵波波速对地震速度场数据进行校正；根据校正后的地震速度场数据和岩石密度计算三维岩石力学参数，得到动态三维岩石力学参数值；将得到的动态三维岩石力学参数值转化为对应的静态值。通过该方法所得的岩石力学参数综合考虑了测井数据和地震数据，适应于三维地质体的空间非连续性，适用性更强。

Description

三维岩石力学参数获取方法和系统

技术领域

本发明属于地球物理勘探领域，具体是涉及一种三维岩石力学参数获取方法和系统。

背景技术

准确求取岩石力学参数对于制定完善的钻井、完井与油气开发方案和技术措施，降低作业风险具有重要的作用。岩石力学参数有动静态两种值，静态值是通过室内岩芯样品试验直接求得；动态值是根据地球物理测井资料推求计算出来的。室内实验可参考性强，然而由于岩芯样品的获取有限，不能得到整个目的层的相关参数。声波测井具有测量的连续性和较高分辨率的优点，可以精确地确定探测范围内的物性参数和地层参数，但只能得出井点附近岩体的力学参数。

人们通常根据地球物理测井资料推求岩石力学参数，并根据室内试验结果进行校正，用井位附近岩体的参数代替一定区块的参数来进行相关计算。而这种简化忽略了岩体岩性和物理特性的空间非连续性，往往导致随后的相关计算中误差增加，甚至在岩性突变位置的计算结果完全失真。

发明内容

针对现有技术中存在的问题，本发明提供一种三维岩石力学参数获取方法和系统，用以克服现有技术中获得的三维岩石力学参数准确性较差的缺陷。

本发明提供了一种三维岩石力学参数获取方法，包括：

根据目标工区已知井的声波测井资料获得横波时差和纵波时差，并确定所述横波时差和所述纵波时差间的横纵波时差关系；

根据所述横纵波时差关系，确定与所述横波时差对应的横波波速，并确定与所述纵波时差对应的纵波波速；

根据所述已知井的密度测井数据和所述纵波波速，确定岩石密度和所述纵波波速间的密度波速关系，并根据所述密度波速关系确定与所述纵波波速分别对应的岩石密度；

根据所述纵波波速对所述目标工区的地震速度场数据进行校正；

根据校正后的地震速度场数据和所述岩石密度计算三维岩石力学参数，得到动态三维岩石力学参数值；

建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系，并将得到的所述动态三维岩石力学参数值转化为对应的静态值。

本发明提供了一种三维岩石力学参数获取系统，包括：

第一确定模块，用于根据目标工区已知井的声波测井资料获得横波时差和纵波时差，并确定所述横波时差和所述纵波时差间的横纵波时差关系；

第二确定模块，用于根据所述横纵波时差关系，确定与所述横波时差对应的横波波速，并确定与所述纵波时差对应的纵波波速；

第三确定模块，用于根据所述已知井的密度测井数据和所述纵波波速，确定岩石密度和所述纵波波速间的密度波速关系，并根据所述密度波速关系确定与所述纵波波速分别对应的岩石密度；

校正模块，用于根据所述纵波波速对所述目标工区的地震速度场数据进行校正；

计算模块，用于根据校正后的地震速度场数据和所述岩石密度计算三维岩石力学参数，得到动态三维岩石力学参数值；

转化模块，用于建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系，并将得到的所述动态三维岩石力学参数值转化为对应的静态值。

本发明提供的三维岩石力学参数获取方法和系统，根据目标工区已知井的声波测井资料获得横波时差和纵波时差，进而确定与横波时差对应的横波波速，确定与纵波时差对应的纵波波速，进而，根据已知井的密度测井数据和纵波波速，确定岩石密度。为保证结果的准确性，根据纵波波速对目标工区的地震速度场数据进行校正，使得能够根据校正后的地震速度场数据和之前获得的岩石密度计算三维岩石力学参数，得到动态三维岩石力学参数值。为了便于工程分析应用，将得到的所述动态三维岩石力学参数值转化为对应的静态值。通过该方法所得的岩石力学参数综合考虑了测井数据和地震数据，适应于三维地质体的空间非连续性，适用性更强。

附图说明

图1为本发明实施例提供的三维岩石力学参数获取方法的流程图；

图2是本发明实施例中密度与纵波波速回归图；

图3a、图3b是本发明实施例中动静态弹性模量回归直线图和动静态泊松比回归直线图；

图4a、图4b、图4c和图4d是本发明实施例中目标工区4个层面的弹性模量分布云图；

图5a、图5b、图5c和图5d是本发明实施例中目标工区4个层面的泊松比分布云图；

图6为本发明实施例提供的三维岩石力学参数获取系统的结构示意图。

具体实施方式

图1为本发明实施例提供的三维岩石力学参数获取方法的流程图，如图1所示，本实施例提供的所述方法具体包括：

步骤101、根据目标工区已知井的声波测井资料获得横波时差和纵波时差，并确定所述横波时差和所述纵波时差间的横纵波时差关系；

具体地，若目标工区有偶极声波测井资料，则直接从偶极声波测井资料中获得横波时差；若目标有声波全波列测井资料，则从声波全波列测井资料中分离出横波时差。而一般测井资料中都会包含纵波时差信息，可以从中获得纵波时差。

以对xx区实地应用为例进行说明。xx区具有进行过纵波和横波测井的典型井的测井资料，通过统计分析，拟合纵横波时差得到公式(1)：

Δt_s＝1.777Δt_p+0.546(¹)

值得说明的是，经过长时间的经验积累，已经获得了针对某些储层的经验公式。比如，若目标工区为中砂岩储层时，根据公式(2)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为：

${\mathrm{\Δt}}_s=\frac{{\mathrm{\Δt}}_p}{{[1-1.15\frac{(1/{\mathrm{\ρ}}_b)+{(1+{\mathrm{\ρ}}_b)}^3}{e^{1/{\mathrm{\ρ}}_b}}]}^{1.5}}---\left(2\right)$

其中，Δt_s为横波时差，Δt_p为纵波时差，ρ_b为中砂岩岩石密度；

或者，若目标工区为泥岩储层时，根据公式(3)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为：

Δt_s/Δt_p＝B-0.8(ρ_sh-ρ_shmin)/(ρ_shmax-ρ_shmin)(3)

其中，B为经验常数，ρ_sh为泥岩岩石密度，当ρ_sh≤2.2g/cm³时，B＝2.5；当ρ_sh≥2.65g/cm³时，B＝1.7。ρ_shmin常取2.2g/cm³，ρ_shmax常取2.65g/cm³。

步骤102、根据所述横纵波时差关系，确定与所述横波时差对应的横波波速，并确定与所述纵波时差对应的纵波波速；

仍以上述xx工区为例，在得到横波时差和纵波时差间的对应关系后，根据统计分析结果，根据公式(4)确定与所述横波时差对应的横波波速，并确定与所述纵波时差对应的纵波波速：

V_s＝0.3046×10⁶/Δt_s

V_p＝0.3046×10⁶/Δt_p(4)

其中，V_s和V_p分别为横波波速和纵波波速。

值得说明的是，也可以根据经验公式：来估算横波波速。

步骤103、根据所述已知井的密度测井数据和所述纵波波速，确定岩石密度和所述纵波波速间的密度波速关系，并根据所述密度波速关系确定与所述纵波波速分别对应的岩石密度；

由于要进行岩石力学参数即泊松比和弹性模量的计算，在计算这些参数时，除了需要用到上述横波波速和纵波波速外，还需要用到岩石密度。本实施例中，可以根据已知井的密度测井数据和纵波波速间的对应关系来估计获得其他未知井的岩石密度，从而获得与目标工区中每个纵波波速对应的岩石密度。

具体地，可以采用指数模型对密度测井数据和纵波波速进行统计分析，确定岩石密度与纵波波速间的密度波速关系。该指数模型的方式与现有技术中相似，不赘述。获得的岩石密度和纵波波速的密度波速关系回归图如图2所示，图2是本发明实施例中密度与纵波波速回归图。

步骤104、根据所述纵波波速对所述目标工区的地震速度场数据进行校正；

为了能够利用根据地震谱资料得到的地震速度场数据计算三维岩石参数，需要解决两个问题：①地震速度场数据与声波测井数据相比，准确性较差并存在概率误差，不能直接用于计算；②横波波速和岩石密度是在纵波波速基础上建立的，不能直接用于地震速度场。因此有必要根据声波测井得到的纵波波速进行地震速度场数据的校正。而地震速度场数据离散性较大，很难进行一般的数据回归分析校正其误差。为了解决这个问题，本实施例提出采用距离加权反比插值法对所述目标工区的地震速度场数据进行校正。

具体地，设空间待插值点即地震速度场中的某速度为P(x_p,y_p,z_p)，P点邻域内有已知散乱点即纵波波速Qi(x_i,y_i,z_i),i＝1,2,…,n，利用距离加权反比法对P点的属性值Zp进行插值。待插值点P(x_p,y_p,z_p)的属性值是待插值点邻域内已知散乱点属性值的加权平均，权的大小与待插值点与邻域内散乱点之间的距离有关，是距离的k(k>0，一般取2或3)次方的倒数，即：

$Z_p=\frac{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}\frac{Z_i}{d_i^k}}{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}\frac{1}{d_i^k}}$

其中，d_i为待插值点与其邻域内第i个点之间的距离。

将待插值点与确定的对应的权重值相乘，即得到校正后的速度。

步骤105、根据校正后的地震速度场数据和所述岩石密度计算三维岩石力学参数，得到动态三维岩石力学参数值；

本实施例中的岩石力学参数包括泊松比和弹性模量。根据校正后的速度场数据得到的岩石力学参数为动态值。具体地：

根据公式(5)计算得到动态泊松比和动态弹性模量：

$\left\{\begin{array}{c}{\mathrm{\μ}}_d=\frac{{V\′}_p^2-{2V\′}_s^2}{2({V\′}_p^2-{V\′}_s^2)}\\ E_d={10}^{-3}{\mathrm{\ρV}\′}_s^2\frac{{3V\′}_p^2-{4V\′}_s^2}{{V\′}_p^2-{V\′}_s^2}\end{array}\right.---\left(5\right)$

其中，μ_d为动态泊松比，E_d为动态弹性模量，ρ为岩石密度，V'_p和V'_s分别为校正后的纵波波速和横波波速。

步骤106、建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系，并将得到的所述动态三维岩石力学参数值转化为对应的静态值。

岩石力学参数的静态值和动态值存在着一定的差值，利用声波测井和地震速度场得到的参数是动态值，不能直接用于工程分析中。为此，在三轴应力下对砂、泥岩等岩芯进行岩石力学参数的动静态同步测试，并进行动静态参数的线形回归。通常情况下根据室内岩芯样品三轴试验得到的静态值和相应深度段的计算动态值进行统计分析，获得动静态值的回归关系，如图3a和3b所示。结果表明：岩石的动静态杨氏模量之间有较好的相关性，而动静态泊松比之间的关系不明显。

具体地，上述举例中，对所述动态三维岩石力学参数值和所述静态三维岩石力学参数值进行统计回归分析，得到回归公式(6)：

$\left\{\begin{array}{c}{\mathrm{\μ}}_s={\mathrm{\μ}}_d\\ E_s=0.7639E_d+4447.2\end{array}---\left(6\right)\right.$

其中，μ_s为静态泊松比，E_s为静态弹性模量；

根据公式(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还用于：

对所述横波时差和所述纵波时差进行统计分析，得到所述横纵波时差关系为：

Δt_s＝1.777Δt_p+0.546(1)

或者，若所述目标工区为中砂岩储层时，根据公式(2)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为：

${\mathrm{\Δt}}_s=\frac{{\mathrm{\Δt}}_p}{{[1-1.15\frac{(1/{\mathrm{\ρ}}_b)+{(1+{\mathrm{\ρ}}_b)}^3}{e^{1/{\mathrm{\ρ}}_b}}]}^{1.5}}---\left(2\right)$

其中，Δt_s为横波时差，Δt_p为纵波时差，ρ_b为中砂岩岩石密度；

或者，若所述目标工区为泥岩储层时，根据公式(3)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为：

Δt_s/Δt_p＝B-0.8(ρ_sh-ρ_shmin)/(ρ_shmax-ρ_shmin)(3)

其中，B为经验常数，ρ_sh为泥岩岩石密度，当ρ_sh≤2.2g/cm³时，B＝2.5；当ρ_sh≥2.65g/cm³时，B＝1.7。ρ_shmin常取2.2g/cm³，ρ_shmax常取2.65g/cm³。

进一步地，所述第二确定模块12具体用于：

根据公式(4)确定与所述横波时差对应的横波波速，并确定与所述纵波时差对应的纵波波速：

V_s＝0.3046×10⁶/Δt_s

V_p＝0.3046×10⁶/Δt_p(4)

其中，V_s和V_p分别为横波波速和纵波波速。

进一步地，所述第三确定模块13具体用于：

采用指数模型对所述密度测井数据和所述纵波波速进行统计分析，确定岩石密度与所述纵波波速间的密度波速关系。

进一步地，所述校正模块14具体用于：

根据所述纵波波速，采用距离加权反比插值法对所述目标工区的地震速度场数据进行校正。

进一步地，所述岩石力学参数包括泊松比和弹性模量，所述计算模块15具体用于：

根据公式(5)计算得到动态泊松比和动态弹性模量：

$\left\{\begin{array}{c}{\mathrm{\μ}}_d=\frac{{V\′}_p^2-{2V\′}_s^2}{2({V\′}_p^2-{V\′}_s^2)}\\ E_d={10}^{-3}{\mathrm{\ρV}\′}_s^2\frac{{3V\′}_p^2-{4V\′}_s^2}{{V\′}_p^2-{V\′}_s^2}\end{array}\right.---\left(5\right)$

其中，μ_d为动态泊松比，E_d为动态弹性模量，ρ为岩石密度，V'_p和V'_s分别为校正后的纵波波速和横波波速。

进一步地，所述转化模块16具体用于：

对所述动态三维岩石力学参数值和所述静态三维岩石力学参数值进行统计回归分析，得到回归公式(6)：

$\left\{\begin{array}{c}{\mathrm{\μ}}_s={\mathrm{\μ}}_d\\ E_s=0.7639E_d+4447.2\end{array}---\left(6\right)\right.$

其中，μ_s为静态泊松比，E_s为静态弹性模量；

根据公式(6)将得到的所述动态三维岩石力学参数值转化为对应的静态值。

本实施例的系统可以用于执行图1所示方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。

最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (18)

1.一种三维岩石力学参数获取方法，其特征在于，包括：

根据目标工区已知井的声波测井资料获得横波时差和纵波时差，并确定所述横波时差和所述纵波时差间的横纵波时差关系；

根据所述横纵波时差关系，确定与所述横波时差对应的横波波速，并确定与所述纵波时差对应的纵波波速；

根据所述已知井的密度测井数据和所述纵波波速，确定岩石密度和所述纵波波速间的密度波速关系，并根据所述密度波速关系确定与所述纵波波速分别对应的岩石密度；

根据所述纵波波速对所述目标工区的地震速度场数据进行校正；

根据校正后的地震速度场数据和所述岩石密度计算三维岩石力学参数，得到动态三维岩石力学参数值；

建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系，并将得到的所述动态三维岩石力学参数值转化为对应的静态值。

2.根据权利要求1所述的方法，其特征在于，所述建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系，并将得到的所述动态三维岩石力学参数值转化为对应的静态三维岩石力学参数值之后，还包括：

三维展布所述静态三维岩石力学参数值，并根据所述目标工区的地质特征信息确定所述静态三维岩石力学参数值是否合理；

若不合理，则调整校正参数，重新根据所述纵波波速对所述目标工区的地震速度场数据进行校正。

3.根据权利要求1所述的方法，其特征在于，所述根据目标工区已知井的声波测井资料获得横波时差和纵波时差，包括：

若所述声波测井资料中有偶极声波测井资料，则直接从所述偶极声波测井资料中获得所述横波时差；

若所述声波测井资料中有声波全波列测井资料，则从所述声波全波列测井资料中分离出所述横波时差；

从所述声波测井资料中获得所述纵波时差。

4.根据权利要求3所述的方法，其特征在于，所述确定所述横波时差和所述纵波时差间的横纵波时差关系，包括：

对所述横波时差和所述纵波时差进行统计分析，得到所述横纵波时差关系为：

Δt_s＝1.777Δt_p+0.546(1)

或者，若所述目标工区为中砂岩储层时，根据公式(2)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为：

${\mathrm{\Δt}}_s=\frac{{\mathrm{\Δt}}_p}{{[1-1.15\frac{(1/{\mathrm{\ρ}}_b)+{(1+{\mathrm{\ρ}}_b)}^3}{e^{1/{\mathrm{\ρ}}_b}}]}^{1.5}}---\left(2\right)$

其中，Δt_s为横波时差，Δt_p为纵波时差，ρ_b为中砂岩岩石密度；

或者，若所述目标工区为泥岩储层时，根据公式(3)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为：

Δt_s/Δt_p＝B-0.8(ρ_sh-ρ_shmin)/(ρ_shmax-ρ_shmin)(3)

其中，B为经验常数，ρ_sh为泥岩岩石密度，当ρ_sh≤2.2g/cm³时，B＝2.5；

当ρ_sh≥2.65g/cm³时，B＝1.7。ρ_shmin常取2.2g/cm³，ρ_shmax常取2.65g/cm³。

5.根据权利要求4所述的方法，其特征在于，所述根据所述横纵波时差关系，确定与所述横波时差对应的横波波速，并确定与所述纵波时差对应的纵波波速，包括：

根据公式(4)确定与所述横波时差对应的横波波速，并确定与所述纵波时差对应的纵波波速：

V_s＝0.3046×10⁶/Δt_s

V_p＝0.3046×10⁶/Δt_p(4)

其中，V_s和V_p分别为横波波速和纵波波速。

6.根据权利要求1至5中任一项所述的方法，其特征在于，所述根据所述已知井的密度测井数据和所述纵波波速，确定岩石密度和所述纵波波速间的密度波速关系，包括：

采用指数模型对所述密度测井数据和所述纵波波速进行统计分析，确定岩石密度与所述纵波波速间的密度波速关系。

7.根据权利要求6所述的方法，其特征在于，所述根据所述纵波波速对所述目标工区的地震速度场数据进行校正，包括：

根据所述纵波波速，采用距离加权反比插值法对所述目标工区的地震速度场数据进行校正。

8.根据权利要求7所述的方法，其特征在于，所述根据校正后的地震速度场数据和所述岩石密度计算三维岩石力学参数，得到动态三维岩石力学参数值，包括：

所述岩石力学参数包括泊松比和弹性模量，根据公式(5)计算得到动态泊松比和动态弹性模量：

$\left\{\begin{array}{c}{\mathrm{\μ}}_d=\frac{V_p^{\′2}-{2V}_s^{\′2}}{2(V_p^{\′2}-V_s^{\′2})}\\ E_d={10}^{-3}{\mathrm{\ρV}}_s^{\′2}\frac{{3V}_p^{\′2}-{4V}_s^{\′2}}{V_p^{\′2}-V_s^{\′2}}\end{array}\right.---\left(5\right)$

其中，μ_d为动态泊松比，E_d为动态弹性模量，ρ为岩石密度，V'_p和V'_s分别为校正后的纵波波速和横波波速。

9.根据权利要求8所述的方法，其特征在于，所述建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系，并将得到的所述动态三维岩石力学参数值转化为对应的静态值，包括：

对所述动态三维岩石力学参数值和所述静态三维岩石力学参数值进行统计回归分析，得到回归公式(6)：

$\left\{\begin{array}{c}{\mathrm{\μ}}_s={\mathrm{\μ}}_d\\ E_s=0.7639E_d+4447.2\end{array}\right.---\left(6\right)$

其中，μ_s为静态泊松比，E_s为静态弹性模量；

根据公式(6)将得到的所述动态三维岩石力学参数值转化为对应的静态值。

10.一种三维岩石力学参数获取系统，其特征在于，包括：

第一确定模块，用于根据目标工区已知井的声波测井资料获得横波时差和纵波时差，并确定所述横波时差和所述纵波时差间的横纵波时差关系；

第二确定模块，用于根据所述横纵波时差关系，确定与所述横波时差对应的横波波速，并确定与所述纵波时差对应的纵波波速；

第三确定模块，用于根据所述已知井的密度测井数据和所述纵波波速，确定岩石密度和所述纵波波速间的密度波速关系，并根据所述密度波速关系确定与所述纵波波速分别对应的岩石密度；

校正模块，用于根据所述纵波波速对所述目标工区的地震速度场数据进行校正；

计算模块，用于根据校正后的地震速度场数据和所述岩石密度计算三维岩石力学参数，得到动态三维岩石力学参数值；

转化模块，用于建立所述动态三维岩石力学参数值与岩芯试验得到的静态三维岩石力学参数值间的对应关系，并将得到的所述动态三维岩石力学参数值转化为对应的静态值。

11.根据权利要求10所述的系统，其特征在于，还包括：

展示模块，用于三维展布所述静态三维岩石力学参数值，并根据所述目标工区的地质特征信息确定所述静态三维岩石力学参数值是否合理；

调整模块，用于若所述静态三维岩石力学参数值不合理，则调整校正参数，重新根据所述纵波波速对所述目标工区的地震速度场数据进行校正。

12.根据权利要求10所述的系统，其特征在于，所述第一确定模块具体用于：

若所述声波测井资料中有偶极声波测井资料，则直接从所述偶极声波测井资料中获得所述横波时差；

若所述声波测井资料中有声波全波列测井资料，则从所述声波全波列测井资料中分离出所述横波时差；

从所述声波测井资料中获得所述纵波时差。

13.根据权利要求12所述的系统，其特征在于，所述第一确定模块还用于：

对所述横波时差和所述纵波时差进行统计分析，得到所述横纵波时差关系为：

Δt_s＝1.777Δt_p+0.546(1)

或者，若所述目标工区为中砂岩储层时，根据公式(2)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为：

${\mathrm{\Δt}}_s=\frac{{\mathrm{\Δt}}_p}{{[1-1.15\frac{(1/{\mathrm{\ρ}}_b)+{(1+{\mathrm{\ρ}}_b)}^3}{e^{1/{\mathrm{\ρ}}_b}}]}^{1.5}}---\left(2\right)$

其中，Δt_s为横波时差，Δt_p为纵波时差，ρ_b为中砂岩岩石密度；

或者，若所述目标工区为泥岩储层时，根据公式(3)确定根据所述横波时差和所述纵波时差间的横纵波时差关系为：

Δt_s/Δt_p＝B-0.8(ρ_sh-ρ_shmin)/(ρ_shmax-ρ_shmin)(3)

其中，B为经验常数，ρ_sh为泥岩岩石密度，当ρ_sh≤2.2g/cm³时，B＝2.5；

当ρ_sh≥2.65g/cm³时，B＝1.7。ρ_shmin常取2.2g/cm³，ρ_shmax常取2.65g/cm³。

14.根据权利要求13所述的系统，其特征在于，所述第二确定模块具体用于：

根据公式(4)确定与所述横波时差对应的横波波速，并确定与所述纵波时差对应的纵波波速：

V_s＝0.3046×10⁶/Δt_s

V_p＝0.3046×10⁶/Δt_p(4)

其中，V_s和V_p分别为横波波速和纵波波速。

15.根据权利要求10至14中任一项所述的系统，其特征在于，所述第三确定模块具体用于：

采用指数模型对所述密度测井数据和所述纵波波速进行统计分析，确定岩石密度与所述纵波波速间的密度波速关系。

16.根据权利要求15所述的系统，其特征在于，所述校正模块具体用于：

根据所述纵波波速，采用距离加权反比插值法对所述目标工区的地震速度场数据进行校正。

17.根据权利要求16所述的系统，其特征在于，所述岩石力学参数包括泊松比和弹性模量，所述计算模块具体用于：

根据公式(5)计算得到动态泊松比和动态弹性模量：

$\left\{\begin{array}{c}{\mathrm{\μ}}_d=\frac{V_p^{\′2}-{2V}_s^{\′2}}{2(V_p^{\′2}-V_s^{\′2})}\\ E_d={10}^{-3}{\mathrm{\ρV}}_s^{\′2}\frac{{3V}_p^{\′2}-{4V}_s^{\′2}}{V_p^{\′2}-V_s^{\′2}}\end{array}\right.---\left(5\right)$

其中，μ_d为动态泊松比，E_d为动态弹性模量，ρ为岩石密度，V'_p和V'_s分别为校正后的纵波波速和横波波速。

18.根据权利要求17所述的系统，其特征在于，所述转化模块具体用于：

对所述动态三维岩石力学参数值和所述静态三维岩石力学参数值进行统计回归分析，得到回归公式(6)：

$\left\{\begin{array}{c}{\mathrm{\μ}}_s={\mathrm{\μ}}_d\\ E_s=0.7639E_d+4447.2\end{array}\right.---\left(6\right)$

其中，μ_s为静态泊松比，E_s为静态弹性模量；

根据公式(6)将得到的所述动态三维岩石力学参数值转化为对应的静态值。