用改进的局域非线性迭代方法计算三维井间电磁场

提出了用改进的局域非线性迭代（MLNI）计算体积分方程的方法，并计算了三给井间电磁场。将井间大尺度电导率异常体分子为近场区域和远场区域两部分，它们的位置和尺寸均随场点位置的变化而变化。对近场区域影响的计算采用局域非线性近似；将远场区域的影响作为外部激励源，采用迭代方法进行计算。由于不必进行直接的大型矩阵求逆运算，因而与体积分方程的直接解法相比，所需的机时更少，并减少了对内存的要求。数值计算结果显示，该方法也适用于高电导率对比地层，且计算精度与体积分方程直接解法的精度相当，是一种计算井间大尺度异常散射场的有效方法。     

采用轴对称二维模型对实际地层井间电磁响应进行成像的精确性研究

采用二维积分方程和体积分方程分别模拟轴对称二维井间地层和三维井间地层模型中异常体的响应,并分析三维异常体沿垂直于发射—接收剖面方向的厚度和偏心程度以及在发射—接收剖面的横向位置对其响应的影响。采用Born迭代方法,通过反演不同模型的模拟数据获得发射—接收剖面的二维电导率成像。结果表明,在将三维地层中得到的井间测量数据在发射—接收剖面采用轴对称二维模型进行成像时,其成像质量与异常体沿垂直于发射—接收剖面方向的厚度和偏心程度直接相关。只要异常体沿垂直于发射—接收剖面方向的厚度足够大且其中心点距离发射—接收剖面较近,就可以采用轴对称二维井间地层模型实现高质量成像。异常体厚度越小、偏心程度越高,成像效果越差。     

三维井间电磁场的联合反演方法

考虑到井间电磁测量信息量的局限性,基于积分方程在发射接收剖面对三维电导率异常体进行近似成像。成像过程中将沿垂直于发射接收剖面方向的成像区域限制在一定范围并假设电导率在该方向不变。成像实例显示,这种剖面成像的精度取决于垂直于发射接收剖面方向的成像范围与异常体在该方向实际范围的接近程度。考虑到实际情况下背景层电导率值往往是未知的或虽已知但并不准确,还提出了联合反演方法,以一维反演结果作为初始值,既反演一维背景地层的电导率,也反演三维异常体在发射接收剖面的电导率,两个过程交替进行。反演实例显示,经过几轮迭代后,背景层电导率分布及异常体在发射接收剖面的电导率分布均大大改善,说明了联合反演的有效性和必要性。     

三维井间电磁场的联合反演方法

考虑到井间电磁测量信息量的局限性，基于积分方程在发射-接收剖面对三维电导率异常体进行近似成像。成像过程中将沿垂直于发射-接收剖面方向的成像区域限制在一定范围并假设电导率在该方向不变。成像实例显示，这种剖面成像的精度取决于垂直于发射接收剖面方向的成像范围与异常体在该方向实际范围的接近程度。考虑到实际情况下背景层电导率值往往是未知的或虽已知但并不准确，还提出了联合反演方法，以一维反演结果作为初始值，既反演一维背景地层的电导率，也反演三维异常体在发射接收剖面的电导率，两个过程交替进行。反演实例显示，经过几轮迭代后，背景层电导率分布及异常体在发射-接收剖面的电导率分布均大大改善，说明了联合反演的有效性和必要性。     

二维井间电磁场的非线性反演方法

采用高精度的非线性反演方法,对轴对称二维井间模型的电导率分布进行了成像。该成像方法既考虑到雅可比矩阵的线性项,也考虑了其偏微分项,因而反演精度高且稳定性强。采用矩阵求逆技术计算了纵向成层背景地层Green函数的系数,可以很方便地得到Green函数对空间坐标的积分,从而使雅可比矩阵的计算更方便、准确。反演实例表明,井间电磁成像技术更适宜于对相对高电导率异常区域进行成像。由于井间电磁成像技术对异常区块的垂向探测范围大于径向探测范围,测量数据中包含的垂向信息多于径向信息,因而垂向分辨率高于径向分辨率。     

地震波旅行时非线性/线性联合插值法

地震波旅行时计算是射线追踪方法的核心问题,其精度直接影响着路径的精度和正演模拟的准确程度.为了提高地震波旅行时的计算精度,基于Asakawa的旅行时线性插值法(LTI)提出了一种改进的射线追踪法--非线性/线性联合插值法,即在震源的近场采用非线性插值计算旅行时,在远场仍用LTI算法.非线性/线性联合插值法提高了近场旅行时的计算精度,进而提高了全场(尤其是远场)旅行时的计算精度.数值模拟和模型试算的结果表明,本文提出的非线性/线性联合插值法较传统的LTI算法更为精确.     

一个大尺度和非线性积分方程的数值解法

讨论了在t夸克凝聚模型中,利用梯形近似下的Schwinger-Dyson方程计算t夸克动力学质量的方法。研究了这个大尺度和弱非线性的积分方程的“精细调节”问题,并给出了一种新的迭代方法。     

二维井间电磁场的非线性反演方法

采用高精度的非线性反演方法，对轴对称二维井间模型的电导率分布进行了成像。该成像方法既考虑到雅可比矩阵的线性项，也考虑了其偏微分项，因而反演精度高且稳定性强。采用矩阵求逆技术计算了纵向成层背景地层Green函数的系数，可以很方便地得到Green函数对空间坐标的积分，从而使雅可比矩阵的计算更方便、准确。反演实例表明，井间电磁成像技术更适宜于对相对高电导率异常区域进行成像。由于井间电磁成像技术对异常区块的垂向探测范围大于径向探测范围，测量数据中包含的垂向信息多于径向信息，因而垂向分辨率高于径向分辨率。     

层状介质二维井间电磁成像

根据分层介质中波的传播规律,利用均匀介质中的三维张量格林函数推导出二维情况下分层介质的张量格林函数,它是一个对kx的一维Sommerfeld积分。准确求取这一积分后,可以方便地利用迭代等方法反演出井间的电导率。给出了利用上行波和下行波方式表达的界面反射和折射张量函数及其各分量。通过积分方程方法对埋置于有8个分层的两个垂直方向电导率异常体模型进行了计算,反演结果与原模型基本吻合,有效地实现了层状介质中二维井间电磁成像。     

层状介质二维井间电磁成像

根据分层介质中波的传播规律，利用均匀介质中的三维张量格林函数推导出二维情况下分层介质的张量格林函数，它是一个对kx的一维Sommefield积分。准确求取这一积分后，可以方便地利用迭代等方法反演出井间的电导率。给出了利用上行波和下行波方式表达的界面反射和折射张量函数及其各分量。通过积分方程方法对埋置于有8个分层的两个垂直方向电导率异常体模型进行了计算，反演结果与原模型基本吻合，有效地实现了层状介质中二维井间电磁成像。     

局部Petrov-Galerkin方法在非线性边值问题中的应用

把一种真正的无网格局部Petrov-Galerkin方法用于求解非线性边值问题.为了克服一般局部Petrov-Galerkin方法计算工作量较大的问题,选择一个分段函数作为加权残值法的加权函数,简化了非线性问题中刚度矩阵的域积分.基于局部Petrov-Galerkin积分方程逐点建立的思想,推导了一种直接插值法用于施加本质边界条件.通过算例表明,这种局部Petrov-Galerkin方法是一种具有收敛快、精度高的方法.     

快速多极边界元方法在二维声散射问题中的应用

快速多极算法(FMM)是求解边界元方法(BEM)在大尺度情况下的一种非常有效的算法.研究了快速多极算法在二维声散射问题的边界积分方程求解中的应用.给出了积分核函数以及其共轭积分算子核函数的多极展开式,局部展开式以及相应展开系数之间的转化关系.分别应用两种不同的层级树结构的FMM来进行求解,并对两种树结构下的求解效率进行了对比.数值算例表明用快速多极算法求解该问题时在存储量和计算量上比直接求解方法效率更高.     

二维流场模拟的矩阵标识法

采用二分步法，从积分型方程出发，在有限控制体上建立守恒型差分格式，对二维浅水波方程进行求解。求解线性代数方程组采用矩阵标识法，提高了计算效率，特别适合于流域大尺度水流的模拟。计算网格采用无序编码方法，灵活方便。算例表明，本文方法能够模拟宽广湖泊一类大尺度的二维水流。     

新的函数变换法与一类非线性波方程的精确孤波解

采用新的函数变换法求出了一类非线性演化方程的两类显示精确孤波解.作为该方程的特例,如K le in-Gordon方程、Landau-G inburg-H iggs方程、Duffing方程和4方程等也都获得了相应的精确孤波解.这种方法也适用于求解具有更高次非线性项的其他非线性波方程.     

非线性热传导问题的基于滑动Kriging插值的MLPG法

利用基于滑动Kriging插值的无网格局部Petrov-Galerkin(MLPG)法来求解二维非线性稳态和瞬态热传导问题,Heaviside分段函数作为局部弱形式的权函数,并通过加权余量法推导相应的离散方程.该问题考虑了材料热传导系数随温度的线性变化,并通过拟线性法来求解非线性问题的解,时间域的离散通过向后差分法来实现.基于滑动Kriging插值构造MLPG中的形函数由于满足克罗内克δ性质,因此可以直接准确地施加本质边界条件.在构造刚度矩阵过程中,只涉及边界积分,不涉及区域积分和奇异积分.将数值计算结果与有限元法得到的结果加以对比可以看出,基于滑动Kriging插值的MLPG法能够很好地解决此类热传导问题.     

具有5次强非线性项波方程的精确解

采用新的函数变换法,并与直接积分法相结合简便地求出了具有5次强非线性项的导数Schr dinger方程四类显示精确孤波解。这种方法同样也适用于求解具有更高次非线性项的其他非线性波方程。     

一种求解大规模方程组的修正Dai-Yuan共轭梯度法

设计一种针对大规模非线性方程组的修正DY共轭梯度算法.该算法的搜索方向不仅自动满足充分下降条件,而且属于信赖域.在适当条件下,可以证明新算法是全局收敛的.初步的数值实验表明新算法可以有效求解大规模非线性方程组.     

随机大系统的Lyapunov函数构造与符号计算

讨论了如何利用随机第一积分为构造非线性随机动态系统的Lyapunov函数的基本原理与方法,并在这个基础上建立孤立子系统的第一积分与随机大系统Lyapunov函数的关系,然后利用这个基本原理为构造随机大系统的备选Lyapunov函数,这一方法即可以得到一般随机系统的Lyapunov函数的结构,又便于作计算机的自动推导与符号计算。     

分析三维随机介质目标散射问题的SMCG方法

应用稀疏矩阵规则网格(sparse matrix canonical grid，SMCG)法分析了三维随机介质目标的电磁放射问题.用矩量法术解介质放射体的体积分方程时，根据放射体离放单元间场相互作用的强弱，将阻抗矩阵分解为近区强相互作用的稀疏矩阵和远区弱相互作用的补充矩阵、在用共轭梯度法迭代求解矩阵方程时，将格林函数在规则网格点上进行泰勒级数展开，进而可利用快速傅里叶变换计算弱相互作用矩阵与待求向量的乘积，而强相互作用矩阵与待求向量的乘积可以直接计算、文中对几种不同情况的随机介质目标的远区放射场进行了计算，结果表明SMCG法的计算结果与满阵矩量法的计算结果吻合良好，而所需的计算机内存和计算量却大为减少．