Radon变换在高维颠簸方程反问题中的应用.pdfVIP

摘要 Radon变换在高维波动方程反问题中的应用 摘要 波动方程声速的反演传统解法是依据散射和逆散射理论，在一维空 间该问题已经有了许多成功的结果，但是在高维的情况许多理论和算法 还有待去研究．和传统方法不同的是我们根据解的奇性演化规律用摄动 方法导出了奇点走时和速度系数的Radon变换之间的关系．从而速度系 数的反演转化为Radon反变换问题．而后者的理论和算法是相当成熟的． fRadon变换是对图象沿直线的积分过程，即对分布在有限区域QI的 图象f(x∽，过点‰帅)并与_y轴夹角为殖线的投影，在复数坐标牙什涉 下表示为 吧 (1) R(xo，Yo；刃=R(x’，刃=If(Xo+抄o+jexp(iODd／ 其中工’是投影直线到坐标系原点(o，O)的距离，一=xsinO+ycosO． 利用傅立叶变换得到 f(x，y)=f(rexp(i纠) (2) 5南P曹喇舢加xp(ira_，cos(妒一69) 其中k(co，力是R(x’，刃关于一的一维傅立叶变换．对绷旅次积分 就可以得到常用的滤波反投影和反投影滤波重建方法． 对厂(x∽进行傅立叶交换得夕@，y)，其@定义coexp(i61)=“+iv，有 厶r 说明k(co，回是，(“，v)的极坐标表示，因此可以直接利用盖(脚，功向 李晓峰：Radon变换在高维波动方程反问题中的应用 v 直角坐标映射来构建，(“，v)．由于离散情况下，极坐标的取样点和直角 坐标的取样点不完全对应，采用了加大直角坐标系的网格密度的方式提 高数据的利用率，并减少误差．这种重建方法的速度能够提高一个数量 级． 另外由于Radon变换是以他∽为未知量的不适定的线性方程组，可 以利用代数迭代方法重建(ART)． 波在介质中传播时，由于受介质结构的影响，它的振幅、相位和速度 都受到相应影响，类似于沿传播路径对介质的投影过程，因此有可能把 Radon变换推广应用到波动力学反问题上．矿√ 建立一个与CT相似的二维模型：在X方向上【-1，1]之间且Y方向充 分长的区域Q内充满已知速度系数的均匀介质，位于以原点为圆心的单 位圆内的介质Q，的速度系数未知，但与周围的均匀介质相近．同时在x 方向的两个边界为吸收型边界条件．旋转Q-，在X=一1边界施加一平面 冲击，在X=1边界测量到达波． f振幅“(x，Y，t)与介质速度系a(x，y)满足双曲型方程 旦：!!；i；匕尘=a(工，y)V-(a(x，y)V“(x，y，f)) (4) 在区域为Q的波动场中，速度系数a(x，y)的分布为： 2 {。(J，y1a。，当?：：?2冀-． (5)。 【a(x，y)=ao+Ⅱ’，当(工，_y)∈Q1 其中a。为常数，a’未知．分析满足al=国，的摄动情况，即 0占l，-且al(x,y)=O(ao) 当s=O时．波在均匀介质中传播，则可以得到平面波解 “o(x，Y，t)=5(aot—X—1) (6) 当蒯时，分析双曲型方程的奇性，把方程的解可以写为以下的形式 (7) u(x，y，t)=以烈工，Y，f))+y(x，y)日(妒(工，Y，f))+··‘ vI 摘要 这里奇性函数妒(J，y，t)的零点的传播反映着波的传输特性． 把上式代入双曲线方程，两边扩’(西和扩(力的系数(奇性)相同，得 譬=±Ⅱ((经2+(争2)…