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

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

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

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

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

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

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

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

地面一井孔电磁系统成像方法研究

在地面一井孔（ＳＴＢ）电磁成像系统中,将发射器放置在地面距裸眼井一定距离和方位处,将接收器置人裸眼井的目的层段。使用这种电磁系统,只需一口裸眼井,现场测量灵活方便。在成像过程中,应用基于逐次逼近解法的迭代反演算法对目的层段的二维电导率剖面进行成像。该方法用一阶Ｂｏｍ近似将积分方程线性化,在选代反演中应用高阶Ｂｏｍ近似。为了提高成像效果,采用多种发射频率并扩大接收器的垂向接收范围和接收密度,但是由于发射源固定在地面,成像效果低于跨井孔电磁成像效果。数值计算结果表明,ＳＴＢ电磁成像系统能较准确地确定出电导率异常体的垂向位置,但水平方向分辨能力差,精确地描述地层中电军率异常体的分布还需要其它测量信息。     

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

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

井间电磁场的一维反演

开发了一种快速高效的井间电磁场一维反演方法。通过矩阵求逆计算纵向成层地层中的Green函数并同时得到雅可比矩阵的元素，加快了计算速度。反演时采用正则化最小二乘法求解超定方程，由于发射－接收的组合数目大大超过待反演的地层参数的数目，因而反演结果精度高。数值计算举例表明，无论是低电导率对比地层还是高电导率对比地层，反演结果均与模型完全吻合。井间电磁场一维反演可以对实测数据进行处理，为高维成像奠定基础。     

基于开放图形库的储层剖面图的绘制

基于开放图形库(OpenGL)技术,实现了垂直于地层走向沿x或y方向以及任意方向的储层剖面图的绘制.首先确定剖切面及被剖切体,并由鼠标确定视图工区中的某一点,再通过该点确定二维平面上的直线方程,有了直线方程就可以确定与地层走向垂直的储层剖面.这样就把削切面和储层体求交面的复杂问题转换成剖切面和多个储层层面求交线的简单问题,计算出所有层面的所有网格线与剖切面的交点集就构成了要求的剖切面.以某油田106井区S1和S2小层间沉积相的剖面图为例,验证了储层剖面图绘制的正确性与直观性.     

大地电磁资料三维傅里叶变换偏移方法及其应用

利用傅里叶变换及内插外推方法对三维大地电磁资料偏移成像问题进行了研究。首先将大地电磁场满足的三维亥姆霍兹方程降维分裂成两个二维响应方程，然后对每一个二维方程进行傅里叶变换，使二维方程变为一维方程，通过大地电磁场垂向方向的导数在地面上所满足的边界条件，确定外推方程中的两个待定系数，从而得到一维外推方程的解。对横向非均匀地电介质，利用分段傅里叶变换方法进行偏移，同时将每一段内电阻率的平均值作为背景值。为了解决纵向介质的不连续问题，取相邻两频点的电阻率值的外推结果进行加权内插。对于横向各段的衔接，采用内插平滑处理。应用该方法对伊朗项目中T01—010测线资料进行了偏移处理，偏移结果除了反映浅层和基底明显的分界面之外，还能较好地反映局部地电异常体引起的异常。     

轴对称异质体引起波散射问题的边界元解

通过将全空间内轴对称异质体边界载荷沿周向作Ｆｏｕｒｉｅｒ级数展开，然后利用边界元法在母线方向求解边界积分方程，使弹性波散射问题的维数由三维降到一维．该方法充分利用旋转体轴对称的几何特点，采用环壳单元使计算量较普通边界元法更小，并且收敛速度较其他边界元法更快．     

螺旋芯棒式吹膜模头熔体流动特性的数值模拟

对吹塑薄膜螺旋芯棒式模头内熔体的流动现象进行了三维数值模拟。结果表明,随着螺槽的逐渐消退,螺槽内熔体的速度由入口处的沿螺槽方向流动逐渐变为沿模头轴向流动,而环形间隙内熔体的速度在入口处存在较小的环形分量,随间隙厚度的增加逐渐变为沿轴向流动。沿挤出方向,螺槽中熔体的压力等值线由垂直于螺槽逐渐变为垂直于轴向,环形间隙中熔体的压力等值线由倾斜于轴向逐渐变为垂直于轴向。     

变截面波纹钢腹板连续刚构桥拟平截面假定试验研究

设计制作了3跨单箱单室变截面波纹钢腹板连续刚构桥试验模型,结合试验和有限元数值模拟,对变截面波纹钢腹板连续刚构桥的"拟平截面假定"进行研究.结果表明:当波纹钢腹板箱梁受弯矩作用时,波纹腹板的轴向抗力可以忽略不计,弯矩基本都由顶底板共同承担;试验梁中跨3/8、1/4截面沿梁高度纵向应变分布趋势与中跨1/2截面基本一致;中跨支座截面沿梁高度纵向应变方向与中跨1/2、3/8、1/4截面相反,且顶底板应变基本相同;中跨1/2截面和中跨支座截面应变分布满足"拟平截面假定",中跨3/8、1/4截面不符合该假定,但从中跨1/4截面到1/2截面的应变分布趋于吻合该假定.     

帽形加劲复杂卷边槽钢柱的稳定性能

为进一步了解帽形加劲复杂卷边槽钢柱的受力性能,利用ANSYS有限元分析软件对帽形加劲复杂卷边槽钢柱在轴压和偏压状态下的稳定性能进行分析.研究了立柱长度、板件宽厚比、截面类型以及偏心距对帽形加劲复杂卷边槽钢柱的稳定承载力、破坏模式及变形等性能的影响.研究结果表明: 在板件厚度不变时,C2截面构件的承载效率相比C1截面构件提高了10%~20%,C3截面构件相比C1截面构件提高了40%~50%.腹板帽形加劲复杂卷边槽钢柱的承载力最大发生在负偏心一侧,而腹板及翼缘帽形加劲复杂卷边槽钢柱在正偏心处承载力最大.     

随钻方位电磁波仪器补偿测量方法研究

基于层状各向异性电磁场的解析解,对带有横向天线的方位电磁波仪器进行快速响应数值模拟。以胜利钻井工艺研究院的多频方位电阻率(AMR)仪器和贝克休斯方位电磁波电阻率(APR)仪器为例,对采用轴向发射、横向接收的测量单元分别进行单发单收、单发双收、双发单收以及双发双收的测量组合响应模拟,根据模拟结果对补偿测量的效果进行评价分析。模拟结果表明,单发单收测量单元能实现基本的方位电磁波的测量功能,但易受到地层环境的影响,使其应用受到一定限制。通过发射或接收补偿可以不同程度地消除仪器偏心、相对倾斜角以及电阻率各向异性的影响,使方位电磁波测量资料更好地应用于地质导向中。     

海洋石油立管的三向大位移分析

本文对截面具有两个不相等主惯性矩的海洋石油立管的三向大位移分析作了理论推导。在推导中,全面考虑了立管弯曲、扭转和拉伸变形。通过综合立管轴向的平衡条件和能量泛函的驻值条件,得出了独特的混合方法。平衡状态时立管中心线的孤长s取为自变量,而立管中心线的坐标X、Y和截面扭转角F以及管内拉力T则为孤长s的函数。在上述理论推导的基础上,形成了求解具有二相等主惯性矩的立管三向六位移分析的有限元方法。根据此方法编写了计算机程序,并给出了算例。     

水平井段岩屑床厚度模式的建立

基于环空固液两相流水动力学理论建立了水平井段岩屑床厚度的理论计算模式.采用正交实验设计来安排实验.在大量实验基础上.分析了环空流体返速、钻具旋转速度、钻具与井眼间偏心度、岩屑注入速度、泥浆密度、泥浆流变性及岩屑尺寸等因素及其交互作用对无因次岩屑床厚度的影响.运用多元统计分析的方法建立了综合上述因素的岩屑床厚度经验模式.并对两个模式进行了验证.结果表明。两模式都可用于预测水平井段岩屑床厚度.使用经验模式求解了环空不淤积临界返速.     

螺旋桨初生空化湍流的多相流数值模拟

摘要：
同时采用修正剪切应力输运(SST)湍流模型和Baseline雷诺应力模型（RSM）求取了E779A螺旋桨在无空化状态和初生空化状态下的梢涡运动轨迹,分析了涡核最小压力系数、湍动能、轴向速度分量和涡核半径沿运动轨迹的变化,并从模拟得到的梢涡卷曲起始和梢涡涡束的角度阐述了梢涡形成机理.空化模型采用改进Sauer模型,考虑了非凝结性气核质量分数、体积分数和气泡初始半径以及湍流脉动的影响,并针对轻度、中度和重度空化面积进行了可信性校验.当空化数σ>初生空化数σi时,叶梢截面压力系数分布相对不再改变的判定准则来确定.涡核中心位于螺旋线垂向截面上最小压力点,涡核边界由湍流涡频率峰值决定.数值模拟结果表明,RSM模拟梢涡路径较修正SST湍流模型稍长、局部梢涡空化范围略大、叶梢最小压力系数和轴向速度分量要小,涡核湍动能分布更为合理.但两者模拟得到的涡核运动轨迹几乎重合,并且初生空化状态下的涡核运动轨迹、最小压力系数和轴向速度分布均与各自无空化状态下非常接近,表明了初生空化状态判定的正确性和改进数值模型对梢涡运动轨迹模拟的适用性.     

序列切片图像中斑块状污染的消除

为消除序列切片显微图像中的斑块状污染噪声,提出了一种三维截面图像顺序滤波的处理算法。根据斑块状污染噪声随机出现在单个层片的特点,通过数据重构,得到三维图像体数据,并针对斑块状污染噪声在不同观察方向上的形态退化,分别进行特定方向的中值滤波。算法可有效消除原始图像中的斑块状污染噪声,使图像分割和三维图像重构组织的连通关系得以保证,该算法等价于采用特殊的序列三维中值滤波器对图像进行处理。