有限元法逼近局域大地水准面的数学模型

顾及到大地水准面的物理特性,依据其在局部所应满足的数学物理方程,讨论了有限元法逼近局域大地水准面的数学思想,给出了相应的数学模型.     

点力源在横向各向同性介质中激发的弹性波

点力源在均匀横向各向同性介质中激发的弹性波可以用透射-反射矩阵法进行求解。假定点力源位于一假想的各向同性层中,并利用各向同性介质中等效于势函数的过源面间断,从而易于求出介质中各个点上的弹性波场;再令假想的各向同性层的厚度趋近于0,则可求出原始问题的解,即横向各向同性介质中的格林函数。最后得出的结果和假想的各向同性层的性质无关,表明这一方法用于求解点力源在各向异性介质中激发的弹性波场是可行的。     

基于梯度塑性模型的多孔介质应变局部化分析的基本理论

讨论了引入梯度塑性模型情况下的含液多孔介质双重内尺度律的特征以及相互作用问题。在引入梯度塑性本构模型的情况下,一方面材料应变局部化分析数值结果的正则性得到保证;而另一方面,含液多孔介质会出现双重内尺度律参数的共同作用问题。本文给出了此时内尺度律预测的一个基本方法,并对稳定性问题进行了分析,讨论了不同情况下实波速存在的条件,给出了对于给定的渗透系数情况下实波速存在的波数区间,并对相关现象进行了解释。     

An Application of the Adjoint Method to a Statistical-Dynamical Tropical-Cyclone Prediction Model （SD-90） Ⅱ： Real Tropical Cyclone Cases

In the first paper in this series, a variational data assimilation of ideal tropical cyclone （TC） tracks was performed for the statistical-dynamical prediction model SD-90 by the adjoint method, and a prediction of TC tracks was made with good accuracy for tracks containing no sharp turns. In the present paper, the cases of real TC tracks are studied. Due to the complexity of TC motion, attention is paid to the diagnostic research of TC motion. First, five TC tracks are studied. Using the data of each entire TC track, by the adjoint method, five TC tracks are fitted well, and the forces acting on the TCs are retrieved. For a given TC, the distribution of the resultant of the retrieved force and Coriolis force well matches the corresponding TC track, i.e., when a TC turns, the resultant of the retrieved force and Coriolis force acts as a centripetal force, which means that the TC indeed moves like a particle; in particular, for TC 9911, the clockwise looping motion is also fitted well. And the distribution of the resultant appears to be periodic in some cases. Then, the present method is carried out for a portion of the track data for TC 9804, which indicates that when the amount of data for a TC track is sufficient, the algorithm is stable. And finally, the same algorithm is implemented for TCs with a double-eyewall structure, namely Bilis （2000） and Winnie （1997）, and the results prove the applicability of the algorithm to TCs with complicated mesoscale structures if the TC track data are obtained every three hours.     

利用卫星测高资料确定海洋重力场

卫星测高是20世纪末发展起来的新型星地观测技术。它是利用星载雷达测高仪测量卫星到海面的高度、有效渡高和后向散射系数，从而大规模探测海洋重力场、海洋环流、海洋风场、海底地质和地貌等海洋环境信息。西安测绘研究所在国内率先利用Topex／Poseidon、     

柯西波场及在角度域成像条件中的应用

角度域成像道集是叠前深度偏移的重要输出结果,它是偏移速度分析、各向异性分析和AVA分析等研究工作的基础.目前存在的角度域成像道集的生成方法受计算效率或角度分辨率的影响,仍然满足不了实际生产的要求.角度域成像道集的生成方法可以大致分为直接法和间接法两大类.在直接方法中,波矢量方向计算和局部平面波分解是两个最重要的内容,它们共同决定角度域成像条件的实现效率和角度域成像道集的质量.为了完善现有的角度域成像道集生成方法,本文提出一种新的波矢量计算方法和局部平面波分解方法.本文先用波动方程任意时刻的柯西条件构造一个只含原波场负频率成分的柯西波场,然后根据柯西波场在时间波数域的振幅谱计算波场的波矢量方向.该方法仅在需要计算波矢量方向的时刻合成柯西波场,不需要增加额外的数据读写操作,是一种高效的波矢量计算方法.本文还以柯西波场为基础提出一种高效的局部平面波分解方法,保证角度域成像条件的实现效率.结合柯西波场和局部平面波分解方法,本文最后给出一种新的角度域成像方法.文中最后的数值实验证明该方法得到的角度域成像道集具有理想的角度分辨率,可以反映地下构造的角度照明情况.     

拟牛顿修正法解算不等式约束加权总体最小二乘问题

根据总体最小二乘准则,可以将附有不等式约束的变量误差（errors-in-variables,EIV）模型转化为标准最优化问题,并运用有效集法、序列二次规划法等优化方法求解。已有算法在涉及计算目标函数的Hesse矩阵（二阶导数）时,存在计算量较大的缺陷。针对上述问题,利用基于拟牛顿法修正Hesse矩阵的序列二次规划算法解算附有不等式约束加权总体最小二乘问题,新算法减小了计算量,可以提高收敛速度。通过实例,证明了该算法具有很好的适用性和计算效率。     

伪逆矩阵在基准点稳定性判定中的应用

基准点决定控制网的位置、尺度和方位。监测网平差中,基准的定义相比普通平差的基准有着更高的精度,对于监测网的基准稳定性,需要通过多期复测的观测值估计的位移加以检验。检测基准点稳定性,一种"绝对"的位移或接近"绝对"的位移,通过经典平差模型中伪逆矩阵来检验后续测量中基准点是否存在动点,并用平均间隙法检校模型正确性,分析伪逆矩阵网点稳定性分析与平均间隙法平差检校进行单独平差,证明伪逆矩阵来判断检测基准网中基准点是否变动方案可行。     

基于局部加密非结构化网格的三维电阻率法有限元数值模拟

总结了目前常用的结构化网格及其局限性,分析了非结构化网格对复杂地质体边界的适应性和Delaunay三角化算法及其网格加密策略,提出了一种新的局部节点加密方法,实现了对复杂模型的完全非结构化四面体全自动剖分,给出了三维直流电阻率模拟中四面体网格的质量评价标准和最优指标.计算和分析表明,数值解在点电流源及附近的奇异区精度最低,网格加密策略可以有效地减少其影响,极大地提高数值解的精度.本文提出的局部节点加密策略计算量最小,对精度的改善也优于经典方法.在精度要求苛刻或模型十分复杂时,局部体积加密策略和二次单元是高精度模拟的可靠保证.