共查询到20条相似文献，搜索用时 140 毫秒
1.
为了克服现有数值计算方法在计算复杂断裂面接触滑移过程中效率低的问题，通过吸收边界元的建模思想，提出了新的数值计算方法。在该算法中，所有网格均位于岩块的边界，利用显示差分方法计算得到岩块之间的相互作用力以及岩块的位移量。利用“小球在抛物面上滑移”、“哑铃在不同倾角的斜面上”两个模拟试验对文中的算法进行验证，同时利用文中算法与有限元算法对直剪试验进行模拟。结果表明：算法能够准确描述物体之间的相互接触，准确计算得到接触块体之间的法向位移。算法能够准确判断出接触物体的“滑移状态”和“稳定状态”，摩擦力的数值计算结果与解析解的误差小于10?10。在计算断裂面接触和剪应力变化规律时，算法与有限元算法的计算结果一致，但文中计算方法的效率显著提高。 相似文献
2.
利用差分进化算法反求含水层参数 总被引：1，自引：0，他引：1
为更好地估算含水层参数,引入了差分进化算法。该算法是一种智能优化算法,其特点是借助现有近似解群体间的距离及方向指导未来的搜索行为,兼具有较好的全局寻优能力和较快的局部搜索能力。用此算法求解了均质各向同性及均质各向异性条件下Theis公式中的含水层参数,所用降深数据是在设定参数后利用解析解计算得到的。计算结果表明,该算法较传统方法计算精度高,不受人为因素干扰。 相似文献
3.
数值积分是基于节点计算的数值分析方法的关键技术之一。针对完全无网格法的计算要求，提出了一种能自动适应任意计算域上各种节点分布方式的数值积分算法。该算法能随计算点的位置不同，自动确定积分域及积分域内的求积点，计算出数值积分结果，给出了临空边坡滑移体和圆孔方板2个实例。结果表明，按文中算法编制的程序具有良好的边界再生能力和较高的计算精度，可作为基于节点数值分析方法的一个程序模块，用于实际计算。底层算法思想，使得程序可方便地转换为其他计算语言。 相似文献
4.
快速精确计算走时在层析成像反演中起重要作用，本文采用有限差分法来计算。首先通过变换，将程函方程化成一守恒形式；然后再差分离散，得到一阶和二阶精度的差分格式，并证明了差分格式条件稳定。模型试算表明：本文算法简明稳定，能快速高精度地计算地震走时，是一种值得注意和采用的走时计算方法。 相似文献
5.
6.
井间电磁波层析成像中的高精度时间域正演计算 总被引：1，自引：0，他引：1
为提高井间电磁波层析成像的分辨率, 正演计算中需保持较高的计算精度.本文正演计算中采用时间域伪谱(PSTD) 法模拟井间电磁波的传播.该算法在时间域有限差分(FDTD) 法的基础上采用快速傅立叶变换(FFT) 计算麦克斯韦方程中的空间导数.该算法除了具备时间域有限差分法的优点外, 在计算条件完全相同的情况下, 计算精度明显高于时间域有限差分法.时间域伪谱法的正演计算为井间电磁波高分辨率层析成像奠定了重要基础. 相似文献
7.
用容许应力法计算钢筋混凝土圆（管）桩的截面强度时，如果桩处于大偏心受压状态，其ｅ－α关系曲线上存在一个无穷间断点，在设计计算机程序时必须设法避开这个间断点。本文分析了该间断点的情况并提出一个利用一维变步长搜索法对偏心角α进行二次搜索的算法。用该算法编制的程序在ＣＯＭＰＡＱ－３８６型微机上通过，具有稳定性好，收敛速度快的特点。 相似文献
8.
9.
正交各向异性介质中地震记录的有限差分数值模拟 总被引：9，自引：3，他引：6
本文首先介绍了正交各向异性介质弹性波动方程的建立过程,然后讨论了高阶有限差分算法的实现原理、算法的稳定性和收敛性以及吸收边界条件的选取,并对正交各向异性介质作了正演模拟。通过波场快照、共炮点地震记录和ＶＳＰ模拟记录等正演分析,结果表明该算法在模拟各向异性介质中弹性波的传播是有效的 相似文献
10.
基于改进差分进化算法的阵列侧向测井反演 总被引：1，自引：0，他引：1
针对阵列侧向测井反演受给定初值影响的问题,将阵列侧向测井反演问题转化为无约束的非线性全局优化问题,提出了一种基于改进差分进化反演的新方法。在30次随机计算中改进差分进化算法与差分进化算法相比,其寻优成功率可以达到93.3%,平均进化代数为差分进化代数的1/8,其最优适应度也仅为差分进化最优适应度的一半。通过对阵列侧向模型进行反演计算,表明该算法满足电测井反演的需求。改进差分进化方法可有效地解决传统优化方法反演结果对初值的依赖性问题。 相似文献
11.
《International Geology Review》2012,54(12):1073-1086
The geologic structure of the Earth's upper crust can be revealed by modeling variation in seismic arrival times and in potential field measurements. We demonstrate a simple method for sequentially satisfying seismic traveltime and observed gravity residuals in an iterative 3-D inversion. The algorithm is portable to any seismic analysis method that uses a gridded representation of velocity structure. Our technique calculates the gravity anomaly resulting from a velocity model by converting to density with Gardner's rule. The residual between calculated and observed gravity is minimized by weighted adjustments to the model velocity-depth gradient where the gradient is steepest and where seismic coverage is least. The adjustments are scaled by the sign and magnitude of the gravity residuals, and a smoothing step is performed to minimize vertical streaking. The adjusted model is then used as a starting model in the next seismic traveltime iteration. The process is repeated until one velocity model can simultaneously satisfy both the gravity anomaly and seismic traveltime observations within acceptable misfits. We test our algorithm with data gathered in the Puget Lowland of Washington state, USA (Seismic Hazards Investigation in Puget Sound [SHIPS] experiment). We perform resolution tests with synthetic traveltime and gravity observations calculated with a checkerboard velocity model using the SHIPS experiment geometry, and show that the addition of gravity significantly enhances resolution. We calculate a new velocity model for the region using SHIPS traveltimes and observed gravity, and show examples where correlation between surface geology and modeled subsurface velocity structure is enhanced. 相似文献
12.
介绍了在倾斜均匀各向同性层状介质情况下,根据地震反射走时反演二维地层结构和地层速度的地震射线层析成像方法——走时反演。正演模型射线追踪是根据Fermat原理实现的,即求解满足该原理的非线性方程组得到射线与界面的交点,进而计算相应的走时。反演是先假设一初始模型,用最优化方法使射线追踪走时与观察走时的残差极小。最后计算了分辨矩阵和信息密度矩阵,以评价反演结果。对有噪情形也进行了反演。 相似文献
13.
投影解析速度建模 总被引：1，自引：0，他引：1
当今地球动力学的研究已逐步引起重视 ,尤其是高新探测技术的发展 ,将会大大加快其研究步伐。地球物理勘探作为重要的探测技术之一 ,已在该领域做出了重要贡献 ,比如能源探测、环境保护、防灾减灾工作中都需要物探技术。文中就地震勘探领域中地震波成像前沿技术———深度偏移的核心问题———速度模型的建立展开深入的研究。在分析现有方法缺陷的基础上 ,提出投影解析速度建模方法。该方法建立在两个重要事实依据之上 :(1)地震波旅行时误差是由描述地下反射界面的深度误差和描述地层的速度误差共同引起的 ;(2 )在反演中要保证零炮检距旅行时不变。在局部平面波前假设和拉东投影原理下 ,导出由旅行时误差计算界面深度误差和速度误差的解析公式。该方法可保证界面深度和层速度同步快速收敛。与现有方法相比 ,该解析法可以明显提高计算效率 ,且应用效果相当明显。 相似文献
14.
15.
井间地震走时波形层析成像方法 总被引：2，自引：1，他引：2
提出了井间地震走时波形顺序反演方法。该方法先利用井间地震走时反演得到速度模型的低频成分 ,然后用井间地震波形反演获得速度模型的高频成分。数值模型试验和实际应用结果表明 ,该方法反演稳健 ,提高了走时成像的分辨率 ,克服了波形成像易于陷入局部极小的缺陷 ,实现了快速高分辨率成像。 相似文献
16.
局部异常体井间地震走时层析成像研究 总被引：1，自引：0，他引：1
通过数值模拟试验 ,验证了利用井间地震走时层析成像方法技术可以分辨线性尺度为1个像元宽度的目标 ;讨论了表层低速带对基岩中异常体地震走时层析成像分辨率的影响 ;给出了高分辨率探测局部异常体井间地震走时层析成像方法技术 ,即最佳的观测系统、有效的数据采集技术、精细的数据预处理和优化的正反演算法。通过实际工程应用表明 ,该方法技术具有良好的效果 相似文献
17.
18.
The Ou Backbone Range strikes northwards through the central northeastern Japan arc and is bounded on both sides by the active reverse Uwandaira and Sen'ya faults. We have applied a traveltime inversion method (seismic tomography) with spatial velocity correlation to active and passive seismic data in order to investigate a three-dimensional (3-D) velocity structure. The data set contains 33,993 P- and 18,483 S-wave arrivals from 706 natural sources and 40 blasts, as well as 2803 P-wave traveltime data from 10 explosions detonated during the 1997 controlled source experiment. The traveltime inversion reveals a zone beneath the Ou Backbone Range in which P-wave velocities (VP) are approximately 6–8% lower than the average velocity at equivalent depths. The low VP and a low VP to S-wave velocity (VS) ratio (VP/VS) of about 1.65 suggest the presence of aqueous fluids in the middle crust. 相似文献
19.
地震波旅行时计算是射线追踪方法的核心问题,其精度直接影响着路径的精度和正演模拟的准确程度。在简要分析目前最典型的几种地震波旅行时计算方法(如有限差分法、线性插值法等)局限性的基础上,提出了一种基于非线性插值的地震波旅行时计算方法,并推导了其计算公式。通过模型的对比试算,不仅验证了所提出的旅行时非线性插值算法的正确性,而且还表明本算法计算旅行时的精度,比Asakawa[5]提出的线性插值算法的精度更高。 相似文献
20.
小波变换与地震奇异性属性 总被引：1，自引：0，他引：1
地震属性是由地震数据提取的有关地震波的几何形态、运动学特征、动力学特征和统计学特征等组成。这里主要讨论了一种新的地震属性(奇异性属性)计算方法。由于地震数据携带着大量的奇异性信息,因此,这一新的地震属性是在偏移地震数据的小波变换和奇异性分析的基础上形成的。小波变换使检测数据中的局部奇异程度成为可能。为能够估计出偏移地震数据每一采样点的局部奇异规律,从小波变换系数的角度来提取地震信号的奇异指数。这种新属性实行单道处理,不需要子波和速度信息,它给出了地下分层情况和断层位置。在地震数据奇异性属性的基础上,用小波变换来划分地层旋回,从而提高了小波变换在地震解释中的应用。 相似文献