共查询到20条相似文献，搜索用时 62 毫秒
1.
由投影重建图像的对称网格迭代算法 总被引：1，自引：0，他引：1
本文对于工业CT检测中常用的代数迭代重建算法提出了改进,利用投影射线之间存在的几何对称结构,提出了图像重建的对称网格迭代算法(简写为SM-IRT).该算法简化了投影系数矩阵的计算,调整了迭代算法逐线校正的迭代顺序.对模拟数据和工业CT实测数据进行了重建图像的数值实验,结果表明:与常规算法比较,本文提出的新算法重建速度快,成像精度高. 相似文献
2.
冻土应力-应变曲线的分形逼近 总被引：6，自引：0，他引：6
基于冻土应力 应变曲线的细观结构具有分形性质的事实 ,提出了逼近冻土应力 应变曲线的分形方法 .利用分形几何中的线性双曲迭代函数系统 (LHIFS)理论 ,选取垂直压缩因子为参数建立了逼近冻土应力 应变曲线的LHIFS ,其次讨论了使LHIFS生成的线性分形插值函数 (LFIF) (即LHIFS的吸引子 )成为冻土应力 应变曲线的最佳逼近时垂直压缩因子应取的最佳点的计算方法 ,然后给出了利用最佳点求冻土应力 应变曲线分形维数的方法 ,最后给出了简单例子 .提出的方法为编程模拟具有分形性质的地质材料的应力 应变曲线及计算其分形维数提供了一般方法 . 相似文献
3.
针对层析成像模型矩形网格剖分存在的一些问题,提出了复杂结构三角网最小走时射线追踪层析成像方法。以Delaunay三角剖分的优化准则,根据模型点、线、面的几何结构关系,进行三角网格剖分;采用三角网波行面扩展法计算声波初至走时,追踪出任意接收点至激发点处的射线路径,射线路径包含了从接收点至源点的坐标、所在三角单元等信息,形成矩阵方程。应用稳定性较好的联合迭代重构技术求解矩阵方程,得到模型的速度分布。数值模拟结果表明,三角网射线层析成像方法分辨率高,成像结果更接近实际结构形态,有效地解决了矩形网剖分对复杂区域网格参数化灵活性差,速度界面描述精度低等问题。 相似文献
4.
5.
在工业CT成像应用中,由于复杂结构件形状、结构、组分等特征的限制,使得射线能量与有效厚度不匹配,常规的扫描方式易使投影数据缺失,投影数据完备性较差。为此,本文通过改变轨迹来弥补投影数据缺失,推导任意轨迹的成像模式。首先,在任意轨迹的基础上进行投影几何关系建模,结合空间几何知识完成投影矩阵的刻画;其次,基于迭代重建算法对轨迹的无约束化,进行任意轨迹成像模式下的迭代重建;最后通过小锥角、大锥角以及厚度差异比较大的构件成像仿真实验,验证算法的正确性和可行性。对于厚度差异比较大的工件,该方法相对于传统的圆轨迹扫描,投影数据更完备,CT重建质量高。 相似文献
6.
本文将剪力墙结构根据精度需要划分成若干6节点矩形单元,每个单元由2根斜置正交弹簧及联结的上下水平刚性横杆组成,2根斜置弹簧方向分别与外力作用下的主拉、压应力方向一致。本文引入单元位移场线性假定及斜压场理论计算弹簧内力,根据变形协调条件由单元结点位移计算弹簧变形,由弹簧刚度及倾角推导了该宏模型单元刚度矩阵,迭代计算使结构内力及位移逼近真值。算例分析表明,该模型能较好地反映剪力墙的受力性能。 相似文献
7.
本文通过柯尔-柯尔模型研究地下不同几何参数的极化椭球体在地表的中梯双频激电响应,从而了解不同产状极化体产生极化率的特征,便于后期交互建模。由理论推导可知,任何具有激发极化特征的地质体在地表产生的激电异常都可近似为一个地下电偶极子的三维空间积分,这个积分形式与磁场正演计算公式一致,这意味着我们可以用解释磁法异常的方法解释激电异常,即磁类比法。只要根据先验信息实时改变模型,使激电正演得到的曲线与实测激电曲线逼近,当拟合差小于均方误差时,我们将最后一次模型近似看作是真实模型,由此获得地下极化体的顶埋深、倾角和几何大小等空间参数。实测数据测试表明,这种磁类比人机交互方法加入矿体的一些已知信息后,其拟合反演结果比自动反演更具合理性。 相似文献
8.
9.
利用X射线计算机断层成像（CT）方法对页岩的扫描成像是一种无损的，对研究页岩微纳孔隙结构有重要意义的方法.传统的CT重构通常使用的是显式的滤波反投影（Filtered Back Projection，FBP）方法，该算法具有较快的成像速度，但常伴随有伪影或不稳定等情况.对于纳米CT而言，可以通过迭代优化的方法对投影数据进行成像，传统的迭代成像有收敛速度慢导致的计算时间长等缺点.Kaczmarz算法作为一种重要的代数重建技术（ART），由于其几何意义明显，操作容易等优点，在CT重构中起着重要的作用，我们可以通过块状迭代或随机迭代的方式对其收敛速度进行改进.对于所求解问题的不适定性，代数重建过程中需要引入正则化的技巧来改善解的稳定性.本文根据实际问题的需要，使用页岩数值模型，验证了正则化Kaczmarz方法的有效性，并对重庆漆辽龙马溪组页岩样品的实际数据进行了处理，得到了较好的效果. 相似文献
10.
三维场地波动传播的快速射线追踪法 总被引：7，自引：0，他引：7
本文基于Snell定律和Fermat原理对三维任意界面情况下的两点间射线追踪问题进行了研究，从同一条射线满足相同的射线参数出发，推得一个适用于任意界面情况下计算反（折）射点的一阶近似公式。结合迭代技术，给出了三维场地条件下射线追踪迭代算法的计算格式，并进行了三维场地射线追踪模拟计算。计算表明：计算速度相当快，且其计算精度可以根据需要满足要求。 相似文献
11.
12.
地震直达波走时层析成像可归结为求解一个大型的、稀疏的、常常是病态的线性方程组.求解方程组常用的迭代法,需要一个比较合理的初始猜测解,也即是初始速度模型.初始模型关系到反演的效率甚至成像的正确性.本文在前人研究基础上提出一种生成模型网格节点初始速度方法,假定震源到检波点路径为直线,记录每条射线穿过的单元和统计每个网格单元穿过的射线数目、自动拾取网格节点所在单元的数目等.实例中,由原始数学模型的正演旅行时资料生成节点初始速度模型,效果可以.最后,分别采用均匀模型和本文方法生成的初始模型进行迭代反演,通过比较,证实该自动生成节点初始模型的可行性和可靠性,并对存在的问题进行讨论和解释. 相似文献
13.
《地球物理学进展》1989,(3)
本文提出了由地震折射波初至走时确定二维速度结构的一种迭代层析成象反演方法。该方法适宜源检距密于常规剖面的折射剖面,反演方法基于线性化问题的迭代解,可用来确定连续变化的速度和亚水平(Subhorizontal)介面的几何形状。每次迭代中,都用试射法,进行两点射线追踪以构成线性方程组。速度场用速度梯度为常数的三角单元的速度来确定,射线路径可解析计算。根据正问题的各别的线性化公式,考虑了两种不同的反演方法。用线性化走时-速度雅可比行列武反演比用类似于普通级数展开法的慢度公式得到的结果要好。有关分辨率的实例表明,由于射线几何特性、分辨率随深度而下降,引起了水平模糊以及源- 相似文献
14.
在研究地壳结构的人工源宽角反射地震资料解释中,常规宽角反射波走时和射线路径计算大都假定地壳模型为层状块状均匀介质.为了逼近实际地壳结构模型,要求模型尺度较大,为了提高地震资料解释的可靠性,须减小模型离散单元的尺寸,但同时计算量大大增加,使资料解释的效率较低.为此,本文尝试同时提高宽角反射地震资料解释效率和可靠性的方法,即使用双重网格计算宽角反射地震波走时和射线路径的最小走时树方法.双重网格法在均匀介质内部仅计算大网格节点,在速度变化点、震源点和检波点区域,同时计算小网格节点;在界面边界点使用比介质内部节点更大的子波传播区域.模型计算结果表明,对于大尺度的层状块状均匀介质模型,在保证精度的条件下,本文所提出的双重网格射线追踪方法的计算效率比单网格方法显著提高. 相似文献
15.
对已知和未知数字图像分别沿射线作Radon变换(线积分),用同一种布置射线特殊方法布置射线,每条射线积分步长易知皆为常数。从而,得到已知数字图像每条射线的线积分近似值(“观测值”);易建立未知数字图像线性代数方程组。应用“代数重建法”加法修正迭代编制Matlab程序,用此程序处理“观测值”数据,重建未知数字图像,其数值计算结果较好,与已知数字图像相对误差不超过2%。本文是ART加法修正迭代的基础工作,可为有关部门提供研究“代数重建法”实际应用参考。 相似文献
16.
最短路径射线追踪方法是计算地震波走时的主要方法之一,该方法基于惠更斯原理和费玛原理,具有稳健、适于复杂介质模型的优点.为处理方便,最短路径方法中的介质模型通常以规则网格进行剖分,界面节点(界面与网格的交点)以其邻近的模型单元节点(即边界单元节点)近似表示.界面近似将导致计算误差,对于反射波尤为严重.反射波的走时精度可通过减小网格的尺寸提高,但这样会大大增加计算时间,为高精度和高效率地计算地震反射波走时,我们提出了一种基于规则网格的走时校正技术.地震波传播至或起始于边界单元节点的走时校正为地震波传播至或起始于该边界单元节点所对应的界面节点的走时.数值模型计算结果表明,走时校正方法可使反射波的走时精度提高约1～2个数量级,而其计算时间则和常规算法基本上在相同量级. 相似文献
17.
18.
提出了位场曲化平的新方法. 给定观测曲面S上的位场、S对下方水平面P的相对高程,确定P上的位场. 利用由P向上延拓到S的积分式,建立这两个面上位场及相对高程三者所满足的方程,它是第一类Fredholm积分方程. 用Fourier逆变换式把这一空间域积分式化为波数域积分式,再由指数函数的Taylor展开进一步化为级数式. 积分方程的解采用逐次逼近法迭代计算,即用S上的位场观测值作为P上位场的初始迭代值,用导出的级数式求得S上的位场计算值、由S上的位场观测值与计算值之差校正P上的位场,多次迭代,直到满足迭代终止准则. 我们还给出该积分方程的波数域迭代计算方法. 模型算例表明,重力异常曲化平的均方差和磁异常曲化平的均方差分别为0.0008 mGal和0.0019 nT,在主频为2.26 GHz的笔记本电脑运行,2048×2048数据量,计算时间是975 s. 野外磁场实际资料处理也证实这种方法的有效性. 相似文献
19.
《地球物理学进展》1990,(3)
本文讨论了射线和射束理论地震图对速度模型和弯曲界面变化的敏感性.Nowack和Lutter(1989)以前的成果给出走时和射线振幅相对于光滑变化的速度变化导数,在实施构造的线性化最大似然反演时这些导数是需要的.在射线近似中,局部应用Snell定律、校正波前曲率以及使用局部平面波反射/透射系数来体现光滑界面,利用费马原理可沿着原始射线轨迹直接计算走时的偏导数.对于反射/透射射线的振幅扰动,必须考虑扰动后的两点射线轨迹的射线偏移.本文遵循的方法是利用没有额外射线追踪的扰动理论近似计算扰动后的两点射线,然后沿着这条近似的两点射线直接计算扰动的射线振幅,包括修正过的反射/透射系数和几何扩展.进行了若干次数值实验,利用走时和振幅反演速度和界面形状来检验所得出的偏导数算子、走时和振幅反演结果,还与对大尺度性质不太敏感而对不均匀曲率较敏感的振幅作了比较. 相似文献
20.
本文采用木板、有机玻璃板制作二维物理模型，并在模型中加入三个大小不一、位置各异、用石膏充填的方洞。对两个物理模型布置相同的观测系统，形成相同的射线方式，在等间距矩形像素条件下，每条射线的积分步长（射线穿过单个像素的射线段长度）皆为常数。由超声仪器沿指定射线采集各条透射波射线的走时数据（即Radon变换线积分的“观测值”）。应用“代数重建法”（ART加法修正迭代法）求解Radon变换的离散化方程组，重建二维物理模型的CT数字图像。在重建的模型CT数字图像中，能够有效地识别出“石膏洞”的存在和位置，验证了“代数重建法”CT成像技术及其对模型结构识别的有效性。本文是ART加法修正迭代的基础性应用工作，可为有关部门提供“代数重建法”CT成像实际应用参考。 相似文献