模糊积分在海洋标准体系评价中的应用研究

采用模糊积分综合评价模型对海洋标准体系进行评价,在指标评分过程中,无选择地邀请多位专家参与评分,利用Sugeno模糊积分和Choquet模糊积分计算,求算评价指标的近似客观值来进行综合评价,将其应用到海洋观测预报及防灾减灾标准体系评价中,得到了与标准审查会专家结论基本一致的结果。     

基于参考线圈和变步长自适应的磁共振信号噪声压制方法

核磁共振地下水探测仪的灵敏度高,接收到的纳伏级磁共振探测信号极易受到强工频谐波噪声的干扰,导致信号特征参数提取的准确度降低,影响反演解释的水文地质参数结果.为了解决这一难题,基于相关抵消的原理,针对全波磁共振信号,设计带有参考线圈的90°移相自适应噪声抵消系统,理论计算了参考线圈相对于探测线圈的距离,提出变步长LMS算法进行噪声压制.仿真结果表明,在不同的信号强度及不同的信噪比下,当信号与工频谐波干扰频谱不重合时,采用设计的自适应噪声抵消系统和变步长算法,信噪比可以提高到5.94 dB以上,初始振幅、弛豫时间特征参数的拟合误差在2.8%以内;当信号与工频谐波干扰频谱重合时,采用双向自适应滤波算法,信噪比可以达到5dB以上,初始振幅、弛豫时间特征参数的拟合误差在10%以内,可以满足实际应用的要求;实测数据处理进一步证明了方法的有效性.     

地震散射波的高精度数值模拟与振幅分析

随着地震勘探要求和技术的不断提高,散射波地震勘探方法逐渐引起人们的重视,尤其在构造复杂地区的高精度勘探中,散射波理论具有广泛的应用前景.本文在讨论了现阶段几种散射波数值模拟方法的利弊后,设计了一种基于微扰论的FK域积分法,在散射场的二次震源和空间能量衰减处理两方面做了改进,其计算精度和效率均得到了提高.此外,采用该算法针对散射波动力学进行了正演模拟试算,从扰动介质空间分布和扰动量两个角度定量地分析了散射波振幅动力学方面的特征.在此基础上得到了五个经验常数,并给出了对应的计算公式和物理意义,通过这五个常数可以划分散射波不同形态的动力学特征.     

地震分析软件震相走时表文件结构研究

通过对现有地震分析软件所配置的震相走时表文件的结构类型及优劣进行分析,将震相走时表分为4种组合结构并分析其实用性。在IASPEI91震相走时表文件结构的基础上进行修改、补充和完善,从而提出了一种安全性高、适用面广、灵活性强、能存储各种结构走时表的新型走时表文件结构。应用这种新型结构震相走时表,地震分析软件可以灵活配置不同的走时表,实现精细定位,并为数字化测震台网快速利用测震资料编制和应用高精度的区域震相走时表创造了条件,对于地球深部构造研究具有一定的促进意义。     

大地电磁测深三维正演模拟及异常特征研究

为了实现大地电磁测深三维高效快速正演算法,本文研究了大地电磁测深积分方程三维正演模拟的理论与方法.首先详细地推导了Fredholm型积分方程;然后给出积分方程的离散方法、并矢格林函数积分的求解方法以及大型矩阵方程的存储与求解方法;最后,通过理论地电模型的计算检验了方法和程序编制的正确性,并展示了几种理论地电模型的数值计算结果,分析其异常特征.算例结果表明,该正演算法在CPU主频为1.81GHz、内存为1GB的计算机上能正确地模拟多个异常体的异常特征;地表电性不均匀体会使地表视电阻率曲线产生十分严重的畸变;TE模式相对TM模式纵向分辨率比较好,而TM模式的横向分辨率较好;随着异常体埋深的增加,反映异常的分辨能力进一步降低.     

GM(1,1)模型的改进及其在变形预测中的应用

为了提高GM(1,1)模型预测精度,采用积分优化、二次拟合优化以及残差改化方法,分步对GM(1,1)模型进行改进,建立灰色多重修正模型.具体改进步骤为:首先,利用积分优化方法对背景值进行纠正,减小模型误差并提高预测精度;接着,对模型参数(发展系数和灰作用量)进行二次拟合优化,使参数更加接近理论真值;然后,根据预测结果进行适当的残差改化,提高模型整体的预测精度;最后,建立根据GM(1,1)模型改进的灰色多重修正模型.以重庆南川地区甄子岩崩塌为例,建立灰色多重修正模型对危岩裂缝累计位移值进行模拟和预测,并与GM(1,1)模型进行对比.精度检验结果表明:灰色多重修正模型后验差比值(0.082 39)明显好于GM(1,1)模型(0.192 67),平均相对残差比(0.073 9)更远好于GM(1,1)模型(0.259 6),表明灰色多重修正模型在预测精度上有较大提高,可靠性更好.     

黄土高原面积-高程分析及其侵蚀地貌学意义

黄土高原是中国乃至世界最为著名的水土流失区之一,黄河的泥沙主要来源于此.面积-高程积分(Hypsometric integral,HI)分析常被用于侵蚀地貌区的地貌发育阶段判定.通过探讨HI与DEM分辨率、分析窗口之间的尺度效应,分析了黄土高原HI空间格局和发育阶段性,并在此基础上揭示了侵蚀产沙和地貌发育关系.研究表明:利用不同分辨率DEM计算的HI具有很好的稳定性,HI对DEM分辨率不具有依赖性;然而,利用不同分辨率的分析窗口所计算的黄土高原原HI值,表现出随着分析窗口的增大,HI表现出按照乘幂函数关系下降的趋势.黄土高原HI具有显著的空间分异特征,即汾渭谷地、宁夏内蒙沿黄冲积平原等地堑凹陷区是HI低值区,而黄土丘陵沟壑区是HI的高值区;黄土高原丘陵沟壑区地貌发育处于幼年期、壮年期和老年期的面积分别占总面积的7.9%、83.4%和8.7%,整体来看处于地貌发育壮年期;黄土高原丘陵沟壑区侵蚀产沙强度和地貌发育过程关系密切,HI与输沙模数存在显著的正相关.当前,黄土高原强烈的侵蚀产沙是其地貌发育处于壮年期的自然表现,这将是一个非常漫长的自然过程.     

非饱和渗流问题的自适应欠松弛变量变换方法

在使用有限元方法求解非饱和土渗流问题时,土-水特征曲线和渗透率函数的强烈非线性经常会造成计算中出现迭代不收敛、计算误差大等问题。基于变量变换的思想,结合时间步长自适应技术提出了一种求解非饱和渗流问题的新方法--欠松弛RFT变换方法（ATUR1）。ATUR1方法通过变量变换,大大降低了Richards方程中未知数在空间和时间上的非线性程度,从而改善这种非线性所带来的计算收敛困难和精度差等问题。欠松弛技术的引入减少了迭代过程中的振荡现象,进一步提高了非线性迭代计算的效率。时间步长自适应技术则有效地控制整个计算过程的误差。数值算例结果说明,ATUR1可以有效地提高计算效率和精度,是一种准确有效的计算方法。     

亚塑性模型不同积分算法的数值实现

亚塑性模型是以Jaumann应力率张量及变形率张量描述的一种率型本构关系,本构关系在非线性有限元分析计算中具有关键作用,解决用应变增量求解应力响应的问题需要一个时间积分过程。针对亚塑性本构模型发展了自适应隐式和显式两种不同的积分算法,给出了误差控制的方法,同时推导了自适应隐式积分算法所需的一致切向模量,并采用了两个不同的单元,利用ABAQUS平台比较了两种积分算法的数值模拟结果。为了实现从ABAQUS/Standard到ABAQUS/Explicit的过渡,开发了UMAT-VUMAT接口,从而可以使已有的UMAT子程序用于大变形动力问题分析。算例分析证明了研究结果的正确性。     

适于Kirchhoff叠前深度偏移的地震走时李代数积分算法

本文基于拟微分算子理论和李代数积分法,根据程函方程和波场坐标变换,提出一种新的适于横向变速介质Kirchhoff叠前深度偏移的地震波走时算法.该算法与Kirchhoff叠前时间偏移所用李代数时间积分表达相比,差异在于增加了波数一次项,且二次项的系数在求积时亦需进行修正.针对单平方根算子象征、李代数积分、指数映射和走时多项式的求解而言,皆需对以往Kirchhoff叠前时间偏移中所用算法进行深化调整.文中数值算例对比了本文李代数积分表达与时间积分的区别,本算法计算结果与线性横向变速介质中的理论值相当吻合.通过走时多项式中各项对结果的影响分析,可知非对称项使计算精度得到了进一步提高.数值试验表明,本算法对横向变速介质中走时求取是可行的,且不需要存储海量走时表,有利于提高Kirchhof叠前深度偏移的精度和效率.