CN103675926A - 一种航空瞬变电磁资料的电导率-深度转换方法 - Google Patents

Abstract

本发明公开了一种航空瞬变电磁资料的电导率-深度转换方法，该算法分计算视电导率、计算视深度和估算成像深度三个步骤。首先提出了一个关于均匀半空间模型电导率与响应数据的变换之间均呈一一映射的函数公式，依此建立查询瞬时视电导率值的数据函数表，克服了直接由感生电动势查询视电导率出现的非唯一性问题，也没有出现视电导率无解的问题，获得可靠的视电阻率值。其次将均匀半空间中感应电场最大幅值的1.28倍作为视深度，并通过一个预先建立的数据函数表快速查询得到，提出一个成像深度的近似计算公式，用2个相邻延时期间的视深度变化的基础上定义模型的层厚度，从而得到电导率-深度转换结果。

Description

一种航空瞬变电磁资料的电导率-深度转换方法

技术领域

本发明给出一种航空瞬变电磁的电导率-深度转换技术，属于时间域航空电磁资料的资料处理解释方法范畴。

背景技术

航空电磁勘查方法具有探测深度大、测量精度高、勘查效率高、工作成本低的技术特点，在找矿、大面积地质调查、水文地质探测与含水岩层推测、非金属矿物勘查等领域中有广泛的应用。由于航空瞬变电磁资料量大质差、对处理速度的要求很高，而电导率-深度成像(Conductivity-Depth Imaging，CDI)方法能够快速地将观测数据处理为地下电导率-深度剖面，可以快速圈定目标异常或作为反演的初始模型.与反演方法相比(Ellis，1998；Christensen，2009；Cox，2010；毛立峰，2011)，CDI方法计算速度快，不需要初始模型，是时间域航空电磁资料解释的第一步和必不可少的一步。目前已有多种时间域航空电磁的CDI方法，如Macnae等(1991)提出的像源扩散成像方法，适用于发射脉冲是阶跃电流脉冲的航空瞬变电磁资料解释。基于此法，Wolfgram和Karlik(1995)将观测数据处理成阶跃响应数据再做上述成像处理，从而实现了发射波形为半正弦脉冲的GeoTEM系统资料的近似成像。Liu等(1993)的CDI方法应用了基于水平电导薄板的近似反演算法，他来进行相邻时间道的响应数据进行的一维线性反演，其中的正演用薄板的阶跃响代替，从而提高了计算速度。Fullagar(2001)用均匀半空间中最大电流深度近视为成像深度，半空间模型的电导率取为视电导率，这种方法估算成像深度较直接，但结果不是很准确，对低阻层的深度估计值偏低，而对高阻层的深度估计值偏高。Huang(2008)提出了AeroTEM系统资料的假层半空间模型CDI方法，这种方法的成像速度快，通过插值和查询的计算方法，提出的τ-α数据表查询视电导率的方法很有效，可以完成在发射线圈高度有变化下的电导率-深度成像。但由于成像深度估算方法很困难，Huang的成像深度以经验关系曲线图形式给出，并没有给出成像深度的算法，另外对本发明针对的航空瞬变电磁系统参数，由Huang的τ-α数据表查询视电导率并不唯一，在高电导率段不同电导率对应的τ-α曲线是交叉的，在模型电导率值较高时得到的视电导率值就不可靠了。航空瞬变电磁法，还有地面瞬变电磁法中，由感生电动势数据直接计算视电导率值时，感生电动势是关于均匀半空间模型电导率是双值函数，如何得到唯一可靠的视电导率是很困难的，目前还没有好的解决办法；另外由于介质非常复杂，成像深度不能仅仅采用视深度，如上一延时的视深度可能大于下一延时的视深度，无法得到合理的成像结果。

本发明提出自己的算法克服了上述固有的困难问题，提出了一个关于电导率和观测数据的变换之间的一一映射函数公式，从而建立查询表获得可靠的视电导率，在此基础上，申请人给出自己的视深度定义方法，并提出简单有效的成像深度计算公式，形成电导率-深度转换结果，并将这些新的发明申请专利保护。

发明内容

本发明的目的在于提供一种航空瞬变电磁多分量资料的快速成像技术，将观测的x和z分量感生电动势响应资料快速转换成电导率-深度剖面电性分布结果。计算简单，不需要反演迭代过程，可以在飞行作业时进行实时成像。

为达到上述目的，本发明采用的方法包括以下步骤：

(1)建立视电导率查询数据函数表

提出了一个关于均匀半空间模型电导率与响应数据的变换之间均呈一一映射的函数公式：

${\mathrm{\γ}}_i\left(\mathrm{\σ}\right)=\left[\ln \right|\frac{V_x^i\left(\mathrm{\σ}\right)-V_x^{i+k}\left(\mathrm{\σ}\right)}{V_z^i\left(\mathrm{\σ}\right)-V_z^{i+k}\left(\mathrm{\σ}\right)}|-2.5*\ln \frac{V_x^i\left(\mathrm{\σ}\right)*V_x^i\left(\mathrm{\σ}\right)+V_z^i\left(\mathrm{\σ}\right)*V_z^i\left(\mathrm{\σ}\right)}{V_x^{i+k}\left(\mathrm{\σ}\right)*V_x^{i+k}\left(\mathrm{\σ}\right)+V_z^{i+k}\left(\mathrm{\σ}\right)*V_z^{i+k}\left(\mathrm{\σ}\right)}]/\ln \frac{t_{i+k}}{t_i}---\left(1\right)$

其中，i是采样时刻点的序号，i＝1，2，...，n-k，n是采样时间道总个数，正整数k是给定的时间道间隔个数，取值1～4。是

表示电导率为σ的均匀半空间的第i个断电后的采样时刻点t_i的x分量感生电动势响应。给定10^-4～10²S/m范围内的呈等比分布的m个电导率值，分别计算各电导率值对应均匀半空间模型响应分量和(i＝1，2，...，n；j＝1，2，...，m)。代入公式(1)得到数据函数表γ(t_i，σ_j)，(i＝1，2，L，n-k；j＝1，2，L，m)。

(2)计算视电导率

设观测的航空瞬变电磁响应数据是

和

(i＝1，2，L，n)，对它进行变换：

其中，i＝1，2，...，n-k，根据Υ_i在数据函数表γ(t_i，σ_j)中反插值计算得到的电导率值即为第i个采样时间道对应的视电导率值。

(3)建立视深度数据函数表

将延时t_i的均匀半空间模型中最大电场对应的深度的1.28倍定义为视深度d(t_i)。给定10^-4～10²S/m范围内的呈等比分布的m个电导率值，分别模拟航空瞬变电磁场在各电导率值对应均匀半空间地下介质中的传播，并记录最大电场值对应的深度，从而获得电导率-视深度数据函数表

，(i＝1，2，L，n；j＝1，2，L，m)

(4)计算成像深度，形成电导率-深度转换结果

根据步骤(2)得到各个采样延时的视电导率，并由数据函数表

插值计算各延时对应的视深度

，为得到更为准确的成像深度D_i，对视深度做变换，第k个延时的D_i的定义如下：

$D_i=\left\{\begin{array}{cc}d\left({\mathrm{\σ}}_a^i\right)& i=1\\ D_{i-1}+[d\left({\mathrm{\σ}}_a^i\right)-d\left({\mathrm{\σ}}_a^{i-1}\right)]*\frac{{\mathrm{\σ}}_a^{i-1}}{{\mathrm{\σ}}_a^i}& i>1\end{array}\right.---\left(3\right)$

这样由所有延时的视电导率

和成像深度D_i构成的层状模型形成电导率-深度转换结果曲线。

本发明利用公式(1)建立具有一一映射关系的电导率-均匀半空间x、z分量的感生电动势响应数据的变换数据函数关系，从而由观测数据变换公式(2)通过反插值方法获得可靠的视电阻率。然后通过查询视电导率对应的视深度，由公式(3)快速得到对应的成像深度，从而将观测数据快速转换为电导率-深度剖面数据。

以发射梯形脉冲的航空瞬变电磁资料为例，系统参数假定如下：飞行高度ht为120m，接收线圈高度hr为70m，收发线圈水平距r为120m。发射磁偶极子源的磁矩M为47.124万Am²，接收线圈有效面积为100m²。系统发射周期性双极性脉冲，其梯形单脉冲的供电时间为10ms，上升沿和下降沿时间均为2ms，梯形的峰值平稳时间是6ms，断电时间长度为10ms，峰值电流强度为300A。在断电期间测量x和z分量的感生电动势响应数据。本发明的具体原理如下：

先计算用于视电导率查询的数据表的均匀半空间模型正演公式，假定离散化的均匀半空间电导率为σ_i，i＝1，2，...，m，则对应的均匀半空间模型的x、z分量航空瞬变感生电动势下阶跃脉冲响应为：

$V_x({\mathrm{\σ}}_i,t_k)=-{\mathrm{\μ}}_0{S}_r\frac{M}{4\mathrm{\π}}{L}^{-1}\left\{{\∫}_0^{\∞}\right[e^{-\mathrm{\λ}(\mathrm{ht}-\mathrm{hr})}-\frac{\mathrm{\λ}-u}{\mathrm{\λ}+u}{e}^{-\mathrm{\λ}(h0+\mathrm{hr})}\left]{\mathrm{\λ}}^2{J}_1\left(\mathrm{\λr}\right)\mathrm{d\λ}\right\}---\left(4\right)$

$V_z({\mathrm{\ρ}}_i^0,t_k)=-{\mathrm{\μ}}_0{S}_r\frac{M}{4\mathrm{\π}}{L}^{-1}\left\{{\∫}_0^{\∞}\right[e^{-\mathrm{\λ}(\mathrm{ht}-\mathrm{hr})}+\frac{\mathrm{\λ}-u}{\mathrm{\λ}+u}{e}^{-\mathrm{\λ}(h0+\mathrm{hr})}\left]{\mathrm{\λ}}^2{J}_0\left(\mathrm{\λr}\right)\mathrm{d\λ}\right\}---\left(5\right)$

其中，μ₀为真空磁导率，

κ为波数，满足κ²＝-sμ₀σ_i，L^-1为逆拉普拉斯变换算子，J₀和J₁分别为第一类0阶和1阶贝塞尔函数。

设给定的2分量观测数据是和

(i＝1，2，L，n)，对它进行变换：

${\mathrm{\γ}}_i\left(\mathrm{\σ}\right)=\left[\ln \right|\frac{V_x^i\left(\mathrm{\σ}\right)-V_x^{i+k}\left(\mathrm{\σ}\right)}{V_z^i\left(\mathrm{\σ}\right)-V_z^{i+k}\left(\mathrm{\σ}\right)}|-2.5*\ln \frac{V_x^i\left(\mathrm{\σ}\right)*V_x^i\left(\mathrm{\σ}\right)+V_z^i\left(\mathrm{\σ}\right)*V_z^i\left(\mathrm{\σ}\right)}{V_x^{i+k}\left(\mathrm{\σ}\right)*V_x^{i+k}\left(\mathrm{\σ}\right)+V_z^{i+k}\left(\mathrm{\σ}\right)*V_z^{i+k}\left(\mathrm{\σ}\right)}]/\ln \frac{t_{i+k}}{t_i}---\left(6\right)$

其中，i＝1，2，...，n-k，根据Υ_i在数据函数表γ(t_i，σ_j)中反插值计算得到的电导率值

即为第i个采样时间道对应的视电导率值。

按上式的定义，对两分量感生电动势观测数据v_x、v_z做变换，计算各校正时刻T_k时的幅度比值数据：

在数据函数表γ(t_i，σ_j)中反插值计算得到的电导率值

即为第i个采样时间道对应的视电导率值。由于γ_i(σ)是σ的一一映射，反插值计算视电导率得到的结果是唯一的，克服了常规方法由观测感生电动势数据查询视电导率出现非唯一值问题。为说明这一问题，图1为断电后16个延时的均匀半空间模型的垂直感生电动势响应与模型电导率之间的函数关系曲线，其中图a和图b分别是x和z分量感生电动势响应Vx和Vz对应的曲线，图c是两分量数据得到响应幅度

对应的曲线，而式(6)得到的响应数据的变换曲线是电导率的单值函数，使电导率分布范围很宽，反插值查询视电导率值就无非唯一性问题，故可得到可靠的视电导率值。σ-γ关系曲线见图2所示。

可靠的视电导率也是获得较合理的成像深度的前期条件，有了各延时的视电导率，接着要计算各延时的成像深度。这里提出一种简单的深度估算方法，先定义视深度，可用数据表查询得到，再由视深度定义各延时的成像深度，形成层状介质模型，得到最终的电导率-深度转换结果。

计算电导率取值为视电导率值

的均匀半空间模型的地下介质中的电场，将延时t_i的地下最大电场对应的深度的1.28倍定义为视深度

。建立一个视深度数据查询表σ_a-d来加快计算速度。图3是10个采样延时的视电导率与视深度之间的关系曲线，各视深度曲线均有随电导率增加而单调下降的特点，故容易通过视电导率值查询视深度值。由数据函数表

插值计算各延时对应的视深度

，为得到更为准确的成像深度D_i，对视深度做变换，第k个延时的D_i的定义如下：

$D_i=\left\{\begin{array}{cc}d\left({\mathrm{\σ}}_a^i\right)& i=1\\ D_{i-1}+[d\left({\mathrm{\σ}}_a^i\right)-d\left({\mathrm{\σ}}_a^{i-1}\right)]*\frac{{\mathrm{\σ}}_a^{i-1}}{{\mathrm{\σ}}_a^i}& i>1\end{array}\right.---\left(8\right)$

这样由所有延时的视电导率

和成像深度D_i构成的层状模型形成电导率-深度转换结果曲线。

第一个例子是为了检验模型电导率超过1S/m的模型的电导率-深度转换效果，设三层介质电导率分别为2S/m、5S/m和2S/m，厚度分别为10m和5m，为了获得有效探测深度，降低系统的基频，将断电时间设为20ms。图4是真实模型、电导率-深度转换结果曲线，视电导率的最大值为3.17465S/m，最小值为1.77979S/m，可见，视电导率查询结果是可靠的，查询中也没有出现无解的情况。电导率-深度转换结果已经较好地反映出中间高导层模型特征了，视电导率最大值对应的深度为13.9504m，接近高导层下界面深度。

第二个例子是二个四层KH模型的电导率-深度转换试算，模型的电导率分别为0.05、0.002、0.5和0.002S/m，模型第一层厚度为40m，第三层是厚度为20m的薄层，其深度分别是100m和200m，数据的采样时间是从0.01ms至10ms。图5为两者的电导率-深度转换结果曲线，均出现KH型特征，薄层段的成像结果的电导率值高于其它层段，第二层的电导率解释结果低于两侧的低阻层段，但高于基底层段，这是它受两侧影响的结果。由于第二层是高阻层，层厚度大于薄层深度为100m的模型的层厚度，所以薄层深度为200m的模型的解释电导率值相对低一些。

第三个例子为均匀半空间中低阻凸晶体状目标体模型的电导率-深度转换，用一维响应数据合成剖面响应数据作为成像数据。设围岩电导率为0.02S/m，目标体的电导率为0.1S/m。凸晶体垂直断面的中心深度及最大厚度分别为200m和25m，凸晶体的两个角点的x坐标分别为-20m和20m。采样时间为从0.01ms至4ms。图6为模型、几个延时的响应幅度A曲线和电导率-深度转换结果。电导率-深度转换结果图上的异常体存在的横向位置上，均出现高电导率值异常区域，最大异常电导率值的成像深度在200米深度左右，与目标体的真实深度很接近。

本发明的核心在于：(1)提出一个关于电导率呈一一映射函数关系的响应数据变换的函数公式γ(σ)，建立数据表快速查询视电导率值，用观测的x、z两个分量的感生电动势的变换查询得到可靠的视电导率值。(2)定义视深度为电场在下半空间模型中的最大电场值对应深度的1.28倍，并建立视电导率-视深度的数据查询表，加快成像的计算速度。(3)根据视电导率值查询视深度值，提出了成像深度的计算公式，形成CDT结果模型。基于上述特点，本发明具有如下优点：

(1)出了一个能够在很宽的电导率范围(10^-4～10²S/m)内快速查询得到可靠的视电阻率值的方法，用了一个关于均匀半空间模型的x、z分量响应数据的变换公式γ(σ)，它与均匀半空间电导率σ之间具有一一映射函数关系，从而克服了直接用响应数据进行视电导率查询出现的非唯一性问题。

(2)在获得可靠的视电导率前提下，用简单有效的方法实现困难的成像深度估算，先用数据表查询视电导率值对应的最大电场深度值，视深度定义为它的1.28倍，CDT结果模型的层厚度由提出的基于相邻的采样时间段内的视深度变化的简单公式给出，将困难的成像深度计算简单化，而且效果较好；

(3)该算法理论较简单，成像思路可以扩展到其它发射波形的系统甚至其它航空瞬变电磁系统的资料解释中，并无本质上的差异；

(4)计算过程仅有简单的查询、插值和一次正演计算，计算速度快，符合航空瞬变电磁要求，使得航空电磁在飞行中的实时处理成为可能。

附图说明

图1为均匀半空间模型的垂直感生电动势响应与电导率关系曲线，其中图a和图b分别是x和z分量感生电动势响应Vx和Vz对应的曲线，图c是两分量数据得到响应幅度A对应的曲线。

图2为14个时间道对应的γ与均匀半空间模型电导率的关系曲线(k＝2)。

图3为10个延时的视深度与半空间模型电导率之间的关系曲线。

图4为高导三层状模型的电导率-深度转换结果。

图5为四层模型响应数据的电导率-深度转换结果。其中，图a是薄层深度为100m的模型的电导率-深度转换结果，图b是(a)均匀半空间中低阻凸晶体目标体电导率分布。

图6为二维低阻菱形目标体模型的电导率-深度成像。其中，图a是均匀半空间中低阻凸晶体目标体电导率分布，图b是几个时间道的剖面响应曲线，图c是电导率-深度转换结果。

具体实施方式

本发明分计算视电导率、视深度和计算电导率-深度转换的层厚度三个步骤。首先用提出的关于电导率呈一一映射的2个分量响应数据的变换的函数公式建立数据表，由观测数据查询得到可靠的视电导率值；其次定义均匀半空间中最大电场对应深度的1.28倍为视深度，建立视深度的数据查询表，用视电导率值快速查询得到视深度值，提出的基于相邻的采样时间段内的视深度变化的简单公式计算成像模型的层厚度，形成最终的电导率-深度转换结果模型。具体包括如下步骤：

(1)输入计算参数：航空瞬变电磁系统参数：飞行高度h0、校正高度h1、发射磁矩M、接收线圈有效面积Sr、接收线圈的高度hr、收发线圈的水平距r，接收信号的延时采样个数n及采样时间序列t_i(i＝1，2，...，n)；计算视电导率时用的两个采样数据的时间间隔个数k；接收x、z分量感生电动势信号数据v_x、v_z；电导率查询表的半空间模型的电导率计算范围，这里设计为10^-4～10²S/m，电阻率个数m、对数等间隔的查询电导率序列σ_j，j＝1，2，...，m；

(2)建立视电导率查询数据函数表。分别计算查询表各半空间模型电导率值对应的响应分量和

(i＝1，2，...，n；j＝1，2，...，m)。代入公式：

${\mathrm{\γ}}_i\left(\mathrm{\σ}\right)=\left[\ln \right|\frac{V_x^i\left(\mathrm{\σ}\right)-V_x^{i+k}\left(\mathrm{\σ}\right)}{V_z^i\left(\mathrm{\σ}\right)-V_z^{i+k}\left(\mathrm{\σ}\right)}|-2.5*\ln \frac{V_x^i\left(\mathrm{\σ}\right)*V_x^i\left(\mathrm{\σ}\right)+V_z^i\left(\mathrm{\σ}\right)*V_z^i\left(\mathrm{\σ}\right)}{V_x^{i+k}\left(\mathrm{\σ}\right)*V_x^{i+k}\left(\mathrm{\σ}\right)+V_z^{i+k}\left(\mathrm{\σ}\right)*V_z^{i+k}\left(\mathrm{\σ}\right)}]/\ln \frac{t_{i+k}}{t_i},$ 得到数据函数表γ(t_i，σ_j)，(i＝1，2，L，n-k；j＝1，2，L，m)。

(3)建立视深度数据查询表。由预先模拟计算的各延时的均匀下半空间中最大电场对应深度值的1.28倍(即视深度)，建立视电导率-视深度数据函数表

。

(4)对观测数据做如下变换：

用各延时的变换后的观测数据，通过数据表γ(t_i，σ_j)，以γ_i(σ)为变量，反插值计算得到的电导率值为对应的瞬时视电导率值。计算中采样三次样条插值。

(5)根据步骤(4)得到的视电导率和视电导率-视深度数据表，插值计算视深度值，并由数据函数表

插值计算各延时对应的视深度

，插值算法按三次样条插值方法进行。最后按如下公式计算第i个延时的成像D_i：

$D_i=\left\{\begin{array}{cc}d\left({\mathrm{\σ}}_a^i\right)& i=1\\ D_{i-1}+[d\left({\mathrm{\σ}}_a^i\right)-d\left({\mathrm{\σ}}_a^{i-1}\right)]*\frac{{\mathrm{\σ}}_a^{i-1}}{{\mathrm{\σ}}_a^i}& i>1\end{array}\right.---\left(3\right)$

这样由所有延时的视电导率和成像深度D_i构成的层状模型形成电导率-深度转换结果曲线。

Claims (1)

1.一种航空瞬变电磁多分量资料电导率-深度转换技术，将观测的x和z分量感生电动势数据快速地转换为地下介质的电导率-深度的成像剖面，其特征在于该方法包括以下步骤：

(1)建立视电导率查询数据函数表

提出了一个关于均匀半空间模型电导率与响应数据的变换之间均呈一一映射的函数公式：

${\mathrm{\γ}}_i\left(\mathrm{\σ}\right)=\left[\ln \right|\frac{V_x^i\left(\mathrm{\σ}\right)-V_x^{i+k}\left(\mathrm{\σ}\right)}{V_z^i\left(\mathrm{\σ}\right)-V_z^{i+k}\left(\mathrm{\σ}\right)}|-2.5*\ln \frac{V_x^i\left(\mathrm{\σ}\right)*V_x^i\left(\mathrm{\σ}\right)+V_z^i\left(\mathrm{\σ}\right)*V_z^i\left(\mathrm{\σ}\right)}{V_x^{i+k}\left(\mathrm{\σ}\right)*V_x^{i+k}\left(\mathrm{\σ}\right)+V_z^{i+k}\left(\mathrm{\σ}\right)*V_z^{i+k}\left(\mathrm{\σ}\right)}]/\ln \frac{t_{i+k}}{t_i}---\left(1\right)$

其中，i是采样时刻点的序号，i＝1，2，...，n-k，n是采样时间道总个数，正整数k是给定的时间道间隔个数，取值1～4是

表示电导率为σ的均匀半空间的第i个断电后的采样时刻点t_i的x分量感生电动势响应；给定10^-4～10²S/m范围内的呈等比分布的m个电导率值，分别计算各电导率值对应均匀半空间模型响应分量

和

(i＝1，2，...，n；j＝1，2，...，m)；代入公式(1)得到数据函数表γ(t_i，σ_j)，(i＝1，2，L，n-k；j＝1，2，L，m)。

(2)计算视电导率

设观测的瞬变电磁响应数据是

和

(i＝1，2，L，n)，对它进行变换：

其中，i＝1，2，...，n-k，根据Υ_i在数据函数表γ(t_i，σ_j)中反插值计算得到的电导率值即为第i个采样时间道对应的视电导率值；

(3)建立视深度数据函数表

将延时t_i的均匀半空间模型中最大电场对应的深度的1.28倍定义为视深度d(t_i)；给定10^-4～10²S/m范围内的呈等比分布的m个电导率值，分别模拟瞬变电磁场在各电导率值对应均匀半空间地下介质中的传播，并记录最大电场值对应的深度，从而获得电导率-视深度数据函数表

，(i＝1，2，L，n；j＝1，2，L，m)

(4)、计算成像深度，形成电导率-深度转换结果

根据步骤(2)得到各个采样延时的视电导率，并由数据函数表

插值计算各延时对应的视深度

，为得到更为准确的成像深度D_i，对视深度做变换，第i个延时的D_i的定义如下：

$D_i=\left\{\begin{array}{cc}d\left({\mathrm{\σ}}_a^i\right)& i=1\\ D_{i-1}+[d\left({\mathrm{\σ}}_a^i\right)-d\left({\mathrm{\σ}}_a^{i-1}\right)]*\frac{{\mathrm{\σ}}_a^{i-1}}{{\mathrm{\σ}}_a^i}& i>1\end{array}\right.---\left(3\right)$

这样由所有延时的视电导率

和成像深度D_i构成的层状模型形成CDT结果曲线；对一条飞行测线的观测数据，逐点进行上面的CDT步骤，形成一个CDT断面的成像结果。

Images