返回
Convolutional equivalent layer for magnetic data processing
DOI:10.1190/GEO2021-0599.1.png)
摘要
En 中文
We have developed a fast equivalent-layer method for processing large-scale magnetic data. We demonstrate that the sensitivity matrix associated with an equivalent layer of dipoles can be arranged to a block-Toeplitz Toeplitz-block (BTTB) structure for the case in which the observations and dipoles are aligned on a horizontal and regularly spaced grid. The product of a BTTB matrix and an arbitrary vector represents a discrete convolution and can be efficiently computed via 2D fast Fourier transform. In this case, the matrix-vector product uses only the elements forming the first column of the BTTB matrix and saves computational time and memory. Our convolutional equivalent-layer method uses this approach to compute the matrix-vector products in the iterative conjugate gradient algorithm with the purpose of estimating the physical-property distribution over the equivalent layer for large data sets. Synthetic tests of total-field anomaly data indicate a decrease in floating-point operations and computation runtime of our method compared to the classical approach of solving the least-squares normal equations via Cholesky decomposition. Faster results are obtained for millions of data, showing drastic decreases in computer memory usage and runtime and allowing for performances of magnetic data processing of large data sets on regular desktop computers. Our results also indicate that, compared to the classical Fourier approach, the magnetic data processing with our method requires a similar computation time but produces significantly smaller border effects without using any padding scheme and also is more robust to deal with data on irregularly spaced points or irregularly observation surfaces. A test with irregularly spaced field data over the Carajas Province, Brazil, confirms the efficiency of our method by estimating the physical-property distribution over the equivalent layer and computing the upward continuation.
Keyword:
INVERSE PROBLEM
CONTINUATION
期刊
IF:
3.2
论文数:
8.4K
被引数:
3.3W
机构
暂无机构信息
引用论文
Fast iterative equivalent-layer technique for gravity data processing: A method grounded on excess mass constraint
GEOPHYSICS
IF3.2
Gate Current Reduction and Improved DC/RF Characteristics in GaN-Based MOS-HEMTs Using Thermally Grown TiO2 as a Dielectric采用热生长TiO2作为介电材料的GaN基MOS-HEMT的栅极电流降低和DC/RF特性改善

