用有限元法计算二维重力场垂直分量及重力位二阶导数

来源 :石油地球物理勘探 | 被引量 : 0次 | 上传用户:xbalr
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
有限元法可以计算密度分布、形态复杂物体的重力场垂直分量 g 及重力位二阶导数 W_x■、W_(yy)。取一个包围密度体的足够大的区域,求解 g 的边值问题可表为▽~2g=-4πK(■ρ)/(■y) 在区域内 (■)g)/(■n)+(sin(θ-α))/(r sin θ)g=0 在边界上与上述边值问题相应的变分问题是泛函F(g)=■[(1)/(2)(▽g)~2+4πK(■g)/(■y)ρ]dS+■_Γ(1)/(2)(sin(θ-α))/(r sin θ)g~2dl取极值。用有限元解上述变分问题时,将区域Ω剖分为三角单元,在单元 e 内进行二次函数插值。首先计算各单元的 F_e(g),然后相加组成总体的 F(g),它是各节点待求的 g 的函数。对 F(g)求极值,得一线性代数方程组。解方程组可得各节点的 g。对 g 进行微商,即可得重力位二阶导数。 The finite element method can calculate the density distribution, the vertical component g of gravitational field and the second derivative W_x ■, W_ (yy) of the gravitational field with complex shape. Taking a sufficiently large area around the density body, the problem of solving the boundary value of g can be expressed as ∇ ~ 2g = -4πK (■ ρ) / (■ y) in the area g) / (■ n) + The variational problem corresponding to the above boundary value problem at the boundary is F (g) = ■ [(1) / (2) (∇g) ~ 2 + 4πK (■ g) / (■ y) ρ] dS + ■ _Γ (1) / (2) (sin (θ-α)) / (r sin θ) g ~ 2dl. When finite element method is used to solve the above variational problem, the area Ω is divided into triangular units and quadratic function interpolation is performed in unit e. First calculate F_e (g) for each cell and then add F (g) that makes up the population as a function of g to be sought by each node. Find the extreme value of F (g) to obtain a linear algebraic system. Solution of equations available for each node g. Derivative of g, you can get the second derivative of gravity.
其他文献
CO_2气体保护焊在焊接过程中,由于电压、电流、焊丝伸出长度和电弧稳定性等因素的影响,一般都要产生飞溅。这些飞溅的渣滓牢固地粘附在喷嘴内壁,堵塞了喷嘴的气流孔道,使保
(89)法国《解放报》(Libération)发表关于第二次世界大战的专题采访,其中讲到诺曼底地区居民回忆当年美国军队在诺曼底登陆的情景。一个退休的技术工人说:Quandjaivuquelaroute,desdeuxctés,nétaitplus... (89) Libération, the
以硫脲为络合剂,在低酸性介质中原子吸收测银已被用于矿石样品分析中、文献曾对矿石中一定量的Fe~(3+)、Cu~(2+)、Mn~(2+)、Au~(3+)、Ni~(2+)、Al~(3+)、Ca~(2+)、Cr~(6+)、
根据大学英语教学改革的需要和众多教师的要求,上海外语教育出版社(以下简称“外教社”)决定,把以往十多年来年年分散在各省市举行的小范围教学研讨会改为集中的、较大规模的教师
一、前言近年来国内外在深孔加工技术上有较大发展。深孔加工(L/D≥10)致力于缩短加工时间,发展了单刃、多刃扩钻、单刃空心钻、深孔钻等,这些钻头成功地解决了深孔加工的排
本文报道了自行建立的用热敏电阻作为测温元件的精密溶解量热装置,并用量热标准物质光谱纯氯化钾和美国国家标准局提供的三羟甲基氨基甲烷(TRIS)对仪器进行了标定,在25℃所得
随着水文地质计算工作的进展和生产实践的需要,参数导水系数T越来越被重视。尤其是在无观测孔的单孔抽水情况下,如何计算T值成了一个重要问题。本文试图通过一些试验资料来
一般来说,贵金属试样的分解比其它常见金属困难。在应用的几种分解方法中,酸溶法一般只能分解银、金、钯、铂及其矿物;所有贵金属试样虽然都可用碱、硫酸氢盐或过氧化物等熔
一.前言利用重力资料,可以获得大面积地壳深部构造轮廓。本文根据重力场的向上解析延拓方法,分析讨论了安宁河地区(北纬26°21′~27°49′、东经100°18′—105′)向上延拓(h
测定锌、铝中痕量杂质是鉴定纯度的基础工作。玻璃碳汞膜电极微分脉冲阳极溶出伏安法(DPASV法)已被用于测定水中铜、铅、镉。但是,纯锌、铝中砷、锑、铋、锡等对测定痕量的