共查询到20条相似文献，搜索用时 15 毫秒
1.
针对地震信号多子波分解与重构技术中匹配追踪算法能够根据地震信号自身特点进行自适应分解、但其计算量庞大的问题,笔者提出一种粒子群快速优化算法,用于快速搜索地震信号稀疏分解的最优匹配原子。即在迭代过程中,将搜索区域确定在高斯函数能量集中的部分,避免了搜索过程的"贪婪性",能有效降低稀疏分解复杂度。同时,在粒子群算法中引入了一种多项式变异算子,可以有效避免搜索最优解的过度集中。实验结果证明,此算法将匹配追踪的分解精度提高了67倍,更使计算效率提高了153倍。 相似文献
2.
Construction of predictive reservoir models invariably involves interpretation and interpolation between limited available data and adoption of imperfect modeling assumptions that introduce significant subjectivity and uncertainty into the modeling process. In particular, uncertainty in the geologic continuity model can significantly degrade the quality of fluid displacement patterns and predictive modeling outcomes. Here, we address a standing challenge in flow model calibration under uncertainty in geologic continuity by developing an adaptive sparse representation formulation for prior model identification (PMI) during model calibration. We develop a flow-data-driven sparsity-promoting inversion to discriminate against distinct prior geologic continuity models (e.g., variograms). Realizations of reservoir properties from each geologic continuity model are used to generate sparse geologic dictionaries that compactly represent models from each respective prior. For inversion initially the same number of elements from each prior dictionary is used to construct a diverse geologic dictionary that reflects a wide range of variability and uncertainty in the prior continuity. The inversion is formulated as a sparse reconstruction problem that inverts the flow data to identify and linearly combine the relevant elements from the large and diverse set of geologic dictionary elements to reconstruct the solution. We develop an adaptive sparse reconstruction algorithm in which, at every iteration, the contribution of each dictionary to the solution is monitored to replace irrelevant (insignificant) elements with more geologically relevant (significant) elements to improve the solution quality. Several numerical examples are used to illustrate the effectiveness of the proposed approach for identification of geologic continuity in practical model calibration problems where the uncertainty in the prior geologic continuity model can lead to biased inversion results and prediction. 相似文献
3.
伪谱法求解波动方程的常规算法,运算量大,运算时间长。笔者综合提高伪谱法计算精度和效率的方法,改进了常规计算流程里计算傅氏变换和反变换中需要多次使用计算的常量,把这些常量提到时间循环体外,从而避免很多重复运算,减少了运算次数,使运算效率提高了一倍以上,并结合目前隐蔽圈闭在油气勘探中的重要性,构建了隐蔽圈闭中一种特殊的构造,即透镜体模型,模拟了高速和低速透镜体模型中的井间波场特征,对比分析了2种模型的合成井间地震记录、时间切片及振幅特征,与理论分析的结果能够很好的相吻合。 相似文献
4.
To improve the resolution of crosshole electromagnetic tomography, high precision of forward modeling is necessary. A pseudo-spectral time domain (PSTD) forward modeling was used to simulate electromagnetic wave propagation between two boreholes. The PSTD algorithm is based on the finite difference time domain (FDTD) method and uses the fast Fourier transform (FFT) algorithm for spatial derivatives in Maxwell's equations. Besides having the strongpoint of the FDTD method, the calculation precision of the PSTD algorithm is higher than that of the FDTD method under the same calculation condition. The forward modeling using the PSTD method will play an important role in enhancing the resolution of crosshole electromagnetic tomography. 相似文献
5.
这里实现了在2.5维各向异性介质中地震波传播的数值模拟。首先从2.5维一阶速度-应力弹性各向异性波动方程出发,得出了该方程的拟谱法数值解法;然后通过对Marmousi速度模型进行数值模拟,表明了拟谱法的有效性;最后计算并分析了均匀和混合各向异性介质中波的传播快照和理论模拟的地震记录,进一步认识了波在各向异性介质中的传播规律。 相似文献
6.
Specific structural–textural and mineralogical–petrographic features of clastic rocks, which make up sparse thin layers within Permian and Triassic cherty sequences in Sikhote Alin and Japan, are discussed. The cherty sequences were retained in the Jurassic accretionary complex as fragments of sedimentary cover of an ancient oceanic basin. They were mainly formed by turbidity currents that originated on intrabasinal uplifts composed of tholeiitic basalts in the Permian. In the Triassic, the currents originated on uplifts consisting of ocean-island alkali basalts along with Permian and Lower Triassic cherty and cherty–clayey rocks. The compositional difference of provenances was apparently caused by tectonic events in the oceanic basin at the Permian/Triassic boundary. 相似文献
7.
井间电磁波层析成像中的高精度时间域正演计算 总被引：1，自引：0，他引：1
为提高井间电磁波层析成像的分辨率, 正演计算中需保持较高的计算精度.本文正演计算中采用时间域伪谱(PSTD) 法模拟井间电磁波的传播.该算法在时间域有限差分(FDTD) 法的基础上采用快速傅立叶变换(FFT) 计算麦克斯韦方程中的空间导数.该算法除了具备时间域有限差分法的优点外, 在计算条件完全相同的情况下, 计算精度明显高于时间域有限差分法.时间域伪谱法的正演计算为井间电磁波高分辨率层析成像奠定了重要基础. 相似文献
8.
大地电磁(MT)数值模拟中通常使用有限单元法,通过伽辽金(Galerkin)法将微分方程转化为与其等价的泛函形式,对泛函求取极值并在单元上定义插值基函数,得到节点上电磁场值的线性方程组,最终形成大型复对称稀疏矩阵。要达到较高的有限元计算精度,一般采用密集的网格或高次插值的方法,这样做大大的减慢了正演的速度。结合两者的优点利用三次插值和h-型自适应相结合的有限元法来实现MT的正演算法。首先从一个粗网格出发并利用三次插值,通过后验误差估计方法局部加密网格,在计算量较小的情况获得较高的计算精度。这种方法可以针对目标区域和介质分界面发生突变处进行网格加密,不需要全局加密网格。最后通过对国际标准模型COMMEMI-2D1的模拟,分别比较二次插值与三次插值的自适应网格数量和数值模拟结果,证明了三次插值自适应有限元算法的可行性。 相似文献
9.
反演问题一直是地球物理的重点研究问题之一。本文以长导线源半航空瞬变电磁为例,对时间域半航空数据进行一维反演研究。将自适应正则化(ARIA)反演方法和阻尼最小二乘(LM)反演方法相结合,提出一种改进的反演方法:自适应正则化-阻尼最小二乘(ARIA-LM)反演算法。本文选用三层理论模型对比阻尼最小二乘算法、自适应正则化反演算法与本文方法,并使用四层理论地层模型来验证算法的有效性和可靠性。结果表明:该算法在一定程度上对高阻分辨率增强,能有效提高反演的精度及稳定性。 相似文献
10.
Explicit dynamic relaxation is an efficient tool that has been used to solve problems involving highly non-linear differential equations. The key feature of this method is the ability to use explicit dynamic algorithms in solving static problems. Few attempts have been made to date to apply this technique in conventional geotechnical engineering. In this study, an algorithm that incorporates the application of a stiffness dependent time step scheme is proposed. The algorithm has been successfully used to solve 2D and 3D non-linear geotechnical engineering problems. To calibrate the developed algorithm, numerical simulations have been conducted for a strip and square footings supported by Mohr–Coulomb material. Performance of four different types of brick elements used in collapse load calculation is examined in terms of convergence speed and accuracy. In addition, the role of employing adaptive time steps in reducing the number of iterations needed for convergence is also evaluated. 相似文献
11.
John F. Pernet-Fisher Geoffrey H. Howarth Yang Liu Peter H. Barry Laura Carmody John W. Valley Robert J. Bodnar Zdislav V. Spetsius Lawrence A. Taylor 《Contributions to Mineralogy and Petrology》2014,167(3):1-17
The Komsomolskaya kimberlite is one of numerous (>1,000) kimberlite pipes that host eclogite xenoliths on the Siberian craton. Eclogite xenoliths from the adjacent Udachnaya kimberlite pipe have previously been geochemically well characterized; however, data from surrounding diamond-bearing kimberlite pipes from the center of the craton are relatively sparse. Here, we report major- and trace-element data, as well as oxygen isotope systematics, for mineral separates of diamondiferous eclogite xenoliths from the Komsomolskaya kimberlite, suggesting two distinct subgroups of a metamorphosed, subducted oceanic crustal protolith. Using almandine contents, this suite can be divided into two subgroups: group B1, with a high almandine component (>20 mol%) and group B2, with a low almandine component (<20 mol%). Reconstructed REE profiles for B1 eclogites overlap with typical oceanic basalts and lack distinct Eu anomalies. In addition, elevated oxygen isotope values, which are interpreted to reflect isotopic exchange with seawater at low temperatures (<350 °C), are consistent with an upper-oceanic crustal protolith. Reconstructed REE profiles for B2 eclogites are consistent with oceanic gabbros and display distinct Eu anomalies, suggesting a plagioclase-rich cumulate protolith. In contrast to B1, B2 eclogites do not display elevated oxygen isotope values, suggesting an origin deep within the crustal pile, where little-to-no interaction with hydrothermal fluids has occurred. Major-element systematics were reconstructed based on mineral modes; group B1 eclogites have higher MgO wt% and lower SiO2 wt%, with respect to typical oceanic basalts, reflecting a partial melting event during slab subduction. Calculated residues from batch partial melt modeling of a range of Precambrian basalts overlap with group B1 trace-element chemistry. When taken together with the respective partial melt trajectories, these melting events are clearly linked to the formation of Tonalite–Trondhjemite–Granodiorite (TTG) complexes. As a result, we propose that many, if not all, diamondiferous eclogite xenoliths from Komsomolskaya represent mantle ‘restites’ that preserve chemical signatures of Precambrian oceanic crust. 相似文献
12.
We consider the inverse problem of permeability estimation for two-phase porous-media flow. The novel approach is based on
regularization by zonation, where the geometry and size of the regions are chosen adaptively during the optimization procedure.
To achieve this, we have utilized level-set functions to represent the permeability. The available data are sparsely distributed
in space; hence, it is reasonable to confine the estimation to coarse-scale structures. The level-set approach is able to
alter the boundaries between regions of different permeability without strict restrictions on their shape; however, when the
data are sparse, a reasonable initial guess for the permeability is required. For this task, we use adaptive multiscale permeability
estimation, which has the potential of identifying main permeability variations. These are described by a piecewise constant
function, where the constant values are attained on rectangular zones. In the current work, we develop a level-set corrector
strategy, assuming adaptive multiscale permeability estimation as a predictor. 相似文献
13.
14.
We focus on the question of whether high phytoplankton production events observed in a United States Pacific Northwest estuary consist of estuarine species blooms fueled by oceanic nutrient input or reflect offshore oceanic blooms that have advected into the estuary. Our approach is to use certain phytoplankton species as indicators associated with water mass origin, either estuarine or oceanic, to help resolve this question in Willapa Bay, Washington. We used species analysis and primary production data from 10 selected dates in May–September of 1998 and 1999, representing periods of high through low productivity. Out of 108 phytoplankton species identified from Willapa Bay, nine were selected and tested as indicators of oceanic species, six as estuarine, and two as surf zone. Our test results demonstrated the oceanic and estuarine species to be satifactory indicators of source waters. The prevalence of these species indicators in our samples revealed that the highest primary production and the appearance ofPseudo-nitzschia spp. were associated with oceanic intrusions of phytoplankton biomass into Willapa Bay. While the largest blooms were oceanic in origin, numerous medium-sized production events were from either oceanic, surf zone, or estuarine sources, indicating a complex situation. 相似文献
15.
16.
基于压缩感知理论的缺失地震数据重构方法 总被引：1，自引：0，他引：1
压缩感知是一种新的理论,它打破常规尼奎斯特-香农采样定理的制约,利用信号的稀疏特性或可压缩特性,用较少的数据即可重构恢复完整的信号。建立了基于压缩感知理论的缺失地震数据重构模型。首先在与稀疏变换不相关的测量矩阵基础上引入一种约束矩阵,使地震数据的缺失满足或接近高斯随机分布;随机缺失的地震数据变换到稀疏域会产生很多与有效信号不相干的随机噪声,接着通过一种新的自适应阈值迭代算法可以很好地消除稀疏系数中的随机噪声干扰,经过逆稀疏变换即得到重构后的地震数据。Marmousi 2模型测试及实际地震资料处理均验证了该方法的可行性和有效性。重构缺失地震数据取得了较好的效果。 相似文献
17.
在三维电阻率的正演计算中往往涉及到快速、准确求解大型线性方程纽Ax=b的问题。通过采用有限差分法来构造出求解点电源三维地电场的大型稀疏对称线性方程组。并引入Lanczos迭代技术，构造出三对角阵方程组，然后采用正交分解法进行求解，它是Krylov子空间方法中的一种。与传统迭代算法相比，它占用内存少，收敛速度快且稳定。针对大型稀疏矩阵及MATLAB语言的特点，采用简单记录矩阵的非零元素值及其所在行、列值的方法存储大型稀疏矩阵，可大大节省机器内存，提高运算速度。理论分析和计算实例显示，此算法是地电三维正演计算的有效方法，为下一步的反演计算打好基础。 相似文献
18.
19.
在中国区域大地构造研究中,对洋岛-海山/洋底高原的识别尚未引起足够重视.为深入研究中国大陆洋板块构造,系统回顾了洋岛-海山/洋底高原的基本概念、基本特征和增生造山过程.洋岛-海山/洋底高原是在海底扩张、大洋壳演化过程中由于地幔热点/柱作用形成的有异常厚度洋壳的区域,是大洋岩石圈的重要组成部分.洋岛-海山/洋底高原在垂向上具有典型的二元结构,下部以镁铁质、超镁铁质岩石为主,上部以碳酸盐岩建造为主.现今大洋盆地中大面积分布着正在演化中和正在俯冲的洋岛-海山,根据比较大地构造学原理,古洋岛-海山的存在指示古大洋盆地的存在,是研究造山带的重要载体.认为地史时期大洋盆地中有相当数量的洋岛、海山,在俯冲增生碰撞造山过程中保留下来的古洋岛-海山残块以构造岩片（块）形式夹持在俯冲增生杂岩中,随大洋盆地关闭;其作为缝合带的重要组成部分,是识别对接带的重要判别依据之一. 相似文献
20.
Noise Removal from Duplicate Paleoceanographic Time-Series: The Use of Adaptive Filtering Techniques
Martin H. Trauth 《Mathematical Geology》1998,30(5):557-574
This study demonstrates that adaptive filters can be used successfully to remove noise from duplicate paleoceanographic time-series. Conventional methods for noise canceling such as fixed filters cannot be applied to paleoceanographic time-series if optimal filtering is to be achieved, because the signal-to-noise ratio is unknown and varies with time. In contrast, an adaptive filter automatically extracts information without any prior initialization of the filter parameters. Two basic adaptive filtering methods, the gradient-based stochastic least-mean-squares (LMS) algorithm and the recursive least-squares (RLS) algorithm have been modified for paleoceanographic applications. The RLS algorithm can be used for noise removal from duplicate records corrupted by stationary noise, for example, carbonate measurements, species counts, or density data. The RLS filter performance is characterized by high accuracy and fast rate of convergence. The modified LMS algorithm out-performs the RLS procedure in a nonstationary environment (e.g., stable isotope records) but at the price of a slower rate of convergence and a reduced accuracy in the final estimate. The application of both algorithms is demonstrated by means of carbonate and stable isotope data. 相似文献