【地球物理中的边界单元法】徐世浙科学.pdf

地球物理中的边界单元法徐世折著

序有限单元法是一种有效的数值方法,但它是区域型的方法,难于处理地球物理中经常遇到的三维问题和无限区域问题,边界单元法只在研究区域的边界上部分单元,使求解问题的维数降低,最后形成的代数方程组的元数大为减少．尤其对于三维问题和无限区域问题,这个优点很突出．所以,边界单元法已成为地球物理数值计算的主要方法之一,地球物理的数值计算中,先后发表了有关论文约30篇,解决了许多地球物理计算中的难题:例如,电阻率法中的三维地形改正问题,国外用有限单元法求解,需解含有50000个未知数的线性代数方程组.

引言地球物理勘探的自的是了解地下地质体的分布,即把地球物理场的观测资料转变为地质资料.其第一步1.作是地球物理反演,即从地球物理场值反演出地下的物性分布,然后,结合地质情况,将物理特性“翻译”成地质特性,从而获得地下地质体分布的资料.可见,反演在地球物理中占有重要的位置.地球物理反演必须以正演为前提所谓正演,就是给定场源和物性的分布,求解场的分布.反演时,将计算值与观测值进行对比,如果两者相差太大,就修改原来的模型,直到计算的场值与实测的场值最佳拟合为止.由此选择出的模型被看作是最佳的反演结果.