逆散射多次迭代毕业论文.doc

多次迭代的散射估计第一部分；理论原因A .J.Berout and D.J.Verschuur 摘要回顾表面相关多次波问题可以通过利用所谓的反馈模型来给定。从结果方程可以总结出，利用提出的方法不需要地下的任何信息。但是资源探测和表面的反射率需要特殊对待。这些信息通过表面因子被定量化，这个操作是多次波去除估计的一部分。表面相关多次波去除算法已经通过诺依曼连续和迭代方程明确地表达出来。诺依曼连续公式需要一个非线性最优化的进程用来表面算子，然而迭代方程用公式需要大量的线性最优化。迭代公式也有利于容易和其它的去除多次波方法相互协调。一个用来去除内部多次波的运算法则也需要被提出。这个运算法则是表面相关多次波方法的扩充。去除内部多次波需要表面和多次波产生层的上边缘大量的模型速率。第二部分（在这个问题上也发表了）这个提出来的运算法则已经被数值模拟实验和野外数据所成功证实。 引言直到现在，地震成像还是基于原始的反射能量。因此，多次波反射的衰减已经成为了一个很重要的课题。在过去，多次波的衰减改善主要是利用通过两个不同的方法来实现的：统计的预测错误滤波和确定的时间差滤波。统计预测错误滤波算法利用合适的多次波尾部可能会被很好的预测从它比更早的到达原始反射通过利用关联函数（Robinson,1957）。对于一维介质，预测错误滤波是很成功的，尤其是在线性拉东域（Taner,1980）.但是，对于更多的复杂介质反映，预测错误滤波是很少有效果的。换句话说，较长的预测误差滤波器可能会扭曲原始的反射事件。第二种方法中的算法，确定性时差滤波，利用原始反射的信息可能有一个有意义的更小时间差与多次波干涉相比较而言。对于介质大积极速度梯度，时间差滤波可能非常成功，尤其是在抛物线拉东域（Hampson,1988）和双曲线的拉东域（Foster and Mosher 1992）。但是在小的速度梯度或速度逆转的情况下，介于基元和多次波之间的时差区分将会很小，特别是存在一个小的补偿，时差滤波同样将会严重减弱原始的反射。也就是说，在复杂介质反映的双曲线或抛物线拉东变换应用的情况下，可能不再有意义。在最近的70年代，相关的较强表面多次波引起了特别的关注。Kennett（1979）描述了一个一维介质关于表面相关多次波的正演模型，并提出了一个一维反演方案。但是，对真实数据而言，他的算法包含了数据的取得和表面信息太过于简单的而不能成功。Riley和Claerbout(1976)描述了二维介质关于表面相关多次波的正演模型，但是他们仍然没有得到一个合适反演方案。一个多维的用来去除表面相关多次波的反演算法在Berkout(1982)被提出。他的公式基本的是任一地下模型都能够被计算，总之，源和接收的信息都被考虑。后者成功倾向于一个绝对必要的野外数据。在他的博士论文，Verschuur(1991)成功的证明了野外数据的逆源小波作为表面相关多次波去除算法的一部分能够通过利用最小二乘法准则被精确预测。对于震源子波估计和表面相关多次波的去除理论能够被用来去除内部的多次波（Berkout，1982)：向下外推后的炮点记录虚拟数据采集表面，所有到达这个表面的相关多次波都可以利用表面相关多次波算法去除。它还提出让内部的多次波叠前偏移的部分去除过程。在这部分（第一部分）多次波的估计和去除（表面相关多次波和内部多次波）理论被提及，并且以一个快速的收敛迭代方式被提出。在第二部分（Verchuur and Berkout ,1997）这个被提到的成功算法将会利用数值模拟和野外数据被证实。 地震数据的正演模型使用一个详细的隐藏的操作算子（Berkout，1982），这个正演震源子波模型能够被下面的公式给出和在连续的情况下，表达式（1a），（1b）和（1c），可用积分方程明确的表达；然后在离散的情况下，表达式（1a），（1b）和（1c）可以用矩阵方程表示，P和P+各自地定义为向上和向下移动的波场，W和W+各自地定义为向上移动和向下移动传播算子，R和R+各自地定义为向上或向下传播的入射波场的反射算子。S+（z0）定义为向下传播在 z0表面获得的震源波场，包括震源模式和虚震源。 表达式（1a），（1b）和（1c）所涉及到的一个暂时傅里叶变换（或者拉普拉斯变换）组成，并且强调了二维和三维波传播。对于一个不断更新的波场理论的算子公式的讨论，读者可以查阅Berkout，（1993）。 如果我们忽略内部多次波的散射，那么 R-(zm)对于m1时被认为是从零开始的，方程（1a），（1b）可能会被替代被给出通过式子（1c）。 如果我们忽略所有的多次波散射（内部的和表面相关多次波），那么R-(z0)也是从零开始的方程（2a）可以简化为 图1，在z=z0表面时地震模型。向下传播的波场组成了震源子波波场S+（z0）和一起向上反射的波场R-（z0）P-（z0），算子D-（z0）增加了直接接收的记录数据：P（z0）=D-（z0）P-（z0）图1，总结了在图表帮助下的正演模型。算子X(z0,z0)在自由压力表面下的假设下，地下（zz0）所呈现的转换函数。算子X（z0,z0)在自由反射界面上的半空间（zz0）所呈现的转换函数，没有内部多次波反射我们可以写成图1，清楚地表达了反射算子R-(z0)的存在引起了环形反馈。这个环形产生了相关多次波。在图1中，算子D-（z0）已经包括了探测模式和检波器鬼波所呈现的影响矩阵P(z0)是在z0的地震测量。每一列的W-和W+定义为每一单炮的反映，其元素可以写成H(w)定义为一个区分的有限簇（任意的层状影响）被认为是从炮点（xj,z0)到反射点（zl,zm)的旅行时，被认为是反射点（xk,zm)到探测器点（xi,z0)的旅行时。这里的，是当可适用时相关的扩展几何因子，应该包括合适地焦散曲线影响。接下来，薄层多次波将包括W-和W+，意思是它们将会按照统计的方法对待。在作者看来，在进行叠前偏移处理时薄层影响就应该被去除，在本文中将不会被讨论。 每一列的R（zm）定义为反射因子R(zm)和反射元素Rkl(zm）。R(zm)在线性拉东变换域中是单独的在网格点（xl,zm)角度反射信息作为函数当层的参数p。如果R（zm）是一个对角矩阵，那么它在zm层的反射信息是单独的角度反射。 最后，值得一提的是，前述的频率域表达式同样可以扩展到时间域。意思是说所涉及到的标量元素可以变为时间函数，例如所有的要素都需要被暂时卷积所取代，考虑到在时间域的一个炮点记录模型，时间域的公式表达*的意思是在一个矩阵里面两个矩阵元素不是相乘，而是褶积。接下来，频率域的公式将会被用到。 多次波的产生和消除 从表达式（1c），（2a），（2c）和（2d）可以很容易地看出我们所要记录的反射数据（图1） P。和M。分别表示了原始的和多次波响应。 表达式（4）表示了震源子波波场S+（z0）产生的原始响应第二次震源波场产生的多次波响应将方程（2d）,(4)和（5）联合起来组成一个重要的结果其中A(z0)是表面算子 表达式（6a）和（6b）在本章中扮演了一个中心作用。接下来所看见的方程（6a）所呈现的是一系列第二类积分方程，这种积分方程的样式是任意传统的物理系统包括环形迭代。在逆散射文章中，第二类积分方程同样以 方程的形式被提及。 从方程（6a）可以看出，多次波的产生和消除的基本表达式可以按照如下的形式表示；多次波的产生（图2）或或或，也就是诺依曼连续第n次的的反射能量是在图1的闭合圆环旅行了n时间，导致产生了第n次的多次波。考虑到所有的多次波能量都是有限的，可以总结到，在物理条件下，诺依曼方程是收敛的。也就是说，对于有限的时间重合，仅仅有一个有限的特征项被涉及。 图2，（a）表面相关多次波的产生流程，（b）表示相关多次波的去除流程多次波的消除（图2b） 或或 或，也就是诺依曼连续 方程（8b）也可以写成其中 考虑到多次波的温纳预测滤波，可以被作为基本模型预测滤波器。对于的应用，仅仅是表面算子需要被估计。按照方程（6b），假如我们忽略直达波，那么可以被认为是一个对角矩阵。最后，我们假设在地震测量的过程中震源和探测表面没有发生任何改变，那么被认为是一个一定比例的单位矩阵，这个比例因子将会由傅里叶变换（拉普拉斯变换）的震源子波倒转的成分所决定。目前的理论包括任何一种可能的情况。接着Verschuur(1991),指出矩阵的元素可以通过下面式子的最小化估计 其中表示所有地震宽带频率的平均值。方程（10a）和（10b）的的估计需要非线性最新优化算法。 利用方程（6a）在物理条件下考虑到其收敛性和唯一性，我们可以总结到 因此，复杂的表面相关多次波，多频段的多次波预测等同于原始响应的展开褶积。方程（10c）表示，尽管可能被许多滤波器系数所定义，最佳化流程仅仅需要对一个相关少数量的算子系数（矩阵元素）用公式表示。这个法则是我们考虑多频段时的基础模型的优势。 迭代法去除多次波从表达式（9a）,(9b)可以看出，表面相关多次波的去除同样可以写成循环的多频段滤波处理。其中和 （11c）方程（11a）很容易的就可以被证明：利用诺依曼连续的递归处理。因此，在方程（11a）可以通过在诺依曼连续中取n次的近似值。方程（11a）,(11b),(11c)有如下的特征；在许多的情况下方程（11c）一般选择能够更接近于解决方法的最初估计。举个例子，在拉东变换基础上去除多次波方法能够用最初的估计。 实际长度不需要被选择为所有路径的长度就像方程（11b）那样。举个例子，表达式（11c）中可以被一个缩短的长度替代，并且多次波的产生仅仅为那一部分。表面算子可以通过重复迭代而最优化 n=1,2,3N （12a）其中方程（12）的优点在于考虑到方程（10），非线性最优化进程都可以被线性的最优化进程所取代。利用方程（12b）表达式（11a）,(11b)和（11c），能够改写成如下的公式形式矩阵方程（13a）和（13b）成了接下来的空间-率域算法的基础或，忽略直接到达表面的或，同样忽略在表面横向的变化在表达式（13d）,(13e)，（13f）中，标量f分别对应了矩阵中的元素。在方程（13c）中，表示了矩阵原始数据中的初始估计。看（11c）的基本表达式，最简单的选择就是通过给出。但是我们同样可以在任意地拉东滤波后让多次波产生部分作为数据的一部分。 理论的介绍让我们先来讨论最简单的情况，即平面波在水平介质中的传播。那么这个矩阵在（13）的迭代算法中将会变为一个标量。对于一个水平面波仅一个反射界面的情况下，我们可以将在表面向上的波场写成如下的形式： 其中是双程旅行时，r是反射界面的反射系数，因此，利用，这个对于多次波的散射便可以通过下面的式子进行证明； 基本关系（见方程（6a）， 总响应的诺依曼连续（见方程（7b）, 对于预测滤波的诺依曼连续（见方程（9b）, 图3证明了我们所选用的最简单迭代算法（13f）在一个简单的例子下（点源，水平反射界面），利用仅仅只是需要很少的迭代。在第二部分，数值模拟实验和野外数据实验将会被讨论。 选用一种参考介质让我们再回顾就转换函数基本表达式（6a） 表示在一个已知的有自由应力表面的介质的转换函数，和表示在虚拟的自由反射界面介质中的转换函数。每一列的和定义为每一单位震源响应（格林函数）。我们将称这个具有自由反射界面的虚拟介质为参考介质。现在，我们再次写出（13a）的转换函数形式，然后我们将发现在我们的迭代方法中参考介质随着每次迭代反复的变化着；或，假设单独角度表面反射，或，在海上的情况下，在表达式（14c）,(14d),(14e)中，标量，分别地对应了矩阵和中的元素。图4，用了一个简单的表格说明了方程（14d）。最近，和在相互定理的帮助下，通过比较波场在参考介质（自由反射）和媒介（自由应力）已经派生出相同的表达式。此外，最近利用方程同样得出了以上的参考介质和媒介同样得到相似的表达式。尽管不同的方法都有同一个目的，提高精度，我们喜欢在图1上的基于反馈图来引出的原始数据因为面临着物理分类的问题；此外它清晰的给我们指出了震源和探测信息的影响。最后，有意思的是，根据提到的方程（13）的算法，随着参考介质的变化迭代的形式同样在变化。方程（13a）,(13b)和（13c）迅速的收敛，即使是我们开始以一个均匀的参考介质，意思是等于零。（a） (b)原始和层间多次波 一次迭代 二次迭代 三次迭代（c） 总响应 一次迭代 二次迭代 三次迭代 图3，简单的地下模型。（b）表面相关多次波的产生的迭代过程（c）在空间-频率域多次波消除的迭代过程。多次波的去除过程包括震源子波的最优化。快速收敛。 层间多次波假如代表的是第一次多次波产生的内部多次波的边界（如，海底），可能广义的认为那么与多次波去除有关的边界可以按照下面的式子定义。仅用一次简单记数选取一次迭代；其中表示的是向下外推预处理表面有因果关系部分（t0）的数据, 在一些稳定的情况下（16a） 在一些稳定的情况下（16b）和符号指的是与有关的没有多次波的反射数据，指的是与有关的没有多次波的反射数据，同样是的。如果我们将方程（15a）和（15c）联立，我们可以写出与边界有关的多次波；或在方程（17a）,一列定义为检波器在，炮点在的单炮记录。一行定义为一个检波器组合在接收，而震源在。因此，通过方程（15）和（17）的组合，与有关的内部多次波的去除算法可以通过下式给出（一次迭代）与表面相关的算法取其最近相似性（一次迭代）因此，如果和大模型之间的去除已知，那么和可以通过向下外推而被计算出来，与表面相关方案类似，的最小平方消除在多次波去除最优化的进程中被证实。有意思的意识是，这个消除理论包括了和之间的色散影响（例如，薄层多次波）。图6，阐述了表面相关多次波和内部多次波之间的区别。在表面相关多次波方法中，一个普通的震源-检测位置在，定义为一个表面多次波的反射点（的一列的特征），在内部的一个普通的震源-检测位置在定义为一个内部多次波的反射点（一列的特征）。最后，图7,表示出了以图3为例子的数值模拟的内部多次波的去除流程。图4，表面相关多次波的产是由表面在k点的反射，。这里的n表明这个过程被认为是第n次时间。真实的介质是地下的自由应力界面。参考介质收敛到相似地下自由反射界面。参考介质随着迭代的每次迭代而变化。 真实炮点记录的多次波去除利用诺依曼连续（8b），可以很容易的证实关于多次波去除的迭代公式可以写成两个版本。在前述中，我们仅仅考虑第一个公式。但是，如果我们利用一个综合的算子将数据转换成一个真实的区域炮点记录（例如，区域震源响应），那么这个迭代过程可以很容易利用区域炮点记录，同样的我们也可以用第二个公式；其中 被认为是综合算子。检验方程（20a）表示了关于区域炮点记录的多次波去除仍然需要利用原始炮点记录。同样地，对于内部多次波（一次迭代）其中在图3中，利用数值模拟，区域炮点记录已经被合成用于斜平面波，并且多次波的去除流程（表面相关多次波和内部多次波）已经被用于区域炮点记录。结果将在图8上显示出来。a,接收和震源在 b,接收和震源在 C,接收在，震源在 图5，根据方程（17a），联合和与反射边界有关的内部多次波 总结多次波的反射已经被细分为表面相关多次波和层间多次波一个反馈模型已经被用于量化就第二类积分方程形式而言的多次波在去除表面相关多次波理论中，两种介质应被考虑自由应力表面的实际介质和自由反射界面的参考介质图6，表面和内部多次波去除过程说明，（a）通过表面数据的褶积来预测多次波（b）通过向下延拓插值的表面数据褶积来预测内部多次波（a）总响应 预测的原始和所有的内部多次波 预测的表面多次波（z=z0(b) 原始响应和所有多次波 预测的原始和内部多次波（zz2 内部多次波（z=z1 图7，表面相关多次波和内部多次波的去除结果说明；（a）表面相关多次波的去除（b）与海底有关的层间多次波去除（z=z1）(a) (b) (c) (d) 图8，根据方程（20a）和（21a）利用区域炮点记录去除表面相关多次波和内部多次波（图3的斜平面波地球模型的响应）a)单炮记录的多次波b)区域炮点记录的多次波c)与z0有关的多次波去除d)与z0和z1有关的多次波去除表面相关多次波的去除过程是将实际介质中的响应转换为在参考介质中的响应。表面相关多次波的去除不需要任何地下信息；仅仅需要表面算子被估计。在海上的情况，表面算子定义为震源-检测波的转置。这个波的估计被认为是表面相关多次波的去除过程中最重要的亮点。表面相关多次波的去除算法可能通过诺依曼连续或迭代方程而公式化。诺依曼公式的表面算子需要非线性为最优化的过程；然而迭代公式需要线性的最优化。这个迭代公式更容易和其它多次波去除方法整合。从表面相关多次波的去除理论算法同样也使用层间多次波的去除。两种算法之间的不同是在表面相关多次波中，多次波的产生的表面算子与多次波的去除算法将被一个向下连续，在内部多次波的去除的多次波的产生边缘算子所替代。意思是对于在表面和多次波的产生边缘的大模型的内部多次波的去除 一定要知道。关于内部多次波去除的表达式，所涉及的最优化过程可能是一个色散的影响（例如，薄层多次波的产生）增加了震源子波从多次波产生的边缘的传播，多次波的去除算法同样适用区域炮点记录。袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇