地震雷达定量整合确定孔隙度和含水饱和度.docxVIP

定量整合地震和雷达反射以得到对土壤孔隙度和含水饱和度的唯一估计整合定量的地震和探地雷达角度依赖反射系数。已经通过数值模拟验证了这种方法，结果显示，此方法能够得到底土边界附近的准确的孔隙度和含水饱和度，而单独的地震或是雷达方法都不能得到准确的结果。此方法还可应用于水文地质环境、农业、岩土工程等领域。1、简介随着硬件以及各种技术的提高，地震分辨率很高，在好的条件下甚至能达到分米级，地震波长很接近雷达波长。此外，浅层地震信号的振幅保真提高了，这对浅层地震资料的AVO分析很有利。雷达数据的AVO分析近年来有很多人在做。用多次偏移的雷达数据进行高频电磁波在浅层的区间速度估计。利用这个速度，在一些假设的基础上，能够估计饱和土层中水的体积，水的体积是由孔隙度Φ以及含水饱和度Sw得出的。然而，单独从雷达数据中是不能分别得到孔隙度Φ以及含水饱和度Sw的。地震和雷达速度之间的经验关系近来被应用于求得岩石物理性质，浅层地震和雷达反射的相关性，被用于尝试提高解释精度，但这些解释和方法的综合都是定性的。在此，我们提出定量整合角度依赖的地震和雷达反射数据的方法，可以唯一估算浅层底土的孔隙度Φ以及含水饱和度Sw。估计浅层底土的孔隙度Φ以及含水饱和度Sw以及他们的横向变化在水文地质和环境岩土工程的应用中至关重要。缺少这两个参数的准确估计会在地质工程的施工中冒很大风险。提出一个共享地球模型，推导出弹性波和电磁波在孔隙度Φ和含水饱和度Sw方面的体积属性，以及已知的固、液、气相成分的属性值。提出了新方法定量合成地震和雷达数据，此方法在合成数据中进行测试。2、地震-地质雷达共享地球物理模型层界面的地震反射系数从物质的拉梅系数（λ、μ）和体积密度（ρ）得到。雷达波反射系数主要通过物质相对介电常数（ε）得到，假定电导率为常数，磁导率为自由空间磁导率。利用布莱格曼混合法则估算三相系统（气、液、固）的拉梅系数。混合相中的拉梅系数表示如下：在对几何信息的缺失的成分的分布， Hashin-Shtrikman（HS）边界，提供最好的（最窄范围）介质上部和下部的预测范围有效模量。从混合规则得出的有效体积特性估计与HS界限相符。在Φ值低，估计值接近HS上部边界值（由固相决定），而Φ增加时，估计值接近HS下界（由液相决定）。对于含水饱和度Sw，变化由流体成分决定。本研究采用混合模型，因此给出了合理的结果。图1显示了孔隙度Φ、含水饱和度Sw与地震和电磁波速度的函数关系。图中显示速度变化是现实的，表明模型很可靠。值得注意的是：孔隙度Φ、含水饱和度Sw是自由变量，用这两个变量完成弹性和电磁属性的合成。3、合成方法雷达和地震数据采集都是在相同的环境和地方进行，在一个给定的层位或几何位置进行数据采集。两个数据集中，波长相似。采用CMP测量，数据已经经过预处理（去噪等），然后将地震数据分解成SH、SV、P模式，将雷达数据分解成TE、TM模式。地震和雷达反射的角度依赖反射系数都与孔隙度Φ、含水饱和度Sw成函数关系。我们选择SH+TE 和SV + TM模式组合，因为在同一测线的情况下，雷达和地震的极化方向一致，能够对地下介质进行相似的探测；SV + TM模式也是如此。假定不考虑损失的影响，地震反射系数与频率无关，角度依赖反射系数的理论表达式都给出了，例如，Aki and Richards[1980] for seismic waves and by Kong [1986] for electromagneticwaves.最后，合成是通过在损耗方程中整合地震和雷达反射系数完成的：用全局非线性左布里兹和菲涅尔反射系数反演孔隙度Φ和含水饱和度Sw。4、数值模拟：结果利用上述公式计算了多个炮检距情况下的地震SV、SH、P-SV反射系数以及TETM模式的电磁反射系数。合成的数据不受界面影响。CMP道集转换成入射角道集，0到45度的范围内有13道。模型分为两层，第一层的参数已知，第二层参数未知，利用给定的孔隙度和含水饱和度进行估算。从图2中可以看出，单独的SH或TE数据能够得到一系列的孔隙度和含水饱和度组合使得方程（3）最优化，得到最小解，这样的不唯一性太严重，并且不能看出图中哪里是最小值。所以单独使用地震SH或者雷达TE数据是不能得到孔隙度和含水饱和度的唯一解的。当利用地震SV数据和雷达TM数据时得到的结果类似。当使用合成数据进行计算时，可以从图3中看出，孔隙度和含水饱和度的组合取得最小值就比较唯一了，两种模式组合取得的结果很相似。本研究中只考虑了孔隙度Φ和含水饱和度Sw两个变量，增加变量此重建方法同样能够适用。将此方法用于实际数据的挑战有：采集完全相同条件下的地震和雷达数据，以及正确的处理流程获得可靠的反射振幅。5、结论