共查询到20条相似文献，搜索用时 159 毫秒
1.
2.
由于离散运算,Radon变换方法中存在截断效应,传统的最小二乘约束反演并不能有效抑制这种效应。介绍了一种频域空间稀疏约束算法,在反演迭代过程中,根据前一次迭代的结果,通过贝叶斯原理将加权矩阵与前一次迭代的结果联系起来,得到新的加权矩阵;然后求解这个加权矩阵方程,得到频率域的稀疏解。比较了用最小二乘反演和压制截断效应的Radon域变换的效果,后者效果明显优于前者。模拟了一个5层的水平地层的地质模型,用这种压制截断效应的Radon变换进行了VSP资料的上下行波场的分离,得到了较好的效果。 相似文献
3.
4.
5.
多层地基非轴对称Biot固结的理论解 总被引：1，自引：1，他引：0
提出了一种新方法来推导多层地基非轴对称Biot固结问题的传递矩阵,相应的理论推导工作量较少。基于柱坐标系下非轴对称Biot固结的基本方程,通过引入中间变量,并对坐标 进行Fourier级数展开,对时间 和坐标 进行Laplace-Hankel变换,得到了6×6阶和2×2阶的两组常微分方程;然后,两组常微分方程进行关于 的Laplace变换和逆变换,得到了单层地基非轴对称Biot固结问题的传递矩阵;结合边界条件和连续条件,运用传递矩阵法得到多层地基非轴对称Biot固结问题在Laplace-Hankel变换域内的解,通过Laplace-Hankel逆变换得到了该问题物理域的真实解答。编制了计算程序,并进行了数值计算与分析。该方法具有计算速度较快的特点,方便工程应用。 相似文献
6.
利用有限正弦和有限余弦变换模拟地震波场:声波方程 总被引：4，自引：1，他引：3
研究了一种利用有限正弦和有限余弦变换解变密度声波方程的方法,它的基本思想有下列4点：(1)对声波方程两端进行关于一个或多个变量的有限正弦或有限余弦变换;（2）将待求波场的Fourier正弦或余弦级数表达式代入到经过变换后的波动方程之中,进而得到级数表达系数所满足的微分方程式;（3）分别对时间和未作变换的空间导数进行有限差分或其它形式的近似,得到变换系数所满足的离散（矩阵）方程;（4）将通过解离散（矩阵）方程而得到的变换系数代入到相应的反演公式中去,用求和的方式得到待求波场的数值解。与其它方法相比,利用有限正弦和有限余弦的优点在于可以无限制地扩展变换方向的计算区间和可以处理任意变化的速度和密度结构而不明显地增加计算量。 相似文献
7.
卷积完全匹配层（convolution perfectly matched layer,CPML）吸收边界是一种高效处理波动方程数值模拟中人工边界反射波的方法。本文基于传统的一阶系统CPML吸收边界条件推广并推导了新的二阶系统CPML边界条件（NCPML）。与常规二阶系统CPML边界条件不同,新边界条件推导的核心思想是在复数-频率域中忽略部分衰减因子空间变化特性,避免其在时间域产生复杂卷积算子,然后反变换至时间域得到基于CPML吸收条件的二阶标量波方程,并应用于二阶标量波方程数值模拟。均匀介质模型测试验证了NCPML吸收条件在内存使用上相对于常规二阶系统CPML与SPML（split PML）吸收条件更少。在对人工边界反射的吸收效果上,NCPML稍逊色于常规二阶系统CPML,但二者均相对于SPML优势明显。最后通过层状模型和Marmousi模型测试验证了NCPML的稳定性及其在效率上的优势。 相似文献
8.
我们建立了一个频率域的位场转换解释系统,由于可以利用褶积定理(空间域的褶积对应于频率域的乘积)的特点,在导出一系列位场转换的频率域滤波算子的基础上,不但可以实现各种单一的位场转换,而且可以很方便地进行综合的位场转换——将若干个频率域的滤波算子连乘(滤波器的串联),从而使这个解释系统具有多功能位场转换的特点,以适应解释的要求。此外,借助于快速富里叶变换的计算机算法,使这个系统在计算机上能省时而准确地实现计算。鉴于上述两方面的优越性,近几年来这个系统在重磁资料解释实践中,特别是在解释宽缓异常时,逐步地得到了推广使用。 相似文献
9.
10.
针对高光谱遥感图像波段多,各波段受云层影响程度不一致,云层本身不均匀等特点,提出了一种频域自适应同态滤波方法.该方法首先检测某个图像云层的范围和对应厚度,根据云层处于低频特点,利用变差函数确定滤波窗口,把云层按照空间邻域大小转换为傅立叶空间,根据云层厚薄调整高通滤波截止频率,然后反变换到空间域.该方法可解决频域空间同态滤波方法破坏无云区的信息,克服空间域滤波方法效果不佳,难以控制和解释的缺陷.结果表明,该方法在保持无云区信息不破坏前提下,较好地改善了不同厚薄云层的质量,计算机内存开销很小,对大数据处理速度较快. 相似文献
11.
The authors propose a semi‐analytical approach to studying wave propagation in multilayered poroviscoelastic grounds due to transient loads. The theoretical development is based on the exact stiffness matrix method for the Biot theory coupled with a matrix conditioning technique. It is developed in the wavenumber frequency domain after a Fourier transform on the surface space variables and the time variable. The usual methods yield a poorly conditioned numerical system. This is due in particular to the presence of mismatched exponential terms. In this article, increasing exponential terms are eliminated and only decreasing exponential terms remain. Consequently, the method can be applied to a large field of configurations without restriction concerning high frequencies, large Fourier transform parameters or large layer thicknesses. Validation and efficiency of the method are discussed. Effects of layering show that the layer impedance influence on solid and fluid displacements. Moreover, this approach can be of interest for the validation of numerical tools. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
12.
矫希国 《吉林大学学报(地球科学版)》1988,(1)
本文主要根据线性方程组的理论给出了求解投影矩阵的方法,再根据投影矩阵与广义逆矩阵之间的关系,进一步又给出了一类广义逆矩阵的求解方法,即满足Penvose’s第一条与第三条广义逆矩阵的一般表达式。 相似文献
13.
The fundamental solutions were obtained for step‐like point forces acting in three orthogonal directions and an instantaneous fluid point source in a fluid‐saturated, porous, infinite solid of transversely isotropic elasticity and permeability. After expressing the governing equations in the form of matrix in the Laplace space, we employed Kupradze's method together with the triple Fourier transforms. This method reduces the simultaneous partial differential equations with respect to three displacement components and a pore fluid pressure to a differential equation in terms of only one potential scalar function, which can be operationally solved in the transformed space. After the Laplace inversion of the potential, the residue theorem was applied to its Fourier inverse transform with respect to one of the transformation variables. The Fourier transforms with respect to two other variables were rewritten into the Hankel transforms. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
14.
多因子线性变换是处理多波段遥感数据的数学方法,应用该方法对于图象主要特征的分离和提取具有很好的效果.但是,遥感图象数据中的岩性信息光谱分布复杂多样,在信息强度和分布方面较弱,使用现有的变换很难将岩性信息提取出来.为了完成遥感岩性填图的任务,本文提出利用目标向量空间线性变换的方法,简化所要提取的目标信息的光谱结构,提取目标信息.根据研究区内岩性目标地物的图象灰度光谱分布选择空间域坐标轴的位置,使需要分离的岩性信息只在一个轴上的投影不为0,从而将岩性目标信息分离出来.该方法在新疆且末卡特里西地区提取铜矿带岩性取得了较好的结果. 相似文献
15.
16.
Isometric Logratio Transformations for Compositional Data Analysis 总被引：37，自引：0，他引：37
J. J. Egozcue V. Pawlowsky-Glahn G. Mateu-Figueras C. Barceló-Vidal 《Mathematical Geology》2003,35(3):279-300
Geometry in the simplex has been developed in the last 15 years mainly based on the contributions due to J. Aitchison. The main goal was to develop analytical tools for the statistical analysis of compositional data. Our present aim is to get a further insight into some aspects of this geometry in order to clarify the way for more complex statistical approaches. This is done by way of orthonormal bases, which allow for a straightforward handling of geometric elements in the simplex. The transformation into real coordinates preserves all metric properties and is thus called isometric logratio transformation (ilr). An important result is the decomposition of the simplex, as a vector space, into orthogonal subspaces associated with nonoverlapping subcompositions. This gives the key to join compositions with different parts into a single composition by using a balancing element. The relationship between ilr transformations and the centered-logratio (clr) and additive-logratio (alr) transformations is also studied. Exponential growth or decay of mass is used to illustrate compositional linear processes, parallelism and orthogonality in the simplex. 相似文献
17.
Spectral unmixing estimates the abundance of each endmember at every pixel of a hyperspectral image. Each material in traditional unmixing algorithms is represented through a constant spectral signature. However, endmember variability always exists due to environmental, atmospheric, and temporal conditions, which leads to poor accuracy of the estimated abundances. This paper proposes a new unmixing algorithm based on a new linear transformation called endmember orthonormal mapping (EOM) to overcome the aforementioned problem. The EOM transformation maps original spectral space to a new EOM space to reduce endmember variability. In the original spectral space, each material is represented by a set of spectra (endmember set) which is extracted using the automated endmember bundles (AEB) method. The EOM transforms each endmember set to a vector in the EOM space so that these vectors are orthonormal. On account of orthonormalized endmembers, the condition number of the mixing matrix in the EOM space reduces. Furthermore, we consider the noise term as an additional virtual endmember set mapped to a vector that is orthogonal to other endmembers. As a result, a promising unmixing accuracy is obtained through applying the least squares abundance estimation in the subspace orthogonal to noise. Experimental results of both synthetic and real hyperspectral images demonstrate that the proposed algorithms provide much enhanced performance compared with the state-of-the-art algorithms. 相似文献
18.
抛物Radon变换法(Parabolic Radon Transform)在地震资料处理中有广泛的应用。PRT可对不同频率的地震数据解耦处理,这一特点使得抛物Radon变换的计算效率比双曲Radon变换有数量级上的提高。在频率域求解时,需要对每一个频率成份求解同样大小的线性方程组。求解抛物Radon正变换的计算方法主要有Levinson递推法、共轭梯度法、Cholesky分解法和直接矩阵求逆法。最小平方抛物Radon正变换所形成的矩阵具有Toeplitz结构,可采用Levinson递推法进行计算。高分辨率抛物Radon正变换所形成矩阵的Toeplitz结构被破坏,一般采用共轭梯度法或Cholesky分解法进行求解。这里详细推导了复Toeplitz矩阵的Levinson递推算法,并分别对求解方程的四种方法进行了讨论,最后给出抛物Radon正变换求解的数值算例,并对所给出的四种方程求解方法的计算效率及计算精度进行了对比。 相似文献
19.
《矿物学报》2013,(Z1)
The interlayer clay-organic complex is an important clay-organic association in sedimentary environments. The interlayer region of clay minerals not only provides storage space for organic matter, but also possesses solid acid sites; and these acid sites were proposed to be involved with the transformation of organic matter into liquid and gaseous hydrocarbons. However, the effect of the organic matter storage in the interlayer space of clay minerals on the hydrocarbons generation has not been made clear. In this study, the interlayer complex of 12-aminolauric acid (ALA) and Na+-montmorillonite (Na+-Mt), labeled as ALAinter-Mt (Na), was synthesized to investigate the role of the interlayer space of montmorillonite in hydrocarbon generation. Simply mixed ALA-Mt complex [ALA-Mt (Na)] was also prepared for comparison. The pyrolysis of ALA and ALA-Mt complexes was studied using thermogravimetry coupled with Fourier transform infrared spectroscopy (TG-FTIR) and a high temperature-pressure system (a confined gold capsule-autoclave system). 相似文献
20.
广义位势理论为岩土材料的本构关系研究指出了更开阔的数学背景。基于广义位势理论可以构造出基于应力空间的模型,也可构造出基于应变空间的模型,利用两者之间的对称关系,把常规基于应力空间的模型直接转换到应变空间中,对本构模型自身无需做任何改动。基于广义位势理论也可以对已有的本构模型进行修正,比如,可基于广义位势理论对剑桥模型进行修正,这种修正可使模型利用三轴试验的部分成果,比单纯使用压缩试验的成果要好一些。另外,也初步探讨了模型的三维化问题。 相似文献