地球物理正演中的无网格法算法研究以二维大地电磁法为例

来源 :长江大学 | 被引量 : 0次 | 上传用户:bell900818
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
大地电磁测深法(Magnetotellurics,MT)以岩石的电性差异为基础和前提,利用天然交变电磁场研究地球电性结构的勘探方法.由于它不用人工供电,成本低,工作方便,不受高阻层屏蔽,而且勘探深度随着电磁场的频率而异,浅可以几十米,深可达数百公里,因此受到广泛的应用.
  在解决人地电磁测深法(MT)正演问题时,有限差分法和有限元法为主要的数值计算方法,但是这些问题都是基于网格的数值方法,如有限元法是将求解域剖分成有限个单元,各个单元通过网格节点连接在一起,基于网格剖分是已有的传统的数值计算方法的基本特征.尽管这些方法取得了巨大成功,但是也有其白身的局限性:1)形成FEM网格时的计算成本高,划分单元是一项十分繁琐的工作,对于复杂的三维问题,其工作量往往比有限元分析本身还大.因而不适应需要不断变动网格的情况.2)有限元(FEM)的计算精度依赖于单元剖分的形状和大小,因而不适应因大变形而发生网格畸变的情况.
  产生上述问题的一个主要原因是这些方法都是基于网格剖分的基础上来建立离散系统方程的,为克服这些困难,无网格方法(Meshfree Method)作为有限元(FEM)的一种重要补充具有重要的研究意义.
  无网格法(Meshfree Method)是近十年米兴起的一种数值计算方法[7],它是在有限元的基础上发展起来的,但是不同于有限元,无网格法的近似函数是建立在一系列离散节点上的.无网格法已在在力学、油藏渗流问题、电磁学等领域得到了广泛的研究,然而在求解地球物理正演问题中的理论和引用研究都还未出现,本文利用无网格法研究二维情况下的MT正演,用Matlab为程序编译工具.本文主要从以下几个方而做了相关研究
  首先从麦克斯韦方程组出发,通过传统的变分原理推导出二维介质大地电磁(MT)边值问题所对应的积分方程,然后用一系列的离散节点离散求解区域,得到建立在离散节点上的无网格形函数,引入了移动最小二乘方法来构造形函数推导出了形函数的求取方法并结合二维大地电磁勘探,给出了相关参数,如影响域、定义域、权函数等的确定方法.
  由于采用的是最小二乘的思想,因而在任意的一点是不满足德尔塔函数的性质,因而说要利用特殊的方法使得在这些点处满足德尔塔,本文采用了罚函数法,拉格朗日乘子法来处理本质边界条件.
  将上述方法中求得的形函数与代入变分问题得到边值问题对应的离散方程,解方程后代入视电阻率计算公式,得到地面上各点的视电阻率.
  最后通过几个典型模型的卡利亚视电阻率拟断面图和相位拟断面图,得出不同模型下电磁场相应特征的分布规律,验证了无网格算法的正确性.
其他文献
随着经济的日益发展,海水利用行业、食品加工业、化学制造业等化学加工业于生产过程中所排放的高盐废水量与日俱增。国内外学者针对高盐废水的处理已开展相关研究,但普遍存在系统稳定性不足、污染物去除效率较低、菌群活性较差等问题。依据含盐废水的实际特性,课题研究不同盐度条件下SBR好氧颗粒污泥反应器的启动过程,优化工艺运行参数条件,考察颗粒污泥系统的抗冲击能力,并利用高通量测序手段分析系统中微生物的多样性、差
学位
昆虫围食膜是中肠上皮细胞分泌的由几丁质和蛋白质构成的薄膜,具有两个重要的生理功能。第一,围食膜是天然的物理屏障,其隔离食物与上皮细胞,可抵御外源有害物质的入侵。干扰或破坏围食膜结构,将促进病原微生物或杀虫剂对昆虫的侵害,影响昆虫正常发育,导致死亡。因此,围食膜是害虫防治中重要的作用对象,参与围食膜形成和维持其结构稳定的关键蛋白被认为是潜在的生物防治靶标。第二,围食膜为中肠不同区室的功能分化提供结构
该论文致力于正丁烷氧化制顺丁烯二酸酐流化床催化剂的制备过程研究,考察了不同工艺条件对催化剂性能的影响.在该实验组前期工作的基础上,以正丁醇同时作溶剂和还原剂,采用统计实验法确定了前驱体制备过程中滴加酸浓度、反应温度、反应时间、溶剂/物料比(丁醇cm/VOg)等因素对催化剂性能的影响.并研究了在较高温度下,不同酸浓度对催化剂性能的影响.确定了最佳的酸浓度范围和最佳反应温度.为在工业上的推广作了一定的
学位
该文采用Cl为氯化剂,在可见光(白炽灯)催化下,气相法合成ClCHSiCl.研究了光源特性(如光源强度,波长范围),反应条件(如反应温度、反应时间,物料分离参数,反应器材质及壁厚,Cl通量),及合成工艺等对光氯化反应的影响,从而确定了氯化反应所需要的能量及反应参数,并且从理论上对反应机理,特性进行了论述.对反应机理的研究,可以使我们明确研究方向,找出影响反应的主要因素.以合成方法的研究,可以提高反
该文报导由银配合物制备AgI均匀球形纳米粒子.AgI纳米粒子粒径均小于100nm.研究反应条件(银离子浓度、配体氨浓度、反应时间等)对AgI粒子形态的影响.并将AgI纳米粒子作为内置核,将Y(OH)CO进行升温水解,形成Y(OH)CO复合粒子,用高浓度NaSO溶解内置核,最终得到Y(OH)CO中空粒子.
学位
该论文利用无定形的混合稀土水合氧化的与丙烯酰胺采用反相悬浮聚合的方法,制备出了REO·2.48HO-PAM球状除氟剂.探讨了影响制备REO·2.48HO-PAM球状除氟剂的因素,如油相种类和油相的粘度,水油比,交联剂的浓度,引发剂的种类及用量,聚合反应温度及体系升温速率大小对造粒的影响和表面活性剂的种类、用量对REO·2.48HO-PAM球状除氟剂粒径,包覆率,颗粒形状等方面进行了细致的研究.最后
该文对超临界CO萃取乙醇水溶液过程进行了一系列的研究.首先针对该方法我们设计了三种实验方法:简单萃取,萃取精馏和近临界萃取.我们将三种方法的实验结果和普通蒸馏方法比较,发现三种方法均优于普通蒸馏,其中近临界萃取和萃取精馏方法优于简单萃取,但是从能耗的角度考虑近临界CO萃取是最具有优势适合工业化的方法.同时利用已有的两相相平衡数据和UNIQUAC方法我们计算了乙醇-水-CO体系的相平衡,并且和实验数
学位
江汉盆地位于湖北省江汉平原中部,面积36360平方公里,是我国东部一个典型的内陆盐湖断陷盆地,可划分为11个凹陷和4个凸起共计15个二级构造单元。江陵凹陷位于江汉盆地西南部,是江汉盆地最大的二级构造单元。以万城断层为界,具有NW分带、NE分块的构造格局,发育众多断陷、洼地,整体呈现多隆多坳的构造格局。江陵凹陷虽经过多年的勘探与研究取得一定的突破,长期以来认为凹陷内的梅槐桥洼陷生烃最好,勘探程度极不
目前,世界范围内的各大油气田,特别是中国东部油田新增油气地质储量多数来源于隐蔽油气藏。其勘探开发难度越来越大、成本越来越高,主要原因在于目前缺乏非常有效的此类隐蔽性储层的预测及描述方法与技术。因此,迫切需要研究和开发适应于当前勘探开发阶段的高精度储层预测方法。  地震岩性参数的准确反演一直是地球物理工作者长期的奋斗目标,在油田开发地震工作中,无论是储层预测还是油气圈定、油藏描述,波阻抗反演都是极其
学位
随着我国各地含有油气藏勘探程度的不断深入,陆地和海上油气藏的勘探幵发研究的不断发展,特别是对岩性油气藏的幵发的概率越来越大,对于特性岩性的研究更是益渐盛。根据目前资料掌握的情况和勘探要求,本次研究地区为位于黄骅坳陷北大港构造带北翼及向北倾没延伸的外围地区,主要目的层为沙河街组多套发育的灰岩层段,并将砂岩类储层列为辅助目的层。  当前,针对岩性地层的油气藏勘探需使用有效的技术与方法的支持。近些年,国