预修正快速傅里叶变换高阶边界元方法在波浪对物体作用问题中的应用

边界元方法被广泛应用于波浪对海上婕筑物作用领域,但由于传统边界元方法的存储量和计算量均为未知量的平方量级,很难满足大范围多未知量问题的计算需要.本文基于高阶边界元方法,应用预修正快速傅里叶变换方法,使计算量与存储量分别降低至O(NlogN)和O(N)量级,并得到一个连续的压强分布以适应结构设计的要求,同时可以通过使用满...     

基于波动方程有限差分算法的接收函数正演与偏移

针对接收函数正演与偏移, 本文采用波动方程有限差分算法. 借鉴成熟的勘探地震学方法, 引入等效速度概念, 建立接收函数转换波与地震勘探反射波的等效走时方程, 实现了基于波动方程有限差分算法的接收函数正演与偏移. 数值计算表明, 波动方程有限差分叠后偏移方法可以对点绕射和穹隆构造模型实现高精度成像. 本文利用数值计算讨论了波动方程有限差分叠后偏移与Kirchhoff叠后偏移对于接收函数偏移的适用性, 还对偏移过程中速度模型的误差进行了分析.     

基于高密度电阻率法对季节性冻土的监测研究

通过建立正演冻土模型,验证了将高密度电阻率法应用于冻土监测的可行性。对高密度电阻率法装置进行特殊处理,将其应用于对季节性冻土的监测,并进行了七个月份的连续监测。结果表明:高密度电阻率法在季节性冻土监测中可行;通过对实测数据进行反演处理,得到电阻率反演剖面,可以看出表层土壤在季节变化中有明显的结冻、解冻行为;并可估算出冬季冻土深度范围。     

中条山地区涑水杂岩新太古代烟庄正长花岗岩年龄及成因:对华北克拉通地壳演化的制约

太古宙末期钾质花岗岩的广泛发育是陆壳成熟和稳定化的重要标志,对了解早期陆壳的形成与演化具有重要的意义.发育于华北克拉通南缘中条山地区涑水杂岩中的烟庄正长花岗岩的形成年龄和成因还没有被很好地限定,构造背景还存在争议.对烟庄花岗岩进行了锆石U-Pb年代学和Hf同位素以及全岩地球化学和Nd同位素研究.烟庄花岗岩锆石SHRIMP U-Pb年龄为2 515±7 Ma.岩石具有高硅(SiO₂=73.05%~74.85%)、高钾(K₂O=4.46%~5.86%)、富碱(ALK=8.32%~9.36%)、贫钙(CaO=0.55%~0.98%)、低TFeO*(0.73%~1.28%)和MgO(0.31%~0.52%)的特征,A/CNK=1.01~1.04,为弱过铝质的钾玄系列.稀土总量变化较大(ΣREE=63.80×10-6~250.02×10-6)),轻重稀土元素分异明显((La/Yb)_N=25.44~92.87),Eu异常变化较大(Eu/Eu*=0.47~0.79).岩石低Sr、Ba,富集Rb、Th、U等元素,亏损Nb、Zr、Y、Yb、Cr、Co、Ni等元素,具有较高的Rb/Sr、Rb/Ba和Sr/Yb比值以及较低Sm/Nd和Nd/Th比值,具有高分异I型花岗岩的特征.烟庄花岗岩具有0附近的全岩ε_Nd(t)值,岩浆锆石具有正的ε_Hf(t)值(2.85~3.66),两阶段Hf模式年龄为2 258~2 883 Ma,多数在2 600~2 883 Ma之间.结合其他方面研究,烟庄花岗岩具有同碰撞和后碰撞花岗岩的特征,推测为新生地壳在由挤压向伸展转换的构造背景下部分熔融所形成,可能有少量地幔物质添加.这期钾质花岗岩的形成,标志着华北克拉通太古宙末期强烈岩浆活动的结束以及稳定陆壳的形成.      

一种层析成像的正演和反演方法

作者在本文中提出了一种层析成像的正演和反演方法，首先对求解区作网络划分，由给定的正演速度模型在网格点上双线性插值，获得整个网格节点的速度分布，然后，对波阵面作一近似，在此基础上，用最小走时最短路径射线追踪法完成正演，反演则是通过正演的方法实现的，理论模型算例表明，本方法效果较好。     

夹层厚度对瑞利面波频散曲线的影响

文章以弹性波横波速度为参数,利用瑞利面波频散曲线的快速标量传递算法进行正演模拟计算,系统研究软弱层厚度对瑞利面波基阶模频散曲线的影响,分析和总结了频率与相速度（f-c）、频率与波数（f-k）、深度与相速度（h-c）等瑞利面波基阶模曲线随软弱层厚度变化的规律。     

全局弱式无网格法在大地电磁二维正演中的应用研究

无网格法作为有限元法的一种重要补充和发展具有重要的研究意义和实用价值,其避免了网格剖分,具有前期处理简单,精度高,自适应分析方便,适用于复杂模型和复杂边界条件等优点.本文重点研究基于全局弱式无网格法的大地电磁二维正演问题.本文首先论述了无网格理论,并推导出基于全局弱式无网格法的大地电磁二维问题对应的最终总体矩阵表达式,通过标准模型数值模拟实验研究了相关重要参数对计算精度及计算效率的影响,并确定其最优值,选取高斯点数量为2～3、无量纲支持域尺寸为1.1～1.3.后续分别用有限元法和全局弱式无网格法计算典型一维模型和二维模型,通过对比计算结果及计算精度,验证本文中提出方法在处理大地电磁二维问题上的有效性及在处理复杂模型上的优越性.     

基本结构有限元算法及大地电磁测深一维连续介质正演

提出了一种新的数值方法——基本结构有限元法. 从经典的伽辽金法（Galerkin method）出发,首先定义了基本结构插值基函数,在此基础上建立了基本结构方程,然后用有限元法进行进一步的具体实现. 该算法被成功应用到大地电磁测深一维连续介质正演计算中. 采用了6种不同的单元插值基函数进行计算,并对计算结果进行了比较和分析.     

我国及邻近地区的嵌套式磁场模型

本文将球谐分析（SHA）和矩谐分析（RHA）依次用于空间范围逐步缩小的区域,得到了嵌套式磁场模型。这个模型可以清楚地反映出大小不同的各级地磁异常的三维结构,为研究地磁异常起源及有关的大地构造问题提供了较好的依据。本文建立的1970年代华北地区地磁场的三级嵌套模型由全球SHA模型、中国RHA模型和华北RHA模型三级组成。该模型指出,华北中部由南向北有一个由北北东逐渐转向北东东方向的Z分量正异常带,这可能与华北的地质构造有关。     

全球月平均高度距平场正斜压性的分析(I)--相关系数分布

利用欧洲中心全球客观分析的7层月平均位势高度距平资料,分别计算了相同格点上12-2月和6-8月各层高度距平之间的同时相关分布。结果表明:全球月平均位势高度距平场的正压、斜压性随地理区域变化很大,并有较明显的季节变化。其中12-2月对流层下部大气环流距平场在海洋上空有较好的相当正压性,在大陆上空正压性较差;对流层上部大气环流异常的正、斜压性分布与下部不同,赤道附近地区是相关系数的高值区,副热带地区主要是低值区;在对流层下部与上部大气环流距平场的相关系数分布图上,海洋上空并不总是高值区,在热带东太平洋区域就出现最强的负相关中心。与6-8月的相关系数分布型比较发现:特别明显的季节变化在非洲大陆和热带印度洋区域,12-2月为强正压区,6-8月为强斜压区。