速度反射率和密度反射率的分解方法
【技术领域】
[0001] 本发明涉及油气勘探技术领域，尤其涉及一种速度反射率和密度反射率的分解方 法。
【背景技术】
[0002] 众所周知，探油目标是低密度、高孔的岩石，但在许多情形下储油岩的速度是高 的，特别是横波速度，从而使其阻抗（密度与速度的乘积）与围岩阻抗的差别较小，降低反 射振幅强度，使之难于探测，解释结果极不确定，导致可能误判油层或者漏判油层。问题的 症结来源于速度反射率（包括纵波速度反射率与横波速度反射率）与密度反射率计算的矩 阵（Aki&Richard'sequation)是病态的或者接近线性相关的，只能通过重组方程和矩阵， 近似解出声波阻抗和横波阻抗反射率二个独立变量，而无法解出三个独立变量。
[0003] 现有的解决方法通常有两种：一种是加大测量的反射角度，改善方程和矩阵使之 由病态变为良性。第二种方法是添加较为严格的约束条件，使矩阵由病态变为良性。
[0004] 第一种方法的缺陷在于，除非角度很大（60度），否则矩阵改善程度有限，即秩仍 然趋向于2而不是3,而且高角度资料不常见，即使有，信噪比太低极不可靠。
[0005] 第二种方法的缺陷在于，较为严格的约束条件的变化很大，取决于时间及位置，添 加不合适的约束条件会解出不可靠甚至是完全错误的解。
【发明内容】
[0006] 在下文中给出关于本发明的简要概述，以便提供关于本发明的某些方面的基本理 解。应当理解，这个概述并不是关于本发明的穷举性概述。它并不是意图确定本发明的关 键或重要部分，也不是意图限定本发明的范围。其目的仅仅是以简化的形式给出某些概念， 以此作为稍后论述的更详细描述的前序。
[0007] 本发明提供一种适用性较强、抗噪音、可靠性较高的速度反射率和密度反射率的 分解方法。
[0008] 本发明提供一种速度反射率和密度反射率的分解方法，包括：
[0009] 在不同探测入射角度测量获取相应的探测参数，根据所述探测参数和第一公式建 立第一矩阵方程；
[0010] 对所述第一矩阵方程添加宽松约束条件，得到第二矩阵方程；
[0011] 对所述第二矩阵方程进行计算求解，得到横波速度反射率视解，并根据横波速度 反射率视解与横波速度反射率真解之间的近似线性关系计算出横波速度反射率近似解；
[0012] 重组所述第一公式得到第二公式，根据所述探测参数和第二公式得到第三矩阵方 程；对所述第三矩阵方程添加所述宽松约束条件，得到第四矩阵方程；
[0013] 对所述第四矩阵方程进行计算求解，得到纵波阻抗反射率近似解和横波阻抗反射 率近似解；
[0014] 根据所述横波速度反射率近似解、所述纵波阻抗反射率近似解和所述横波阻抗反 射率近似解计算得到密度反射率近似解和纵波速度反射率近似解。
[0015] 本发明提供的速度反射率和密度反射率的分解方法通过对测量得到的矩阵方程 添加宽松约束条件，得到与横波速度反射率真解具有近似线性关系的横波速度反射率视 解，从而计算出较为精确的横波速度反射率近似解，再通过重组矩阵方程计算出纵波阻抗 反射率近似解和横波阻抗反射率近似解，最终计算出较为精确的密度反射率近似解和纵波 速度反射率近似解。本发明所添加的宽松约束条件只需要合理的纵波速度反射率与横波速 度反射率的线性关系（可从区域资料或文献中获得），具有普遍的适应性。通过与真实的速 度反射率和密度反射率对比验证，可证实本发明方法所获得的解可靠性高。综上所述，本发 明具有适用性较强、抗噪音、可靠性较高等优点。
【附图说明】
[0016] 参照下面结合附图对本发明实施例的说明，会更加容易地理解本发明的以上和其 它目的、特点和优点。附图中的部件只是为了示出本发明的原理。在附图中，相同的或类似 的技术特征或部件将采用相同或类似的附图标记来表示。
[0017] 图1为本发明速度反射率和密度反射率的分解方法的流程图。
[0018] 图2为本发明速度反射率和密度反射率的分解方法中步骤S10的流程图。
[0019]图3为不添加约束条件直接对第一矩阵方程运算得到的横波速度反射率视解与 横波速度反射率真解的关系示意图。
[0020]图4为添加严格约束条件后运算得到的横波速度反射率视解与横波速度反射率 真解的关系不意图。
[0021]图5为本发明添加宽松约束条件后对第二矩阵方程运算得到的横波速度反射率 视解与横波速度反射率真解的关系示意图。
[0022] 图6为本发明运算得到的密度反射率视解与密度反射率真解的关系示意图。
[0023]图7为本发明运算得到的纵波速度反射率视解与纵波速度反射率真解的关系示 意图。
【具体实施方式】
[0024] 下面参照附图来说明本发明的实施例。在本发明的一个附图或一种实施方式中描 述的元素和特征可以与一个或更多个其它附图或实施方式中示出的元素和特征相结合。应 当注意，为了清楚的目的，附图和说明中省略了与本发明无关的、本领域普通技术人员已知 的部件和处理的表示和描述。
[0025] 图1为本发明速度反射率和密度反射率的分解方法的流程图。
[0026] 如图1所示，在本实施方式中，本发明速度反射率和密度反射率的分解方法包括：
[0027] S10:在不同探测入射角度测量获取相应的探测参数，根据所述探测参数和第一公 式建立第一矩阵方程；
[0028] S20 :对所述第一矩阵方程添加宽松约束条件，得到第二矩阵方程；
[0029] S30 :对所述第二矩阵方程进行计算求解，得到横波速度反射率视解，并根据横波 速度反射率视解与横波速度反射率真解之间的近似线性关系计算出横波速度反射率近似 解；
[0030] S40:重组所述第一公式得到第二公式，根据所述探测参数和第二公式得到第三矩 阵方程；对所述第三矩阵方程添加所述宽松约束条件，得到第四矩阵方程；
[0031] S50:对所述第四矩阵方程进行计算求解，得到纵波阻抗反射率近似解和横波阻抗 反射率近似解；
[0032] S60:根据所述横波速度反射率近似解、所述纵波阻抗反射率近似解和所述横波阻 抗反射率近似解计算得到密度反射率近似解和纵波速度反射率近似解。
[0033] 图2为本发明速度反射率和密度反射率的分解方法中步骤S10的流程图。
[0034] 如图2所示