电容层析成像系统图像重构算法
————————————————————————————————作者：————————————————————————————————日期：
电容层析成像系统图像重构算法
地震层析成像理论及技术-瑞雷面波理论基础与反演成像
瑞雷面波理论基础与反演成像 瑞雷面波是1887年由英国学者瑞雷（Rayleigh ）首先在理论上确定的，这种面波分布在自由表面上。当介质为均匀各向同性介质时，瑞雷面波的相速度和群速度将一致，否则瑞雷波的相速度将不一致，出现频散现象，当介质具有水平层状性质时，瑞雷面波的频散规律与介质的分层结构紧密相关。面波研究的目的是要通过面波信号得到地下介质的结构及其物理力学性质，这就需要进一步反演解释研究。 1. 瑞雷波的理论基础 由于均匀弹性半空间介质的边界附近沿x 方向传播的平面瑞利谐波y 方向的质点位移为零。设半空间充满x-y 平面，z 方向向下为正，坐标原点位于介质的自由表面上，如图所示1-1 为推导方便，引入势函数Φ和ψ来分别表示x 和z 方向的位移（u 和w ），则 ，u w x z z x ?Φ?ψ ?Φ?ψ = -= +???? 1.1 平面瑞利波波前 质点位移随深度增加 而衰减 波的传播方向 图1-1 均匀弹性半空间中的平面瑞利波
由位移表示的二维运动方程为 2222 22u u w w （）（）ερλμμερλμμ??=++?????=++???t x t z 1.2 由此可见，势函数的引入将胀缩波和剪切波区分开来（Φ与胀缩波对应，ψ与剪切波对应）。将式（1.1）代入（1.2）得 22222222 222222x t z t x z z t x t z x ρρλμμρρλμμ??????Φ??ψ?? -=+?Φ?ψ ? ?????????????????Φ??ψ??+=+?Φ+?ψ ? ???????????（）（）-（）（）（）（） 1.3 又有 22222222p s 22 2v v ，λμμρρ ?Φ+?ψ=?Φ=?Φ=?ψ=?ψ??t t 1.4 由于平面瑞利波的位移发生在x-z 平面内，因此由式（1.1）和式（1.4）可知，瑞利波是P 波和SV 波相互作用的结果。 对于一个角频率为ω，波数为k ，沿x 方向传播的瑞利谐波，其势函数可表示为： ()()F ()G ()，ωω--Φ=ψ=i t kx i t kx z e z e 1.5 其中，F()z 和G()z 分别表示瑞利波胀缩分量和旋转分量的振幅随深度变化的函数；波数R 2L k π = ，R L 为瑞利波波长。 将式（1.5）代入式（1.4）并整理得 22222p 2 2 222s F()F()=0v G()G()=0v ωω?? ?-- ? ?????? ?-- ??? ? z k z z z k z z 1.6 上述二阶偏微分方程的通解为 1122F()=A B G()=A B --++qz qz qz qz z e e z e e 1.7
电容层析成像技术在多相流测试中的应用和进展 摘要 电容层析成像技术是最早发展起来的过程成像技术之一，因其具有非侵入式的测量方式、成本低、安全性好、易于实现等优点，具有极好的工业应用前景，获得了很多学者的研究并取得丰硕成果。本文简单介绍了电容层析成像系统的原理和发展历史，总结了研究热点和现状，指出了发展中的存在的问题和发展方向。 关键词：电容层析成像；电容传感器；图像重建 引言 电容层析成像技术（ECT）是根据被测物质各相具有不同的相对介电常数，在流体流动管道的某一截面上，沿管道壁均匀地布置敏感阵列电极，任意两个电极板均可组成一个两端子电容，当各相组分分布或者浓度发生变化时将引起混合流体等价相对介电常数发生变化，从而使测量电极对间的电容值发生变化。采用阵列式电容传感器，各电极之间相互组合的多个测量值便可反映出多相流体的浓度以及在管道截面上的相分布情况，以这些电容测量值作为投影数据，通过一定的图像重建算法，便可重建出被测物场的浓度分布图像。典型的ECT系统由电容传感器阵列、数据采集系统、图像重构计算机三部分组成。电容成像系统通过传感器阵列从不同方向获得管道内介质分布状况的投影信息；测量及数据采集系统收集电容传感器的输出电容值，并进行滤波、变换、放大，将数据通过接口传给成像计算机；成像计算机通过重建算法重构管道横截面的相分布图像，同时向测量及数据采集系统传送控制信号[1, 2]。目前，电容层析成像的研究多围绕这几部分展开，即电容测量电路研究、电容传感器研究及图像重建与数据处理研究。 该技术在上世纪80年代末由英国曼彻斯特大学理工学院（University of Manchester Institute of Science and Technology，简称UMIST）提出的一种新的计算机层析成像技术[3]，主要用于工业管道内的多相流检测，可提供常规仪器无法探测的封闭管道及容器中多相介质的浓度、分布、运动状态等可视化信息，与其他测量技术或仪表相配合还可应用于多相流总质量流量、分相质量流量以及流速的实时检测。比之其它技术，电容层析成像技术具有适用范围广、非侵入式、安全性能好等优点[4-6]，各种工业生产过程中常见的多相流均能应用该技术，并且 成本低廉，更适合我国的国情。电容层析成像技术为从根本上解决多相流的多参数可视化测量问题提供了一条较好的途径，是实现多相流参数在线检测的一种理想手段[7]。因此，完善和发展电容层析成像技术理论是一项很有意义的研究工作。近年来，国内外的研究人员积极探索，做了大量的工作，电容层析成像技术的研究已经取得了很大的进展，理论研究方面也不断推进，取得了一些新的研究成果[8-13]。
技术创新 中文核心期刊《微计算机信息》（管控一体化）２００７年第２３卷第６－３期 ３６０元／年邮局订阅号：８２－９４６ 《现场总线技术应用２００例》 图像处理 电容层析成像系统图像重建算法的研究 ＳｔｕｄｙｏｎＩｍａｇｅＲｅｃｏｎｓｔｒｕｃｔｉｏｎＡｌｇｏｒｉｔｈｍｆｏｒＥｌｅｃｔｒｉｃａｌＣａｐａｃｉｔａｎｃｅＴｏｍｏｇｒａｐｈｙＳｙｓｔｅｍ （江南大学）曹琳琳 ＣＡＯＬＩＮＬＩＮ 摘要：本文利用Ｔｉｋｈｏｎｏｖ正则化和奇异系统理论，分析了引起电容层析成像系统逆问题不适定性的根本原因是由于敏感场 矩阵小奇异值的存在。针对一般Ｔｉｋｈｏｎｏｖ正则化方法将所有的奇异值都采取同一正则化参数修正带来的误差，本文将小奇异值对应的项设定正则化参数，而舍去零奇异值对应向量，既减少了误差又加快了速度。例算结果表明，用本文方法重建图像，比其它如线性反投影算法（ＬＢＰ）、Ｌａｎｄｗｅｂｅｒ迭代法及一般Ｔｉｋｈｏｎｏｖ正则化算法，都有一定程度的改善。关键词：电容层析成像；图像重建算法；Ｔｉｋｈｏｎｏｖ正则化；奇异系统中图分类号：ＴＰ２１２文献标识码：Ａ Ａｂｓｔｒａｃｔ：ＢａｓｅｄｏｎＴｉｋｈｏｎｏｖｒｅｇｕｌａｒｉｚａｔｉｏｎａｎｄｓｉｎｇｕｌａｒｓｙｓｔｅｍｔｈｅｏｒｙ，ｉｔｉｓａｎａｌｙｚｅｄｔｈａｔｔｈｅｃａｕｓｅｏｆｔｈｅｉｌｌ－ｐｏｓｅｄｃｈａｒａｃｔｅｒｉｓｔｉｃｏｆｔｈｅｉｎｖｅｒｓｅｐｒｏｂｌｅｍｉｎｅｌｅｃｔｒｉｃａｌｃａｐａｃｉｔａｎｃｅｔｏｍｏｇｒａｐｈｙｓｙｓｔｅｍｉｓｄｕｅｔｏｔｈｅｓｍａｌｌｓｉｎｇｕｌａｒｖａｌｕｅｓｏｆｔｈｅｓｅｎｓｉｔｉｖｉｔｙｄｉｓｔｒｉｂｕｔｉｏｎｍａ－ｔｒｉｘ．ＤｉｆｆｅｒｆｒｏｍｔｈｅｃｏｎｖｅｎｔｉｏｎａｌｔｅｃｈｎｉｑｕｅｉｎＴｉｋｈｏｎｏｖｒｅｇｕｌａｒｉｚａｔｉｏｎｍｅｔｈｏｄｔｈａｔａｌｌｔｈｅｓｉｎｇｕｌａｒｖａｌｕｅｓａｒｅｍｏｄｉｆｉｅｄｗｉｔｈｔｈｅｓａｍｅｐａｒａｍｅｔｅｒ，ｗｈｉｃｈｗｉｌｌｒｅｓｕｌｔｉｎｍｏｒｅｓｏｌｕｔｉｏｎｅｒｒｏｒ，ａｎｅｗｔｅｃｈｎｉｑｕｅｉｓｐｒｏｐｏｓｅｄｉｎｔｈｉｓｐａｐｅｒ，ｉｎｗｈｉｃｈｓｍａｌｌｓｉｎｇｕｌａｒｖａｌｕｅｓｅｘｃｅｐｔｔｈｅｚｅｒｏｓａｒｅｍａｎｉｐｕｌａｔｅｄｂｙｔｈｅｒｅｇｕｌａｒｉｚａｔｉｏｎ．Ｌｅｓｓｓｏｌｕｔｉｏｎｅｒｒｏｒａｎｄｒａｐｉｄｅｒｓｏｌｖｉｎｇｐｒｏｃｅｄｕｒｅｃａｎａｃｈｉｅｖｅｄｂｙｕｓｉｎｇｔｈｉｓｔｅｃｈ－ｎｉｑｕｅ．Ｎｕｍｅｒｉｃａｌｅｘｐｅｒｉｍｅｎｔｓｓｈｏｗｔｈａｔｔｈｅｐｒｏｐｏｓｅｄｍｅｔｈｏｄｃａｎｐｒｏｖｉｄｅｉｍａｇｅｓｓｕｐｅｒｉｏｒｔｏｔｈｏｓｅｒｅｃｏｎｓｔｒｕｃｔｅｄｂｙｔｈｅｌｉｎｅａｒｂａｃｋｐｒｏｊｅｃｔｉｏｎ（ＬＢＰ），ＬａｎｄｗｅｂｅｒｉｔｅｒａｔｉｖｅｍｅｔｈｏｄａｎｄｔｈｅｓｔａｎｄａｒｄＴｉｋｈｏｎｏｖｒｅｇｕｌａｒｉｚａｔｉｏｎｍｅｔｈｏｄ． Ｋｅｙｗｏｒｄｓ：ｅｌｅｃｔｒｉｃａｌｃａｐａｃｉｔａｎｃｅｔｏｍｏｇｒａｐｈｙ，ｉｍａｇｅｒｅｃｏｎｓｔｒｕｃｔｉｏｎａｌｇｏｒｉｔｈｍ，Ｔｉｋｈｏｎｏｖｒｅｇｕｌａｒｉｚａｔｉｏｎ，ｓｉｎｇｕｌａｒｓｙｓｔｅｍｔｈｅｏｒｙ 文章编号：１００８－０５７０（２００７）０６－３－０２７２－０３ 引言 电容层析成像（ＥｌｅｃｔｒｉｃａｌＣａｐａｃｉｔａｎｃｅＴｏｍｏｇｒａｐｈｙ，简称 ＥＣＴ）中图像重建算法的研究是ＥＣＴ技术和应用的重点环节。 实现图像重建的基本思路是在分析电极激励的静电场问题得到敏感场数据以后，建立被测介电常数与测量电容值之间的关系方程，再运用合适的方法反演截面图像，并要求一定的成像质量和速度。ＥＣＴ图像重建属于逆问题，通常观测数据值远远少于被测数据，而且由于敏感场矩阵本身存在的大条件数，导致求解问题的不适定性，另由于ＥＣＴ系统固有的“软场”性质，待解问题的非线性，使这类问题的求解有一定的困难。 当前存在的ＥＣＴ图像重建算法中，常用的方法有线性反投影算法（ＬＢＰ）、Ｌａｎｄｗｅｂｅｒ迭代法及Ｔｉｋｈｏｎｏｖ正则化算法。ＬＢＰ将问题看成简单线性问题，求解速度快，但是误差较大；而 Ｌａｎｄｗｅｂｅｒ迭代法利用ＬＢＰ得到初始图像，然后计算电容值和 测量电容值之间的误差，反复进行修正，可以得到比较精确的图像，但同时速度慢，不利于实时应用。 Ｔｉｋｈｏｎｏｖ正则化方法用于ＥＣＴ的图像重建，它方法上是引 入一正则化参数试图减小敏感场矩阵的条件数，然后进行求解，但是实质上却对敏感场矩阵所有的奇异值都加上了一个正则化参数，这样对大奇异值项来说，肯定会造成一定的误差，所以本文为了避免这种误差的存在，将解展开为奇异向量的线性组合，通过分析小奇异值对应的项对计算结果产生的影响，给出了选择合适的正则化参数的方法，可使图像重建达到比较理想的结果。１基本理论知识 设存在一病态线性方程组 Ａｘ＝ｙ（１） 式中Ａ属于ｍ×ｎ矩阵，ｘ为ｎ维向量，ｙ为ｍ维向量。标准 Ｔｉｋｈｏｎｏｖ正则化方法将问题转化为求下列的范函最小值问题： （２） 式中λ为正则化参数，该范函极值问题的正则化解ｘα也是 下列方程的唯一解： （３） 设Ａ的奇异系统为 ，即满足 （４） 则可得到 （５） Ａ的奇异值为μｉ，则ＡＴ Ａ的特征值为μｉ２。那么根据特征值 理论满足 （６） 从而使得 （７） 代入（５）式则可得到正则化解为 （８） 可以看出，方程（１）的解可以看作是奇异向量ｘｉ和系数 的线性组合。但是如果系数矩阵Ａ的性态不好，存 在相对很小的奇异值，则相对小奇异值的某些组合分量即具有 很大的系数。此时如果已知向量ｙ存在误差或噪声，并且该噪 曹琳琳：硕士研究生 ２７２－－
No.13,2010 现代商贸工业 Modern Bus iness Trade Industry2010年第13期 地震波层析成像反演方法及其研究综述 冯 微 (长江大学物理科学与技术学院,湖北荆州434025) 摘 要:通过研究利用初至波走时的层析反演方法建立近地表速度模型,提供近地表地下介质的速度信息,进一步为静校正或浅层工程勘探服务。 关键词:速度建模;层析成像;初至波 中图分类号:TB 文献标识码:A 文章编号:1672 3198(2010)13 0368 01 地震勘探是利用人工在地表激发和接收地震波,再对地震波作分析处理以及解释而得到地下构造信息和岩性信息的一种方法。在整个地震勘探过程中,精确的求取地震波在地下介质中的传播速度,一直是地震勘探的核心问题之一。尤其在地表条件较复杂的区域,地表速度的横向剧烈变化会严重影响中深层目的层的成像效果。近地表速度不准确,将会直接影响到速度分析、偏移成像的质量以及静校正的精度等地震勘探的各个环节和最终的勘探成果。 1 地震面波及波形反演 利用面波进行结构反演一直是了解地球介质结构的重要途径。近几年来,在面波理论和面波反演方面做了大量工作。陈蔚天和陈晓非(2001)提出了一种求解水平层状海洋-地球模型中面波振型问题的新算法,它简洁、高效,彻底消除了高频情况下数值计算的精度失真问题。张碧星等(2000,2002)对瑞利波勘探中 之字形频散曲线形成的物理机理和多模性问题进行了理论分析,研究了诸波模的传播特性及相互关系,以及地表下低速层介质的位置、厚度及其它参数对 之字形频散曲线的相互影响.在面波反演理论方面,朱良保等(2001)通过保角变换,把面波群速度的反演变成了球谐系数的线性化反演,使其计算速度快,等值线光滑,构造界限清晰。众多研究者根据从面波资料求出的频散曲线,对不同地区的地下速度结构作了反演,揭示了横向结构差异的广泛存在。 根据走时反演地下结构是获取结构信息的经典做法。刘伊克等(2001)根据三维地震观测的初至走时数据,利用最小平方与QR分解相结合的算法,在三维空间重建近地表低降速带速度模型。同时,采用分形算法克服了初至波波形差异以及折射波相位反转导致的拾取误差,实现了三维初至拾取的大规模全自动化运算。李录明等(2000)针对地震勘探中的复杂地表问题,提出了一套地震初至波表层模型层析反演方法.它利用地震直达波、回折波、折射波以及三者组合的初至波和层析反演方法具有的纵、横向变速优势,实现适应速度任意变化的复杂表层模型反演。 在利用远震体波接收函数反演地下结构方面。钱辉等(2001)对接收函数反演地壳结构速度的算法作了分析,使之适应正演参数的变化,并利用天然地震接收函数揭示了青藏高原东部地壳结构。 近年来,非线性反演越来越受到重视,许多研究者把新的最优化理论引入地震学反演中。孟洪鹰和刘贵忠(1999)提出了多尺度地震波形反演的小波变换方法。对于一维非线性地震波形反演问题,此方法和已有的简单迭代法及多重网格法比较表明,此方法更为有效。杨峰和聂在平(2000)提出了用于二维轴对称非均匀介质结构的反演和成像的一种新的反演迭代方法变分玻恩迭代方法.与传统的玻恩迭代方法相比,其收敛速度和成像质量均有较大改善。 2 地震勘探、测井问题中的地震波研究及其它 在地震勘探和测井方面,许多研究者针对实际问题,提出了新的方法。沈建国和张海澜(2000)计算了井内靠近井壁的偏心声源激发的声场,得到了在井壁不同位置的接收波形,分析了直达波、井壁反射波、纵波、横波和面波在这些波形中的反映。为了处理横向强变速介质中的深度成像问题,程玖兵等(2001)提出一种基于共炮道集的优化系数的傍轴近似方程叠前深度偏移算子,在基于反射系数估算的成像条件下,可实现叠前深度偏移成像。陈生昌等(2001)实现了一种基于拟线性Born近似的叠张海明等:地震波研究前深度偏移方法,扩大了拟线性Born近似的应用范围,使其能够适应更强的横向速度变化。张美根和王妙月(2001)利用有限元法和最小走时射线追踪的界面点法,实现了各向异性弹性波的叠前逆时偏移.陈志德等(2002)利用叠前深度域地震成像对速度模型变化的敏感性,采用偏移迭代逐次逼近最佳成像速度,研究开发了一套快捷有效的三维叠前深度偏移深度域速度模型建立技术。顾汉明等(2002)在频率-波数域中采用解析法,解出多层条件下海底实测的多分量地震数据分解成上行和下行P波和S波的算法,导出海底各层地震反射系数随入射角变化(简称RVA)的递推计算公式。金胜汶等(2002)给出了一种高效率、高精度的炮检距域叠前深度偏移方法,并得到各个不同照射角下的成像结果。 3 讨论和结论 地震波理论是固体地球物理学研究的重要基础.地震波研究领域的任何实质性进展都会促进固体地球物理学的发展.在过去的4年里,中国地球物理学家在该领域做了很多有意义的研究工作,其中不乏创新性的理论工作.当前地震波研究领域的重要课题包括: (1)复杂地球介质中地震波激发与传播理论; (2)高效计算三维介质中地震波传播的数值方法; (3)利用先进的地震波数值模拟方法,开展设定地震与强地面运动的数值模拟研究,为精细的地震危险分析与预测奠定基础。 参考文献 [1]周庆凡.我国天然气发展前景广阔[J].中国石化,2009. [2]刘英祥.我国天然气价格与天然气发展问题研究[J].企业经济, 2009. [3]牛建娣.我国天然气市场供需状况及发展对策分析[D].对外经济 贸易大学,2007. ! 368 !
层析成像 姓名：李文忠 学号：200805060102 班级：勘查技术与工程（一）班
前言 层析成象是在物体外部发射物理信号，接收穿过物体且携带物体内部信息，利用计算机图象重建方法，重现物体内部一维或三维清晰图象。层析成象技术最大的特点是在不损坏物体的条件下，探知物体内部结构的几何形态与物理参数(如密度等)的分布。层析成象与空间技术、遗传工程、新粒子发现等同列为70年代国际上重大科技进展。层析成像应用非常广泛，如医学层析的核磁共振成像技术、工业方面的无损探伤、在军事工业中，层析成象用于对炮弹、火炮等做质量检查、在石油开发中被用于岩心分析和油管损伤检测等，层析成象是在物体外部发射物理信号，接收穿过物体且携带物体内部信息，利用计算机图象重建方法，重现物体内部一维或三维清晰图象。声波层析成像技术 声波层析成像方法所研究的主要内容，一个是正演问题，即射线的追踪问题，是根据已知速度模型求波的初至时间的问题；另一个问题就是反演问题，即根据波的初至时间反求介质内部速度或者慢度分布的问题。层析成像效果的好坏与解正演问题的正演算法和解反演问题的反演算法都有直接的关系。论文详细研究声波层析成像的射线追踪算法，重点探讨了基于Dijkstra算法的Moser曲射线追踪算法，并用均匀介质模型、空洞模型、低速斜断层等模型使用Moser曲射线追踪时的计算精度与计算效率，发现了内插节点是影响Moser曲射线追踪效果的主要因素，得到了内插节点数为5～7之间，计算速度较快，计算精度较高。模型试算的结果表明，正演采用内插10个节点，
反演过程中采用内插5个节点，效果最佳。在层析成像正演算法的基础上，详细研究了误差反投影算法(BPT)、代数重建法(ART)、联合迭代法(SIRT)；研究了非线性问题线性化迭代的最速下降法、共轭梯度法(CG)；重点推导和建立了层析成像的高斯—牛顿反演法(GN)；详细研究了非线性最优化的蒙特卡洛法(MC)、模拟退火法(SA)、遗传算法(GA)；研究了将非线性全局最优化和线性局部最优化方法相结合的混合优化方法，探讨了基于高斯牛顿和模拟退火相结合(GN-SA)混合优化算法。在此基础上，以速度差为10％的低速斜断层模型为例，详细探讨了线性化算法SIRT、GN；非线性最优化算法SA、GA以及混合优化算法GN-SA五种算法对该模型的计算结果，并探讨了直射线和Moser曲射线追踪的反演效果。数值试验表明，基于Moser曲射线追踪的高斯—牛顿反演法的层析成像效果最佳，计算效率最高。采用基于Moser曲射线追踪的高斯—牛顿法，对速度差为25％的等轴状空洞构造、速度差为33％的不连通空洞模型、速度差为33％的高速岩脉进行了反演试算，对于这些理论模型，高斯—牛顿法均取得了较好的成像效果。为进一步验证各种层析成像法，在实验室制作了水泥台和石膏板实物模型，并分别在水泥台中央制作一个方形空洞，在石膏板中央制作一个倒“L”形空洞。对这两个实物模型进行了实测，对测量的数据，用高斯—牛顿法进行层析成像反演，均取得了较好的成像效果。通过本文的研究和数值试验，得到了以下结论：(1)基于直射线追踪方法，适用较为简单的地质体，亦或是测量精度要求不高的问题。由于直射线追踪方法在成像过程中，只需要追踪一次就可以
2001年9月地球物理学进展第16卷第3期井间地震层析成像的现状与进展 裴正林 （石油大学（北京）物探重点实验室，北京，100083） 摘要：综述了井间地震层析成像研究的现状，给出了小波变换域井间地震层析成像方法的最新 进展，并对井间地震层析成像研究给予展望. 关键词：井间地震层析成像；小波多尺度；研究进展 中图分类号：P315.3+1文献标识码：A文章编号：1004-2903（2001）03-0091-07 1井间地震层析成像的研究现状 井间地震层析成像也称为井间地震CT技术，它能够提供被探测地质体的构造和岩性 分布的高分率图像.井间地震CT技术是从医学CT技术发展起来的，其数学基础是Radon变换.井间地震CT的研究基本始于20世纪70年代初，80年代处于对大量模型数据和少量实 际数据的成像研究阶段，90年代以来，井间地震CT进入实用化阶段，并取得不少可喜成果，同时，也逐渐意识到射线CT所固有的缺点，开始研究波动方程CT. 从地震波的运动学和动力学特征出发，井间地震CT方法可分为两大类：一类是基于几 何光学或射线方程的方法称之为射线CT；另一类是基于波动方程的方法称之为波形CT.当 非均匀体的线性尺度大于地震波长时，射线CT是适用的；而当非均匀体的线性尺度与波长 相近时，衍射和散射就起主导作用了，基于射线理论的成像方法就不再适用，这时候必须用 波动方程CT方法. 井间地震层析成像方法主要包括两部分：正演方法和反演方法.井间地震层析成像的正 演方法可分为两种；一是射线追踪方法；二是波场的数值模拟方法. 射线理论和射线方法是研究地震波传播理论的重要方面之一.用射线理论可以研究地 下复杂构造、横向不均匀介质中的地震波传播问题.经过射线追踪，计算地震波的走时、波前 和射线路径. 70年代以前的各种射线追踪方法一般适合于较为简单模型的射线追踪［1］.由于实际 的介质速度变化较大（速度差大于10%），因此，需要研究复杂结构模型的射线追踪方法. 收稿日期：2001-03-15；修订日期：2001-06-15. 基金来源：“九五”国家科技攻关项目资助（959130602）. 作者简介：裴正林，1962年生，2000年获中国地质大学（北京）地球探测与信息技术专业博士.高级工程师，现在石油大学（北京）从事博士后研究.主要研究方向：信号处理，小波变换、遗传算法及神经网络应用，层析成像理论方法和地震数据 处理、偏移方法等方面研究.E-mail：zhenglinpei@https://www.doczj.com/doc/3f3669419.html,.
地震层析成像——（一）模型参数化 冷独行整理 地震层析成像（seismic tomography）是指利用大量地震观测数据反演研究区域三维结构的一种方法。其原理类似于医学上的CT，但地震层析成像比医学上的CT技术更复杂。大量数据以及其他许多不定因素，包括存在多种数据误差、解的不唯一性在内的地球内部成像问题。 Aki和Lee[3]以及Aki等[4]利用区域台阵的三维成像，以及Dziewonski等[5]对全球大尺度上地幔速度结构的勾画成为成像研究中开拓性的工作 地震层析成像是典型的地球物理反演问题，大多数地震层析成像问题都涉及到以下几个方面：①模型参数化， ②正演（射线追踪）， ③反演， ④解的评价。 一、模型参数化 成像的目的就是要获得接近实际地下结构的模型，所以在成像前必需要建立模型来描述地层结构，而且选取模型的好坏决定了获得地层结构信息能力的好坏。 过于简化的模型可能使结构中有意义的信息被忽略，复杂的模型可能使反演的不确定性增强，同时可能引入虚假信息。。 模型参数化可分为两类。 一类是Tarantola和Nercessian等提出了“不分块”的参数化。不对模型进行离散化，反演完全在泛函空间中进行，只是在最后计算想要的截面时采取离散化。由于反演在泛函空间中进行，理论上可以计算空间任何位置上的速度，结果不受离散化的影响，有利于成像的显示。 另一类是离散化的模型参数化。其优点是数学上容易处理，运算相对简单；缺点是在一般方法中出现的某些简化，在用离散时可能被掩盖掉。现在通用的大都是离散化的模型参数化，通常采用两种方法来表示地层结构。 一种是使用少量参数确定三维解析函数（如，Dziewonski；Spencer和Gubbins）， 例如：Woodhouse、Dziewonski[19]和Su等[20]在全球地震层析成像使用球谐展开来表示模型；Burmakov等将速度扰动展开成一定阶数的切比雪夫多项式，以减少未知量个数，提高求解效率；朱露培提出的频谱参数化法，将待求扰动场按其空间频率展开，反演各阶频率系数。解析函数法优点是模型参数少，最终的模型被压缩且易于被其他研究者利用，还有一些计算方面的优势。其缺点是不易对小尺度局部异常成像、模型边界不稳定，而且观测方程中的满系数矩阵是反演所不希望的。 另一种是采用离散函数，根据赵大鹏教授的综述，使用在区域地震层析成像中，如分块法或节点法，这种方法可模拟局部小异常，产生大型稀疏系数矩阵，利于在现代计算机上计算，所以目前几乎所有局部和区域的体波层析研究都采用分块或节点法。分块法将地球按某种准则划分为有限个三维块体，假定每块内速度是常值。 分块法在进行射线追踪和走时计算时都十分简便，但是在模型中人为的引入了块体间的
毕业设计开题报告 电子信息工程 电容层析成像技术测量电路的设计 1、选题的背景、意义 过程成像(PT:PROCESS TOMOGRAPHY)技术是近年来才发展起来的一种两相或多相流测量技术,其优点是利用被测物体外部的检测信息,获得被测物体内部变化∕高速流状态。过程成像经常使用特殊方法设计的探测器,通过非侵入式的方法取得被测两相流或多相流介质的场(如电磁场)信息,可以根据场的信息和被测物体的作用原理,应用数学的方法重建两相流或多相流在管道内或反应装置的内部的横截面上的动态分布的情况。在我们日常生活中，过程成像可用于研究化工、石油等各种固体、气体的物料输送管道中的气或固两相流和气或固或液多相流得流态化、反应、扩散以及混合等动态过程,以监控反应器中气泡的分布和大小以及反应器中气泡的破碎和合并等过程；通过工业过程中的建立的模型,研究反应器中反应速率、质量传递以及热量传递的关系,提高反应器的选择性、转化率以及安全性等[1]。 电容层析成像技术（ECT）是医学CT技术在工业流动过程上的改革与发展，是目前用来解决多相流参数测量难度大的最新手段。ECT(Electrical Capacitance Tomography)是在应用于多相流参数检测的一种新型技术,原理是依靠检测非导电物场内介质分布变化引起的电容值的变化,通过某种图像重建算法来反演物场内的介质分布,从而实现对两相流参数的测量。工业过程成像技术中，电容的成像技术（ELECTRICAL CAPACITANCE TOMOGRAPHY,简称ECT）以它廉价、高速和非辐射等特点，在近十几年来获得很大发展[2]。 其实，早在二十世纪八十年代中期，以英国曼彻斯特理工大学BECK M S教授为首的研究小组就已经提出了“流动成像”（FLOW IMAGING）得概念，并研制成功了8电极的电容成像系统。在国外，美国能源部MORGANTOWN研究中心几乎与BECK的研究小组同时发明出了一种在线监测流化床中空隙率分布的16电极电容的成像系统（CAPACITANCE IMAGING SYSTEM,简称CIT），该系统可用于对流化床内物料密度三维分布地监测。电容成像的技术应用于工业上的多种需要进行多
2011-8-3 基金项目：国家自然科学基金(60762001)；广西高等学校优秀人才计划（桂教人才0804）。 作者简介：赵进创(1968-)，男，博士，教授，研究方向：电成像技术；嵌入式系统开发等；刘金花(1987-)，女，硕士研究生，研究方向：电成像算 法。 收稿日期： 修回日期： 改进敏感场的电容层析成像图像重建算法 赵进创,刘金花,黎志刚,傅文利,李贤宇 ZHAO Jin-chuang,LIU Jin-hua,LI Zhi-gang,FU Wen-li,LI Xian-yu 广西大学计算机与电子信息学院, 广西 南宁530004 College of computer , electronics and information, Guangxi University, Nanning 530004, China E-mail: zhaojch@https://www.doczj.com/doc/3f3669419.html, Image reconstruction algorithm based on updated sensitivity field for ECT Abstract: The Landweber image reconstruction algorithm based on imaging sensitive field mean filtering method is proposed to solve the so-called “soft -field” characteristic problem of sensitive field of Electrical Capacitance Tomography (ECT) system. The algorithm principle is that the neighborhood pixel sensitivity is averaged by template convoluting method, which can reduce the sensitivity of region near electrodes and improve that of central region of pipe. The algorithm can eliminate partly the affect on the quality of image reconstruction due to uneven sensitivity and improve image reconstruction accuracy. Simulation results indicate that the algorithm is superior to conventional Landweber algorithm in image reconstruction accuracy and convergence speed. Key word: ECT; Landweber image reconstruction algorithm; Sensitivity; mean filtering 摘 要: 针对电容层析成像系统中敏感场的“软场”特性，提出了一种基于成像敏感场灵敏度均值滤波的Landweber 图像重建算法。该算法是通过模板卷积的方式对敏感场灵敏度进行邻域平均，降低靠近极板区域的灵敏度，提高管道中心区域的灵敏度，部分消除了因敏感场不均匀对图像重建质量的影响，提高图像重建精度。仿真结果表明，该算法在图像重建精度和收敛速度上均优于传统的Landweber 图像重建算法。 关键词：电容层析成像;Landweber 图像重建算法;灵敏度;均值滤波 DOI: 文章编号： 文献标识码: A 中图分类号：TP212.9 1 引言 电容层析成像 ( Electrical Capacitance Tomography , ECT) 作为一种非侵入式的流动参数前景。其原理是通过计算机采集安装在封闭的工业管道、容器外壁的传感器阵列在不同观测角度下的投影数据即电容测量值，采用相应的图像重建算法 显示被测物场的二维或三维介质分布图像[6][11][12]。 图像重建算法是ECT 系统的关键技术之一，目前国内外研究ECT 图像重建常用的算法主要分为2 类: 一类是非迭代算法，如线性反投影算法（LBP ），另一类是迭代算法，如Landweber 迭代法[1]-[5][9]。LBP 算法简单、成像较快, 但其重建图像精度低。 Landweber 迭代法是利用LBP 法重建的图像作为迭代过程的初值, 由于初值有时偏离实际值较大, 造成迭代误差累积，影响图像重建质量和算法收敛速度。本文针对此问题提出一种基于灵敏度矩阵均值滤波的Landweber 迭代算法，与传统的Landweber 迭代法相比，该算法成像质量高，收敛速度快。 2. ECT 图像重建模型 ECT 系统的正问题就是由已知的介电常数分布,求出传感器各极板对之间形成的电容值,可表达为如下的数学模型[6][10]-[12]： dxdy y x y x S y x C D ij ij )),(),,((),( （1） 网络出版时间：2011-10-24 10:08 网络出版地址：https://www.doczj.com/doc/3f3669419.html,/kcms/detail/11.2127.TP.20111024.1008.007.html
(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201910127240.0 (22)申请日 2019.02.20 (71)申请人 天津大学 地址 300072 天津市南开区卫津路92号 (72)发明人 董峰 梁光辉 任尚杰 (74)专利代理机构 天津市北洋有限责任专利代 理事务所 12201 代理人 程毓英 (51)Int.Cl. G01N 29/44(2006.01) G01N 27/00(2006.01) (54)发明名称 基于统计逆的电学/超声双模态内含物边界 重建方法 (57)摘要 本发明涉及一种基于统计逆的电学/超声双 模态内含物边界重建方法，该方法采用径向形状 模型对待重建内含物边界即目标边界进行参数 化表征，然后利用电学成像模型构建形状系数估 计的似然模型，利用超声反射成像模型构建形状 系数估计的先验模型，再通过最大后验估计法求 解最佳形状系数， 实现内含物边界的重建。权利要求书1页 说明书6页 附图2页CN 109946388 A 2019.06.28 C N 109946388 A
1.一种基于统计逆的电学/超声双模态内含物边界重建方法，该方法采用径向形状模型对待重建内含物边界即目标边界进行参数化表征，然后利用电学成像模型构建形状系数估计的似然模型，利用超声反射成像模型构建形状系数估计的先验模型，再通过最大后验估计法求解最佳形状系数，实现内含物边界的重建。步骤如下为： (1)利用径向形状参数化模型，在局部极坐标系下将目标边界剖分成一组等极角分布的离散点，记离散点到局部极坐标中心的距离为径向距离，并用r i ,i＝1,2,…,N来表示，目标边界用一组形状系数r＝[r 1,r 2,…,r N ]来表征； (2)基于贝叶斯统计理论，当给定电学测量电压V m 和测量噪声e，形状系数r的后验估计表达如下： π(r|V m )∝πe (V m -V(r))· π(r)式中，π表示概率密度，∝表示正比号，V m 表示测量电压，V(r)表示由电学成像正问题模型计算获得的边界电压，π(r|V m )表示形状系数r的后验概率密度，π(r)表示形状系数r的先验概率密度，πe (V m -V(r))表示形状系数估计的似然模型， 满足如下公式：式中，Γe 表示测量噪声e的协方差，表示Γe 的逆； (3)对于步骤(2)中形状系数r的先验概率密度π(r)，通过超声反射成像技术确定：首先，利用超声反射成像技术获得目标边界上部分离散点位置信息，然后采用等极角插值技术计算目标边界形状系数的初始估计r u ，最后， 构建形状系数的先验概率模型如下：式中，表示形状系数协方差矩阵，表示的逆； (4)结合(2)和(3)，当给定电学测量电压V m 和测量噪声e，形状系数r的后验概率密度表 示成如下形式： (5)对于(4)中形状系数r的后验概率密度，采用最大后验估计法进行求解，形状系数r 的最佳估计值通过求解如下公式获得： (6)由(5)获得的最优形状系数能够计算目标边界上离散点坐标，再采用光滑函数对目标边界上离散点进行拟合来获得目标边界。 权 利 要 求 书1/1页2CN 109946388 A
第84卷 第6期 2010年6月 地 质 学 报 AC TA GEOLOGICA SINICA V ol . 84 N o .6 June 2010 注:本文为国家自然科学基金项目(编号40774051,40404011)、中国地质科学院地质研究所基本科研业务费(编号J0707,J0803)、科技部中美国际合作项目(编号2006DFA21340)和国家专项“深部探测技术与实验研究”(编号SinoProbe -02)资助的成果。收稿日期:2009-07-05;改回日期:2010-01-20;责任编辑:章雨旭。 作者简介:贺日政,男,1973年生。博士,副研究员,主要从事利用天然地震波探测青藏高原深部结构与构造研究。Em ail :herizheng @cags .ac .cn 。 地震波速层析成像方法研究进展 贺日政,高锐,郑洪伟,管烨,李秋生,李文辉,熊小松,邓攻 中国地质科学院地质研究所,北京,100037 内容提要:本文回顾了利用天然地震观测获取地下速度结构的方法。尽管有不同的新方法涌现,天然地震波速层析成像方法,尤其是多震相联合反演的格点层析成像方法,是当今使用广泛使用的层析成像方法之一,是对地球内部成像的最有效方式。波速层析成像方法的未来发展首先是提高第一手的观测资料,即增加接收地震波信息的地震台站分布密度;同时,通过多种地球物理方法联合反演相互约束可以给出较为严格的地球物理模型,并来降低了地球物理反演和解译的多解性,这是地球物理探测研究的趋势,也是天然地震波速层析成像方法的研究趋势。 关键词:地震波;速度;层析成像方法;研究进展 经过近30多年的快速发展,地震学已成为研究地球内部结构的主要手段,是深部地球物理探测技术中首选技术。天然地震波的非凡穿透能力,同时地震波包含着其传播过程中所穿越地球内部结构的丰富信息,使得天然地震学研究是当今地球内部结构的主要方法技术之一。地震波速度成像技术常见有三种,即波形拟合反演,接收函数方法,地震波层析成像。 波形拟合:自Woo dhouse and Dziew o nski (1984)首先利用波形拟合方法分析了全球数字化台站数据后,波形拟合方法得到了广泛运用(Chen ,1993;Song and H elmbe rg er ,1992,1998;Minkoff and Symes ,1997;Pratt ,1999;黄建平等,2009)。目前波形拟合反演技术充分利用从震源至台站间的全波形信息,既可以对震源结构也可以对接收台站区域反演,或二者同时反演获取其目标区域的速度结构特征,甚至还可以模拟地震破裂过程等。波形拟合方式可以直接分析地震波在传播过程中受介质的影响,且直观地给出拟合波形与实际波形记录的对比结果(Aki and Richards ,1980)。尽管最近十年计算机技术的飞速发展,为波形拟合广泛运用提供了基础,但由于计算量非常大,利用波形拟合反演获取区域性的三维速度结构仍不是首选。 接收函数:自Vinnik 于1977年介绍用P -SV 转换波接收函数方法研究地幔结构(Vinnik ,1977) 以来,利用接收函数方法获取接收台站下方的速度结构信息现已成为天然地震学研究中又一手段 (Am mon et al .,1989,1990;刘启元等,1997),特别适合于对台站下方界面的研究。目前这种方法在国内已经普遍运用于小区域布置密集台阵剖面研究当中(刘启元等,2000;吴庆举等,2004;Chen and Ai ,2008;Xu and Zhao ,2009)。但由于台站分布、多次波影响和方法本身限制,接收函数方法获取的速度结构只是台阵下方局部二维结构,除非台站分布密度较大,否则还不大适合于三维结构反演。 地震波层析成像:与上述两种不同的是,自地震波层析成像技术(Aki and Lee ,1976)出现以来,地震波层析成像技术很快成为获取地壳/上地幔速度结构的最有力的技术手段。地震波在传播过程中受到地球内部物性的影响(Shea rer ,1999),记录到的地震波包含了其所穿越地球内部区域的速度结构等信息,据此可以获取大尺度范围内的地球非均匀速度结构,进而研究地球地幔内部物性特征。因此,地震波层析成像是当今研究地球内部基本圈层三维结构最有利的技术手段之一。 1 层析成像方法研究进展 层析成像技术首先由Aki 等提出,并给出了小尺度(Aki and Lee ,1976)和区域尺度(Aki e t al .,1977)远震体波层析成像(Teleseismic Body -w ave
论电学原理在日常生活中的应用 摘要：电学原理的出现，在人类历史上具有非凡的意义，它不仅促进了科学技术的进步，而且在人类日常生活中也得到了广泛的应用。随着电学知识的普及和发展，人们对于电学的掌握和了解也越来越深，通过将其应用到日常生活中，提高了人们生活水平。本文笔者将论述电学原理在日常生活中应用的事例和影响，浅析电学原理，也希望电学原理能够在更多方面造福人类。 关键词：电学原理；日常生活；应用 前言： 随着电在人们生活中的所扮演的角色越来越重要，如何将电学原理更好的融入我们的生活显得格外重要，笔者将在本文中简介一些电学知识，通过列举和介绍生活中的电学应用的实例，突显电力原理在人们日常生活中的重要性。 一、电学内容及其发展历程 1.1电学知识 电学主要研究的内容包括静电、静磁、电磁场效应、电路等。下面笔者将会简单介绍一些电学的有关知识，以便读者更好的理解电学。 首先笔者介绍的是基础电路知识，电路就是电流所流经
的路径，是由电气设备和元器件，按一定的方式联接起来。电路由电源、负载、连接导线和辅助设备组成，而电路的作用是进行电能与其它形式的能量之间的相互转换。其次，笔者将论述一个定律：欧姆定律。在同一电路中，导体中的电流跟导体两端的电压成正比，跟导体的电阻成反比，这就是欧姆定律。最后，笔者带领大家走进电与磁。在历史上，电和磁的研究是分开的。当磁被发现和研究时，人们从中发现了磁场，而从安培定则和丹麦的奥斯特中我们可以得知电是可以产生磁场的，所以后来人们又在两者的研究中发现磁也可以产生电流。说到这，笔者不得不提到电动机，电动机就是利用这一原理，它利用通电线圈产生磁场，然后磁场对电流的受力，使电动机转动。 1.2电学原理发展历程 电学是物理学的一个分支，其发展得到了古今众多物理学家的不断完善。从公元前的琥珀和磁石，到磁和静电、雷和静电，电学其实已经有了概念，但是其真正的开始，始于18、19世纪的欧洲。到后来，伏打电池的利用与电磁学的发展，使得电磁学的发展突飞猛进，到了20世纪，随着原子物理学、原子核物理学、粒子物理学的进步，人们的视野也转入了微观领域，这时候经典电磁学受到了挑战，出现了很多经典电磁学不能解释的问题，因为经典电磁学忽略了粒子性方面，所以在一些物理学家的推动下，量子电磁理论应运
地震层析成像 摘要：层析成像方法是一种公认的基于地震数据的有效方法，近20年来，层析成像方法发展迅速。从原理上讲，层析成像方法可分为两大类，一是基于射线理论走时层析成像，二是基于波动方程的散射层析成像。本文介绍新的层析成像方法及其技术，包括各向异性介质的2D立体层析成像；时移层析成像的超声数据试验；绕射层析成像的迭代方法：真振幅偏移的本质；用于速度模型构建的下行波折封层析成像和反射层析成像；多尺度波动方程反射层析成像，并在后面展开层析成像方法应用于构造速度模型的分析和实例。 关键字：层析成像；偏移成像；速度模型；克希霍夫偏移。 一、引言 偏移成像在地震勘探和开发过程中，已经成为一种关键的地震数据处理技术。成像的精度和可靠性依赖于速度模型的准确与否。 速度分析历来都是地震资料处理的基础工作，从均方根速度、层速度以及叠加速度等，贯穿于地震资料处理的方方面面，速度分析方法丰富多样。迄今，层析成像方法是一种公认的基于地震数据的有效方法，近20年来，层析成像方法发展迅速。从原理上讲，层析成像方法可分为两大类，一是基于射线理论走时层析成像，二是基于波动方程的散射层析成像。后一种层析成像很复杂，正处于理论研究阶段。尽管其实际应用不多，但却是层析成像的发展方向。 走时层析成像比较成熟，有很多的实际应用。它又可细分为初至走时层析成像和反射走时层析成像。初至走时层析成像方法简单直观，稳定性较好，主要应用于井间地震以及近地表的速度分析，但是，初至走时层析成像由于只利用初至走时，所以，得到的速度模型比较粗糙，分辨率也较低。 反射层析成像主要应用于地下速度和反射层深度的反演，以及叠前或叠后偏移的速度分析之中。前者由于速度和深度之间的藕合关系，以及反射波到达时间及其层位难于拾取等，制约了它的广泛应用，但是，这是一种极具价值和潜力的反演方法。后者则是利用经过叠前或叠后CRI道集中同相轴未被拉平的剩余时差，经过层析成像来修正用于偏移的速度模型。这种构建速度模型的方法，目前正广泛应用于叠前深度或时间偏移中。 值得关注的还有，地震资料与其他地球物理资料间的联合反演，其反演结果互为验证、相得益彰，为我们提供了更为可靠的反演结果。 二、新的层析成像方法及其技术 1.各向异性介质的2D立体层析成像 立体层析成像是一种利用局部相关同相轴作为输人的斜率层析成像方