共查询到20条相似文献，搜索用时 528 毫秒
1.
慢度平方三角网格模型VTI介质二维QP波各向异性立体层析反演 总被引：1，自引：1，他引：0
基于慢度平方三角网格剖分模型下的各向异性射线扰动理论,建立了具有垂直对称轴的横向各向同性(VTI)介质中二维QP波立体层析所需的FRECHET导数矩阵,实现了慢度平方三角网格模型的VTI介质二维QP波各向异性立体层析.考虑到各向异性参数是影响地震波运动学特征的次一级因素,制定了先反演慢度再反演各向异性参数的反演策略.由于慢度平方三角网格模型的稀疏性,使得数据空间对模型空间的FRECHET导数矩阵规模被大幅压缩,在降低了计算成本的同时也很好的保证了界面位置、背景慢度、各向异性等各个参数的反演精度.理论数据算例证实了 FRECHET导数的正确性和反演策略的合理性. 相似文献
2.
针对层析成像模型矩形网格剖分存在的一些问题,提出了复杂结构三角网最小走时射线追踪层析成像方法。以Delaunay三角剖分的优化准则,根据模型点、线、面的几何结构关系,进行三角网格剖分;采用三角网波行面扩展法计算声波初至走时,追踪出任意接收点至激发点处的射线路径,射线路径包含了从接收点至源点的坐标、所在三角单元等信息,形成矩阵方程。应用稳定性较好的联合迭代重构技术求解矩阵方程,得到模型的速度分布。数值模拟结果表明,三角网射线层析成像方法分辨率高,成像结果更接近实际结构形态,有效地解决了矩形网剖分对复杂区域网格参数化灵活性差,速度界面描述精度低等问题。 相似文献
3.
为克服矩形网速度间断面及边界描述精度差的缺点,本文提出了一种基于波前最小走时单元的三角网射线追踪算法.该算法用三角单元剖分速度模型,把与波前面相邻、处于振动状态的三角单元称为波前单元,从源点所在单元开始构造由波前单元组成的波前单元域;1)选域中最小走时三角单元为当前单元;2)它先与相邻三角单元振动传递,再进行相邻单元各自内部节点的传递、并将处于振动状态的相邻单元扩展为波前单元;3)然后确认当前单元的振动状态、决定其在波前单元域中的去留,以此三步循环进行波前单元域的扩展及节点的最小走时和次级源搜索,直至波前单元域在剖分区域的消亡;最后利用各节点最小走时和次级源,通过最小走时搜索,拾取从接收点到源点的射线路径.数值模拟表明,该算法的模型剖分灵活性强,对速度间断面的描述准确,追踪结果精度高. 相似文献
4.
针对矩形网格模型的声波层析成像检测结构混凝土存在的主要问题:1)矩形网格无法充分模拟结构混凝土的速度间断面、2)射线追踪面临穿到凹多边形检测区域之外或穿越内部预留空洞的风险,本文研究了三角网声波层析成像检测技术.依据Delaunay三角剖分的优化准则,结合波前法对模型进行三角剖分;采用"基于波前最小走时单元的三角网射线追踪全局算法"计算声波初至时、追踪射线路径;应用代数重构技术求解层析成像反演控制方程,得到模型的速度分布.数值模拟及应用实例表明,三角网声波层析成像,模型参数化灵活、速度间断面描述准确,射线追踪可靠性强、精度高,反演成像分辨率高、更接近实际结构形态,适应具有复杂几何边界和内部结构的结构混凝土的质量检测. 相似文献
5.
田继枫 《地震地磁观测与研究》2017,38(6):41-47
随着计算机技术的飞速发展,大地电磁测深数据的各种反演方法得到空前发展。对于反演过程中的网格剖分问题,前人研究有限。本文就此进行研究,确定影响网格剖分方式的因素,找到一种合适的网格剖分方法。为此,设计一个复杂模型,对正演模拟得到的数据加入随机噪声作为反演的输入数据,选择不同的网格剖分方式,使用DASOCC与NLCG反演方法,以相同参数分别进行反演,对所得结果进行对比分析。结果显示,在设定的模型及参数条件下,针对DASOCC反演方法,纵向及横向网格均不加密,反演效果比较理想,即横向采用一个测点对应一个网格的剖分方式,纵向采用表层加密,100 m以下后一网格为前一网格厚度1.1倍的剖分方式;针对NLCG反演方法,纵向与DASOCC反演方法相同,采用不加密网格的剖分方式得到的结果较好,而横向则建议采用一个测点对应一个网格,并在2个测点之间插入一个网格的剖分方式。 相似文献
6.
7.
在复杂山地和复杂海底条件下,地表和海底的剧烈起伏对地震波数值模拟提出了更高的要求.常规有限差分法采用矩形网格对模型进行网格剖分,由于矩形网格自身的限制,起伏地表或起伏海底只能由一系列阶梯状折线代替,从而引起人为虚假绕射波.此外,在模拟液-固界面的反射波时,如果界面与网格线不一致,则需要更密的网格才能得到精确的结果.为了解决上述问题,本文将自适应网格生成技术引入到起伏海底速度模型的网格剖分中,采用高阶仿真型有限差分法(MFD)对曲线坐标下的声波方程波进行了数值模拟.利用自适应网格生成技术对速度模型进行网格剖分不仅可以准确地描述模型边界,而且可以有效消除虚假绕射波.高阶仿真型有限差分法可以有效压制频散提高计算精度.模型试算结果表明,本文方法对复杂海底模型具有很好的适应性. 相似文献
8.
依场地类别进行了强震记录分组,对模型参数的变化规律进行了统计分析.在模型随机参数向量满足独立性假设的前提下,得到了地震动随机函数模型的联合概率密度函数.引入数论选点方法对地震动随机函数模型的概率空间进行剖分,可以较少的样本点描述概率空间.以所选模型参数代表点代入地震动随机函数模型,即可以得到地震动时程样本集合.在集合层次上对比了模型预测地震动与真实记录的差异,两者在均值谱和标准差谱层次上均吻合较好,证实了模型预测结果的合理性. 相似文献
9.
在许多地震反演和偏移成像方法中,都要涉及到射线路径和旅行时的计算.本文将波前面三角形网格剖分和三维波前重建法射线追踪技术结合使用,实现了射线路径和旅行时的准确快速计算.三维波前重建法射线追踪过程中可以保证稳定合理的射线密度,克服了常规射线追踪方法存在阴影区的问题.波前面三角形网格剖分在描述和拆分波前面时更加准确有效,而且不需太多的网格数目,从而提高了射线追踪的精度和效率.该方法在三维复杂构造成像方面有独特的优势,目前在实际的Kirchhoff 偏移中的已经有相关应用. 相似文献
10.
11.
12.
13.
An unstructured mesh tidal model of the west coast of Britain, covering the Celtic Sea and Irish Sea is used to compare tidal distributions computed with finite element (FE) and finite volume (FV) models. Both models cover an identical region, use the same mesh, and have topography and tidal boundary forcing from a finite difference model that can reproduce the tides in the region. By this means, solutions from both models can be compared without any bias towards one model or another. Two-dimensional calculations show that for a given friction coefficient, there is more damping in the FV model than the FE model. As bottom friction coefficient is reduced, the two models show comparable changes in tidal distributions. In terms of mesh resolution, calculations show that for the M2 tide, the mesh is sufficiently fine to yield an accurate solution over the whole domain. However, in terms of higher harmonics of the tide, in particular the M6 component, its small-scale variability in near-shore regions which is comparable to the mesh of the model, suggests that the mesh resolution is insufficient in the near-coastal regions. Even with a finer mesh in these areas, without detailed bottom topography and a spatial varying friction depending on bed types and bed forms, which is not available, model skill would probably not be improved. In addition in the near-shore region, as shown in the literature, the solution is sensitive to the form of the wetting/drying algorithm used in the model. Calculations with a 3D version of the FV model show that for a given value of k, damping is reduced compared to the 2D version due to the differences in bed stress formulation, with the 3D model yielding an accurate tidal distribution over the region. 相似文献
14.
实现大地电磁法快速三维反演的关键在于找到快速计算灵敏度矩阵的方案. 本文在对大地电磁三维张量阻抗表达式进行深入分析的基础上，获得了三维快速松弛反演算法的灵敏度函数表达式，解决了三维快速松弛反演的核心问题；为了减小反演解的多解性，针对三维问题定义了最小构造函数，实现了求最小构造的三维快速松弛反演算法. 对二维棱柱体、三维棱柱体理论模型的大地电磁合成数据进行了反演试算，反演结果与理论模型相吻合. 对日本Kayabe地区实测资料进行了XY模式反演，反演得到的地电模型较好地反映了地热资源的分布. 反演试算结果表明大地电磁三维快速松弛反演具有计算速度快、结果稳定可靠等特点；使三维快速松弛反演在普通微机上得以实现. 相似文献
15.
We tested the usefulness of acquiring multicomponent GPR data to detect cracks in a historical building, and to monitor their dynamics, caused by a slowly and irregularly moving landslide. We used 2 GHz bipolar antennas in a configuration that allows for acquiring the in-line and cross-line electric field components with x- and y-directed antennas. The 2 × 2 data matrix was collected on a floor in the building along transects at four different times over a period of one year. The data were processed with a standard 2D scalar algorithm and with the latest 3D single component vector algorithm that corrects for antenna effects. We have implemented a 3D single component vector migration algorithm in a 2.5D sense to produce 2D slices of a 3D vector migration image by applying the algorithm on line data. This procedure allows for migrating single component line data taking into account all vector effects as well as three-dimensional wave propagation. We show that the 2.5D vector migration images have a much better in-line resolution than the migration images obtained by applying a standard 2D scalar migration algorithm.The GPR profiles agree with the a priori information about the structure of the floor. In particular, we detected two different types of anomalies, only a few of which can be due to utilities and to metallic mesh. Some shallower anomalies agree well with the cracks visible on the tiling, suggesting that some cracks can be directly detected using GPR. Visually there were no changes in the cracks on the floor and no clear changes in the GPR data could be attributed to possible subsurface changes in the cracks. The variations in the GPR images seemed primarily caused by changes in the coupling of the antennas with the investigated structure (floor) depending on the season when the measurements were made. For this reason the monitoring aspect of the survey is not successful. 相似文献
16.
Fei Li Tao Xu Minghui Zhang Zhenbo Wu Chenglong Wu Zhongjie Zhang Jiwen Teng 《地震科学（英文版）》2014,27(2):127-136
Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks, which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block, and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization, which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models. The damped least squares method is employed in seismic traveltime inversion, which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model. 相似文献
17.
High‐resolution data obtained from airborne remote sensing is increasing opportunities for representation of small‐scale structural elements (e.g. walls, buildings) in complex floodplain systems using two‐dimensional (2D) models of flood inundation. At the same time, 2D inundation models have been developed and shown to provide good predictions of flood inundation extent, with respect to both full solution of the depth‐averaged Navier–Stokes equations and simplified diffusion‐wave models. However, these models have yet to be applied extensively to urban areas. This paper applies a 2D raster‐based diffusion‐wave model to determine patterns of fluvial flood inundation in urban areas using high‐resolution topographic data and explores the effects of spatial resolution upon estimated inundation extent and flow routing process. Model response shows that even relatively small changes in model resolution have considerable effects on the predicted inundation extent and the timing of flood inundation. Timing sensitivity would be expected, given the relatively poor representation of inertial processes in a diffusion‐wave model. Sensitivity to inundation extent is more surprising, but is associated with: (1) the smoothing effect of mesh coarsening upon input topographical data; (2) poorer representation of both cell blockage and surface routing processes as the mesh is coarsened, where the flow routing is especially complex; and (3) the effects of (1) and (2) upon water levels and velocities, which in turn determine which parts of the floodplain the flow can actually travel to. It is shown that the combined effects of wetting and roughness parameters can compensate in part for a coarser mesh resolution. However, the coarser the resolution, the poorer the ability to control the inundation process, as these parameters not only affect the speed, but also the direction of wetting. Thus, high‐resolution data will need to be coupled to a more sophisticated representation of the inundation process in order to obtain effective predictions of flood inundation extent. This is explored in a companion paper. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
18.
各向异性介质模型电性结构复杂,如何进行合理的网格剖分成为获得高精度正演结果的关键,为此本文开展时间域航空电磁各向异性大地三维自适应有限元正演算法研究.通过结合非结构时间域有限元算法和自适应网格优化技术,实现各向异性介质条件下三维时间域航空电磁自适应正演.考虑到时间域航空电磁响应随时间的衰减特性,为了综合评价不同时刻的后验误差,本文将时间作为加权因子,调整各个时刻后验误差的相对权重,进而实现对浅部和深部网格的同步优化.通过与一维解析结果进行对比验证了本文算法的可靠性.数值实验结果显示电导率各向异性对自适应网格影响严重,其最大主轴电导率的数值及其分布特征直接决定了网格加密效果.此外,各向异性对时间域航空电磁三分量响应的分布形态和异常幅值也会产生严重影响,利用全域视电阻率极性图,可以很好地识别各向异性主轴方向. 相似文献
19.
本文基于非结构网格实现了海洋可控源电磁法三维有限元正演模拟.该算法采用完全非结构网格剖分,可以模拟任意起伏地形和复杂地电模型.为了避免场源的奇异性,采用一次场/二次场分解算法,一次场由基于Schelkunoff势函数的一维解析公式得到.为了提高算法的精度和效率,采用对测点附近单元和异常体区域进行体积约束加密的方法,实现了非结构网格的局部加密.一、二维模型计算和分析表明,本文采用的局部加密方法能够明显地改善算法的精度,最大相对误差基本在1%以内.对三维模型计算及对比分析,说明了该算法对三维可控源电磁正演的实用性.复杂海底地形模型的正演模拟表明,海底地形对电磁场的影响很大,在进行海洋可控源电磁资料解释时,地形的影响有必要考虑在内. 相似文献
20.
本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究. 相似文献