体CT锥束精确图像重建Katsevich算法的理论研究

本文主要介绍了关于新型的锥束CT重建算法——Katsevich算法。首先描述了如何引入Katsevich重建公式,同时介绍了Tam-Dannielson窗口和PI线的一系列与其他重建算法不同的重建方式,指出反投影过程与Feldkamp重建算法的反投影过程是一样的。我们可以利用这个重建算法有效地进行锥束螺旋CT重建,并得到效果理想的图像。我们也可以设计相应的并行算法,便于分析这些算法的复杂度,如：分析算法时间并且采用并行算法加快运行速度。这个方法也可以被用来改进成为其他重建公式（例如Zou和Pan的BPF型反转公式）。     

基于迭代法的同轴法X射线相衬成像的相位复原

相位复原问题是指在只有强度可以测量的情况下,从强度数据中提取出物体的相位信息。相位复原是同轴法X射线相衬成像的关键问题之一。本文采用误差递减法,通过限制物体的空间大小,结合测量到的强度数据,在物体的空域和频域都进行约束。通过反复迭代,进行带宽外推,最终实现物面和像面的相位复原。仿真实验证明了此迭代算法的收敛性,并且以塑料管为样品,在上海光源采集到同轴法X射线相衬成像的实验数据,对上述算法进行了验证。     

ANALYSIS OF SEISMICITY CHANGES PRIOR TO THE 2014 YUNNAN JINGGU MS6.6 EARTHQUAKE

According to the Region-Time-Length (RTL) algorithm,the analysis of seismicity changes prior to the 2014 Yunnan Jinggu M_S6.6 earthquake was conducted by using the earthquake catalogues about 6 and 15 years before this earthquake,respectively.When the studied period was nearly 6 years,an enhancement of seismic activity was detected around the epicenter since the beginning of 2013.The anomalies mainly distributed in the region of 22.5°~24.5°N and 99°~102°E.The range and degree of anomalies changed from small to large,and then to small chronologically.As the surface integral in respect to RTL,the physical parameter I_RTL,which could reflect the regional seismicity level,began to increase since August 2013,and then reduced after reaching the peak point.The time length from the peak point of I_RTL curve to the earthquake occurrence was 9 months.When the analyzed catalogue was nearly 15 years,the 2007 Ninger M_S6.4 occurred in the studied region.Seismicity quiescence was detected prior to the Ninger M_S6.4.Before the Jinggu M_S6.6,seismicity quiescence was detected firstly,and then enhanced activity was observed 1 year prior to the earthquake occurrence.The anomalies mainly distributed in the region of 22.5°~24.5°N and 99°~102°E.The time length from the peak point of I_RTL curve to the earthquake occurrence was 7 months.The above study showed that even the earthquakes location was near and the magnitude was close to each other,a big difference in seismic activity before the earthquakes may exist.Before the Jinggu M_S6.6,there was some difference in seismicity changes according to different beginning time of catalogues,but the distribution of anomalies and the time length from the peak point of I_RTL to the earthquake occurrence were uniform.So there was an important significance for exploring the relationship between the distribution of anomalies and the earthquake location,and the relationship between the time of the peak point of I_RTL and the earthquake occurrence time.     

Multiscale anisotropic diffusion for ringing artifact suppression in geophysical deconvolution data

Ringing artifact degradations always appear in the deconvolution of geophysical data. To address this problem, we propose a postprocessing approach to suppress ringing artifacts that uses a novel anisotropic diffusion based on a stationary wavelet transform (SWT) algorithm. In this paper, we discuss the ringing artifact suppression problem and analyze the characteristics of the deconvolution ringing artifact. The deconvolution data containing ringing artifacts are decomposed into different SWT subbands for analysis, and a new multiscale adaptive anisotropic filter is developed to suppress these degradations. Finally, we demonstrate the performance of the proposed method and describe the experiments in detail.     

吉林前郭地区一维P波速度模型研究

基于2013年10月31日至2014年8月25日吉林前郭地区发生的290次地震事件的P波到时数据,利用Kissling提出的VELEST程序,获取了前郭地区最小一维P波速度模型。并将新模型应用到前郭地震重定位和PTD定震源深度中,二者综合显示获取的最小一维P波速度模型对于提高地震定位的精度有一定的现实意义。     

山区谷底沉陷预测模型及其参数反演

我国西部地区地势复杂,沟壑纵横,地下开采极易导致边坡失稳,引发采动滑坡。在地下采动沉降与滑坡体挤压上升的叠加影响下,谷底区域地表沉降值明显小于类似地质采矿条件下的平原地区。为准确预测山区谷底区域地表沉降值,基于简支梁的弹性变形理论,并借助概率密度函数建立了山区谷底区域地表沉陷预计修正模型,明确模型参数物理意义及其取值方法。依据修正模型,以实测值和预测值之差平方和最小为原则构建适应值函数,基于模拟退火粒子群算法提出新的模型参数反演方法,借助MATLAB语言编制了相应的参数反演程序。最后将研究成果应用于山西某矿,得到谷底区域预测结果中误差为73 mm,与实测值基本一致,取得了较好的工程实践效果。     

基于逐减随机震源采样法的频率域二维黏滞声波方程全波形反演

全波形反演方法利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而,巨大的计算量是阻碍其发展的一个瓶颈问题。为此,研究者们提出了震源编码技术来减少计算量,但是此方法在模型更新过程中会引进随机串扰噪声,降低反演结果准确性。所以,在保证计算精度的情况下,本文提出了采用逐减随机震源采样的方法来高效计算全波形反演问题。笔者将此方法应用于频率域二维黏滞声波波动方程全波形反演,开始了在频率域进行随机震源采样类方法的研究,计算过程中共使用了依次增大的8个频率段;并应用Overthrust模型来验证此类随机震源采样法的正确性。实验结果表明：基于逐减随机震源采样法的反演结果与实际Overthrust模型的拟合误差为0.065 65,而应用基于全部震源的全波形反演方法得到的反演结果与实际Overthrust模型的拟合误差为0.064 64,两者差别不大;但计算用时由740 min减少到291.2 min,即计算效率提高了2.54倍。为了更好地确定方法的有效性,将其应用于Marmousi模型进行试算。模型试算结果表明：基于逐减随机震源和基于全部震源得到的反演结果与实际Marmousi模型的拟合误差分别为0.080 12和0.078 97,相差不大;但计算用时由1 218.9 min减少到274.4 min,计算效率提高了4.44倍。综上,在保证反演精度的情况下,基于逐减随机震源采样法的频率域全波形反演方法大大减少了计算量,具有不可替代的计算优势,并且没有引进随机串扰噪声。     

基于改进果蝇算法的非线性模型参数估计方法

在对基本果蝇优化算法的寻优流程进行深入分析的基础上,提出一种单方向搜索处理的改进果蝇优化算法（IFOA）。该方法可以对极值点为非零非负的非线性函数进行优化处理,将其应用于非线性模型参数估计。实例表明,IFOA方法在参数估计精度上优于线性近似法与非线性迭代方法;与以遗传算法为代表的智能搜索方法相比,其估计精度相当,并具有参数设置少、寻优过程简单、易于程序实现等优点。     

利用MODIS数据反演沙漠化地区地表温度的简化模式——以陕西北部地区为例

利用MODIS数据可见光波段、近红外波段和中红外波段,获得提取地表温度所需要的2个基本参数：地表比辐射率和大气透过率;然后对热红外波段b31、b32运用劈窗算法提取出陕北地区地表温度;并利用两种不同经验公式参数得出的结果与相应位置气象观测站观测的卫星过境时刻的地表气温进行比较。结果表明,劈窗算法简化模式能获得较准确的地表温度,符合陕北地区的实际地表状况,印证了应用该简化模式可以在大范围内快速实时监测地表温度。     

导线环垂直激发场高精度数值模拟

导线环垂直激发场是导线环在地下产生的激发场中垂直于地磁场的分量,在地面核磁共振找水方法中,垂直激发场在地下的分布直接决定地面核磁共振响应,探讨影响垂直激发场地下分布的诸因素进而研究这些因素对地面核磁共振响应的作用,对核磁共振找水正确野外施工具有一定的指导意义,有助于采集更高质量的野外数据.利用经典Chave算法对含有贝塞尔函数的强振荡积分核进行了数值积分试验,将其计算结果和对应的解析解进行了对比,试验结果显示该算法用于强震荡函数无穷积分的计算具有良好的收敛性、较高的精度和较快的计算速度,采用Chave算法对含有贝塞尔函数的强振荡函数进行积分是可行的.在此基础上,进一步计算并数值模拟了导线环垂直激发场及其相位,分别从地磁倾角、磁偏角、地下介质导电性、深度等方面探讨了影响垂直激发场及其相位分布的主要因素.