删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

正则化Kaczmarz算法在页岩纳米CT重构中的应用

本站小编 Free考研考试/2022-01-03

唐巍1,2,3,,
王彦飞1,2,3,,
1. 中国科学院地质与地球物理研究所, 中国科学院油气资源研究重点实验室, 北京 100029
2. 中国科学院地球科学研究院, 北京 100029
3. 中国科学院大学, 北京 100049

基金项目: 国家自然科学基金重大研究计划项目(91630202)与中国科学院先导科技专项(XDB10020100)资助


详细信息
作者简介: 唐巍, 男, 1992年生, 博士研究生, 研究方向为勘探地球物理.E-mail:tw9202@163.com
通讯作者: 王彦飞, 男, 博士, 研究员, 主要从事计算及勘探地球物理领域的研究工作.E-mail:yfwang@mail.iggcas.ac.cn
中图分类号: P631

收稿日期:2017-11-02
修回日期:2018-08-13
上线日期:2018-11-05



Application of regularizing Kaczmarz algorithm in shale Nano-CT reconstruction

TANG Wei1,2,3,,
WANG YanFei1,2,3,,
1. Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
2. Institutions of Earth Science, Chinese Academy of Sciences, Beijing 100029, China
3. University of Chinese Academy of Sciences, Beijing 100049, China


More Information
Corresponding author: WANG YanFei,E-mail:yfwang@mail.iggcas.ac.cn
MSC: P631

--> Received Date: 02 November 2017
Revised Date: 13 August 2018
Available Online: 05 November 2018


摘要
利用X射线计算机断层成像(CT)方法对页岩的扫描成像是一种无损的,对研究页岩微纳孔隙结构有重要意义的方法.传统的CT重构通常使用的是显式的滤波反投影(Filtered Back Projection,FBP)方法,该算法具有较快的成像速度,但常伴随有伪影或不稳定等情况.对于纳米CT而言,可以通过迭代优化的方法对投影数据进行成像,传统的迭代成像有收敛速度慢导致的计算时间长等缺点.Kaczmarz算法作为一种重要的代数重建技术(ART),由于其几何意义明显,操作容易等优点,在CT重构中起着重要的作用,我们可以通过块状迭代或随机迭代的方式对其收敛速度进行改进.对于所求解问题的不适定性,代数重建过程中需要引入正则化的技巧来改善解的稳定性.本文根据实际问题的需要,使用页岩数值模型,验证了正则化Kaczmarz方法的有效性,并对重庆漆辽龙马溪组页岩样品的实际数据进行了处理,得到了较好的效果.
X射线层析成像/
CT重建/
Kaczmarz算法/
正则化

X-Ray computerized tomography (CT) is a non-destructive and effective method to study the size, shape, 3D pore structures and interconnections of pores in shale. The main methods of CT reconstruction are FBP (Filtered Back Projection), ART (Algebraic Reconstruction Technique) and methods based on statistic models. FBP is most widely used in industry. The results obtained by this method contains artifact which is difficult to remove. Kaczmarz method, which has explicit geometric meaning, is one kind of ART methods. This method project the solution to the hyperplanes determined by the rows of matrix in a fixed order. Obviously the convergence rate depends on the order of the matrix's row, and it is not so efficient. If we project the solution in each step to the center of the projection points from different hyperplanes, which is known as block projection, we can improve the efficiency of computing greatly. Because the coefficient matrix is highly sparse and the observations are not only incomplete but also contaminated, the problem is not well-posed in the sense of Hadamard. Regularization method is an effective method to solve ill-posed problem. Many forms of regularization methods can be used to solve this problem. We transform the regularized normal equation to the equivalent augmented regularized normal system of equations to solve large dimensionality ill-posed problem. At the same time, we can find some constrains when we deal with the real problem like the absorption of pixels around the sample should equal to zero, the absorption of pixels should not be negative and so on. If we add these kind of constrains, we will get a better result. Our proposed data-constrained regularized Kaczmarz algorithm can perform better than the FBP methods and other projection methods in the numerical test and practical data test, hence the new method is promising for the future practical data processing.
X-Ray computerized tomography/
CT reconstruction/
Phase retrieval/
Tikhonov regularization



PDF全文下载地址:

http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_14765
相关话题/中国科学院 北京 计算 数据 代数