用改进的逐次逼近解法计算二维井间电磁场

提出了一种计算积分方程的新方法——改进的逐次逼近解法 (MSAM) ,并用该方法计算了轴对称二维井间电磁场。与传统的逐次逼近解法 (SAM)相比 ,该方法收敛性强 ,应用范围广 ,可适用于高电导率对比地层。由于不必进行直接的大型矩阵求逆运算 ,因此与积分方程的直接解法 (IE)相比 ,该方法计算速度更快 ,所需内存更少。采用矩阵求逆方法计算了成层介质中的二维 Green函数 ,并对含 Green函数的积分进行了简化 ,从而加快了计算速度。数值计算结果显示 ,在地层电导率对比度达到 2个数量级时 MSAM仍收敛 ,且计算结果与直接求解积分方程的结果一致 ,因此

直流电测井响应的二阶几何因子逼近计算

张庚骥等：直流电测井响应的二阶几何因子逼近计算，测井技术，１９９６（５）２０，３２０～３２４。对直流电测井几何因子理论及逐次逼近解法进行了进一步的推导论证，同时给出了电极条件下二阶直流电几何因子逼近的具体计算方法和数值分析结果。计算结果显示，较一阶几何因子逼近结果有明显改进，已经基本接近真值，进一步证明了直流电几何因子理论在电测井数值模拟方面的有效性。

利用积分方程的加速迭代算法计算随钻电磁波电阻率测量仪器的响应

采用改进型逐次逼近解法(MSAM)和Aitken加速技术相结合的迭代算法计算二维积分方程.该算法将纵向成层原状地层作为背景地层,将计算区域限制在井眼和侵入带内,具有未知量数目少、收敛速度快、计算精度高的优点.利用该算法对随钻电磁波电阻率测量仪器在轴对称二维地层中的响应进行数值模拟.模拟结果显示,幅度衰减曲线和相位移曲线受井眼、侵入和围岩的影响程度不同,二者径向探测深度和垂向分辨率也有差异,利用补偿后的幅度衰减电阻率曲线和相位移电阻率曲线交叉点的坐标可精确确定地层层界面位置.

二维井间电磁场的正反演计算

提出了一种改进的逐次逼近解法 ,并用该方法对轴对称二维井间电磁场进行了正反演计算。与逐次逼近解法相比 ,该方法收敛性强 ,应用范围广 ,可适用于高电导率对比地层。与直接求解积分方程相比 ,由于不必进行直接的大型复矩阵求逆运算 ,该方法计算速度快 ,所需内存少。反演中采用了Born迭代算法 ,将成像区域集中于一定范围内。考虑到信息量和计算机内存的限制 ,采用双重面元分割法 ,将反演区域分割为一定数目的大面元 ,每个大面元被分割为几个小面元。属于同一大面元的不同小面元具有不同的磁矢势 ,但具有相同的电导率 ,从而减少了反演过程中未知量的数目。根据第一次成像结果确定出更准确的成像范围 ,并进行第二次成像。数值计算结果表明 ,改进的逐次逼近解法是计算二维井间电磁场的一种有效方法 ,将该方法用于反演过程能够得到较高分辨率的二维井间电导率图像。

利用积分方程计算阵列感应测井响应

利用积分方程计算了阵列感应测井仪器各子阵列在轴对称二维地层中的响应,计算中将纵向成层原状地层作为背景地层,从而将计算区域限制在井眼和侵入带内。与数值模式匹配法和有限元素法将全部地层作为求解区域相比,未知量范围和计算量均大大减小,计算精度亦可相应提高。采用改进型逐次逼近解法(MSAM)计算了积分方程,该方法所需内存量小,且可适用于高电导率对比度地层。采用递推算法计算纵向成层背景地层中Green函数及其对空间坐标积分的待定系数,从而大大加快了计算速度。数值算例的结果说明了利用积分方程计算阵列感应测井响应的有效性。

张庚骥先生的《电测井算法》

张庚骥先生从事电法测井的研究和教学工作50余年，矢志不渝地钻研电测井方法，培养了一大批测井技术骨干，其中不乏在国内外享有盛名的电法测井专家、教授。中国工程院孙龙德院士高度评价张先生的研究成果：“首次提出了计算非均匀介质电场的逐次逼近解法，……，在国际上被命名为张氏几何因子，大大提升了中国测井界的学术地位。”2009年，在张先生80岁高龄之际，石油工业出版社出版了由他独立著作的《电测井算法》一书。正如中国石油学会石油测井专业委员会主任陆大卫先生在该书《序》中所言，“《电测井算法》集中反映了张庚骥教授电法测井理论和方法的精髓”。为了让更多的电法测井人员了解和读懂张先生的研究成果，用这些理论指导科研和生产实际，我们请张先生的弟子、中国石油大学（北京）的高杰教授为此书写了一个书评，供参考。

一种气动柔性机械手工作的模型研究

提出一种波壳伸缩式的弯曲关节,及其构成的机械手。此关节采用板弹簧为柔性骨架,波壳弹性体受气压后轴向膨胀作为肌肉动力。文中给出了关节的结构,着重分析了气体压力与弯曲变形参数的静态模型。由此导出悬臂端在受力偶、法向力和一大小、方向不断变化的切向力的共同作用下,悬臂梁大挠度变形的数学模型,并采用网格与逐次逼近凑解法求得模型的数值解。