首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 234 毫秒
1.
利用磁异常总梯度模确定磁源边界位置   总被引:14,自引:0,他引:14  
本文介绍了一种利用磁异常总梯度模确定磁源边界位置的方法。根据理论证明与实际模型相结合的方法,分析了磁异常总梯度模极大值与磁源边界位置的对应关系,从而,可利用微机自动搜索边界位置成图。该方法在区域磁异常解释中,可广泛用于划分构造单元、确定构造带的位置、区分不同岩性与地层分布。  相似文献   

2.
台阶是重力梯度带异常解释的主要物理模型。通过模型计算,证实台阶布格重力异常Δg的一阶水平方向导数V xz极大值位置与台阶上沿的地面投影位置相一致,而与倾角关系不大。因此,利用V xz极大值可以方便地确定台阶上沿的地面投影位置。根据V xz曲线极大值两侧的对称性可以判断出台阶的倾向。同时推导出应用1/2与3/4V xz极值点位置反演台阶模型参数的计算公式,并证明这个公式比应用1/2与1/4V xz极值点反演精度更高。  相似文献   

3.
目前,激电数据的二维反演主要针对高密度激电数据,而对于常规的激电测深则很少涉及.据此,针对常规激电测深二维反演所涉及到的一些技术难点进行探讨.在正演过程中,根据实测数据实现有限元网格自动剖分,并采用双网格系统,用相邻节点电导率双线性连续进行有限元模拟,并对模拟的误差进行校正.在反演过程中,利用互换原理和Broyden的方法计算偏导数矩阵;采用与奇异值分解算法精度相当的共轭梯度算法求解偏导数矩阵方程.最后通过对模拟数据和实测数据进行反演试算,验证本计算方法是可行的.  相似文献   

4.
在MapGIS中自动生成地球化学等值线时,用三角剖分方法能最大程度地保留原始信息,但其图面可视化效果差,后期修改整饰的工作量巨大。用Gird网格化数据生成等值线,又会出现元素浓集中心和极值点位置存在一定程度的偏差和位移的情况,对后期的地球化学特征分析和异常查证工作带来了不便。经过研究发现,把离散采样点分析数据中的极值点提取出来,按照一定的阀值给极值点建立缓冲区,通过点对区空间分析,把规则网格数据中位置距离极值点较近,可能影响到三角剖分网优化效果的畸形网格节点剔除;再把离散的极值点数据和剔除了畸形网格节点的规则网格数据合并,进行三角剖分等值线追踪,很好地弥补了三角剖分生成等值线可视化效果差和网格化数据生成等值线出现浓集中心偏移的不足。该方法无需借助其他软件工具,在MapGIS中即可实现极值无偏移的地球化学等值线生成,既兼顾了制图精度,又满足了可视化效果,有很强的实用性。  相似文献   

5.
地质体边缘深度在重、磁位场数据半定量解释中起着至关重要的作用。由于重、磁异常及其各阶导数均满足欧拉齐次方程,tilt-Euler法在边缘深度反演方面备受青睐。然而,当重、磁异常的总水平导数或者总梯度模等于0时,倾斜角的一阶导数无法计算,导致倾斜角不能满足欧拉方程,tilt-Euler法无法使用。为了解决此问题,本文基于正则化思想,对倾斜角的一阶导数进行修改,使得重、磁异常的总水平导数或者总梯度模等于0时,倾斜角的一阶导数依然可以计算,修改后的倾斜角导数依然满足欧拉方程,称改进的方法为rtilt-Euler法;同时利用识别精度更高的归一化总水平导数垂向导数(NVDR-THDR)边缘识别方法对反演结果进行约束,剔除偏离边缘位置的坏点。理论模型试验结果表明,改进后的方法消除了重、磁异常总水平导数或者总梯度模很小或者等于0时,倾斜角导数无法计算以及反演结果不稳定的问题。将该方法应用到澳大利亚奥林匹克坝氧化铁铜金矿床边缘深度反演中,反演结果显示氧化铁铜金矿床边缘深度主要集中在0~100 m和100~200 m这两个深度段内,与沉积物剖面显示的矿床边缘深度0~200 m相符,证明了该方法的有效性。  相似文献   

6.
马国庆  吴琪  熊盛青  李丽丽 《地球科学》2021,46(9):3365-3375
场源中心位置的计算是重磁数据反演的主要任务之一,现主要通过异常与场源位置之间的数学物理方程来估算地质体的位置.为了快速、准确获得地质体的位置信息,提出基于重磁梯度比值的深度学习技术实现场源位置的获取;其利用深度学习技术所建立的重磁梯度比值水平分布与地质体埋深、构造指数的关系,快速实现异常场源位置计算,且提出利用多个值的相互关系来更加准确、稳定地计算出地质体的信息.该方法可以计算复杂地质体的中心位置,且避免了以往线性方程反演方法需对结果进行筛选的复杂过程,对于存在剩磁的磁异常则采用解析信号的深度学习方法来进行位置反演.理论模型试验证明利用梯度比值的深度学习方法可以准确获得地质体的深度,且通过对比更多点的深度学习计算结果发现,采用多个不同比例极值点可以减弱噪声带来的干扰,从而得到更加准确的位置.最后将该方法应用于实测磁异常的反演工作,获得了地下磁性物体的中心位置,且计算结果与欧拉反褶积法相接近,因此该方法具有良好的实用性.   相似文献   

7.
为减少反演的多解性,多数情况下期望电阻率和极化率异常具有相同结构。笔者利用交叉梯度算法实现了联合反演。在计算交叉梯度项过程中,提出一种新的离散化网格边界处理方法,提高了交叉梯度项计算的速度和精度。试算结果表明三维电阻率/极化率联合反演能有效地提高同结构异常体的反演效果。  相似文献   

8.
位场数据异常存在的线性特征往往对应着地下断裂构造、不同岩性地质体的边界接触带或者其他具有一定密度或磁性差异的构造特征。对这些线性特征进行增强、提取并进行半定量地解释是重磁资料处理的主要内容。然而,位场数据中若混入噪声干扰,再利用总水平导数、总水平梯度倾斜角等方法进行高次求导运算会把噪声放大,导致提取的线性构造位置发生偏离甚至出现错误。为了获取较为准确的线性构造及边界位置,选取对噪声干扰不敏感的小波模极大值方法,将该方法用于模型试验和实际矿区数据处理中,都能较好地定位出异常体模型边界和矿体的投影边界,表明小波模极大值是一种有效的重磁异常线性特征增强与提取方法。  相似文献   

9.
传统方法对重力异常的总梯度作归一化计算时,通常是对总梯度的模作算术平均。笔者对归一化函数的分母分别作了几何平均、调和平均、修匀平均和中值平均,计算归一化函数。通过模型试验对比这四种不同的统计平均,其中几何平均重力归一化总梯度在进行反演计算时,抗噪能力更强,反演结果更加集中、准确。将该方法应用于山东招远金矿采空区,改进后的方法可以更加准确地确定采空区的水平位置和中心埋深。  相似文献   

10.
水平梯度法提取重磁源边界位置   总被引:11,自引:0,他引:11  
本文介绍了一种用水平梯度极大值确定密度和磁性体边界的方法,并由计算机自动搜索成图。可广泛用于区域重磁异常的构造解释之中。  相似文献   

11.
A modification of the Elliott grid for plotting ellipse shape data are described. The effects of strain on elliptical markers are easy to visualize when the data are plotted on the new grid, and this allows clear interpretation of displayed data. New graphical methods for manipulating distributions of elliptical markers are directly related to an existing numerical method. When a distribution on the grid is strained, all the points move along straight, parallel lines. An initial distribution in which all the points lie on a straight line is strained into a distribution with the points lying on a hyperbola. Such curves include the analogues of ‘theta curves’. If the points lie on a circle centred at the grid origin, they are strained so as to lie on an ellipse. These are the analogues of ‘onion curves’.  相似文献   

12.
Adopting a random function model {Z(u),u ε study areaA} and using the normal equations (kriging) for estimation amounts to assume that the study areaA is embedded within a infinite domain. At first glance, this assumption has no inherent limitations since all locations outsideA are of no interest and simply not considered. However, there is an interesting and practically important consequence that is reflected in the kriging weights assigned to data contiguously aligned along finite strings; the weights assigned to the end points of a string are large since the end points inform the infinite half-space beyond the string. These large weights are inappropriate when the finite string has been created by either stratigraphic/geological limits or a finite search neighborhood. This problem will be demonstrated with numerical examples and some partial solutions will be proposed.  相似文献   

13.
Y. Jie  G. Guangmeng 《Atmósfera》2013,26(4):473-477
This paper presents a preliminary analysis of thermal anomalies before the Baja California M7.2 earthquake occurred on April 4, 2010. The results show that surface air temperature of the Mexicali station, which is the closest one to the epicenter, reaches its maximum value on March 30, which is also the maximum value from January 1 to April 10, 2010. Weather stations far away from the epicenter (100 km away) show no other thermal anomaly. In particular, the National Centers for Environmental Prediction Final Operational Global Analysis (NCEP-FNL) air temperature data show a maximum temperature change between March 30 and 31 at the epicenter grid point, while the other surrounding 7199 grid points show medium to low temperature change. The microwave brightness temperature measured by the Advanced Microwave Scanning Radiometer on board the EOS Aqua satellite (AMSR-E) also shows an increase on March 31. The results of multiple thermal parameters analysis obey well the DTS criterions: deviation (D), time synchronism (T), and space consistency (S), so they can be considered as convincing seismic thermal anomalies. This analysis might contribute to make reliable earthquake predictions.  相似文献   

14.
The first step in a seismicity analysis usually consists of defining the seismogenic units, seismic zones or individual faults. The worldwide delimitation of these zones involves an enormous effort and is often rather subjective. Also, a complete recording of faults will not be available for a long time yet. The seismicity model presented in this paper therefore is not based on individually defined seismic zones but rather on the assumption that each point in a global 1/2° grid of coordinates represents a potential earthquake source. The corresponding seismogenic parameters are allocated to each of these points. The earthquake occurrence frequency, one of the most important parameters, is determined purely statistically by appropriately spreading out the positions of past occurrences. All the other significant seismicity characteristics, such as magnitude-frequency relations, maximum possible magnitude and attenuation laws including the dependence on focal depth are determined in a global 1/2° grid of co-ordinates. This method of interpreting seismicity data allows us to establish a transparent, sufficiently precise representation of seismic hazard which is ideally suited for computer-aided risk analyses.  相似文献   

15.
地震波走时计算的逆风差分算法   总被引:5,自引:0,他引:5  
快速精确地旅行时计算在地震资料的叠前偏移与层析反演中起重要作用。利用逆足差分格式求解程函方程,首先在波阵面上寻找全局极小点,然后计算该点周围各点的地震波初至走时,同时,采用局部算法以确保计算的稳定性。该算法对于任意复杂地质模型都具有较高的精度。  相似文献   

16.
Kriging in a finite domain   总被引:2,自引:0,他引:2  
Adopting a random function model {Z(u),u study areaA} and using the normal equations (kriging) for estimation amounts to assume that the study areaA is embedded within a infinite domain. At first glance, this assumption has no inherent limitations since all locations outsideA are of no interest and simply not considered. However, there is an interesting and practically important consequence that is reflected in the kriging weights assigned to data contiguously aligned along finite strings; the weights assigned to the end points of a string are large since the end points inform the infinite half-space beyond the string. These large weights are inappropriate when the finite string has been created by either stratigraphic/geological limits or a finite search neighborhood. This problem will be demonstrated with numerical examples and some partial solutions will be proposed.  相似文献   

17.
目前复杂地质条件滑坡抗滑桩工作状态下的受力变形研究仍不够充分。为了研究江顶崖堆积体滑坡抗滑桩的桩身响应,在抗滑桩内布置弦式钢筋应力计和埋入式应变计,对抗滑桩进行长期桩身应力-应变监测,基于监测数据建立计算模型计算抗滑桩全桩弯矩,用有限元软件PLAXIS 3D建立模型对抗滑桩进行研究。结果表明:监测到的单根钢筋最大应力为112.5 MPa,远小于设计值,弦式应变计最大应变为34 με;建立的单筋矩形截面受弯构件计算模型较为可靠,可用于类似工程中;有限元软件计算所得弯矩小于设计弯矩,两者变化趋势相同;抗滑桩目前处于安全状态,滑坡治理效果显著。  相似文献   

18.
有限元法是复杂介质地震模拟的有力工具,它能比较客观地反映地震波传播过程,比较细致地再现地震图像。但由于有限元法地震模型所需网格节点数多、数据规模大,使它的计算受到限制。采用对角化集中质量矩阵、刚度矩阵变二维为一维存贮,并通过使用有限差分吸收边界条件来处理地质模型的人为边界反射问题,从而使计算所需内存降低,计算时间减少,并能比较好地实现地震声波用有限元法计算,理论模型的计算结果验证了算法的正确性和可行性。  相似文献   

19.
Variant是VB2.0及其以后版本的缺省数据类型,它是VCB中的一种特殊的数据类型可以包含数值,字符串或日期数据,还可以包括含特殊数值Empty,Error及Null等。Variant变量能够存储所有系统定义类型的数据,如果把它们赋予Variant变是,则不必在这些数据的类型间进行转换,VisualBasic会自动根据上下文字成任何必要的转换。这样我们可以使用Variant来替换任何数据类型,从而更有适应性。虽然使用方便灵活,但是Variant比其它数据类型所占内存大,并且运行速度慢,应该有限制地使用。  相似文献   

20.
航空物探测量中测线间水平往往存在差异,表现为沿测线方向的条带。调平处理旨在减少或消除测线间的水平差,是航空物探资料处理中的关键,对航空物探资料处理和解释具有重要作用。笔者以实际的航磁资料处理为例,详细介绍了航空物探资料微调平处理的原理与实现。微调平处理中首先通过方向滤波从原始资料中提取噪声网格,并按飞行测线提取噪声网格中的数据并分离有用的地质信息,最终将分离后的测线水平误差从原始资料中去除得到调平后资料,以此实现航空物探资料的微调平处理。同时,针对资料处理中遇到的诸如滤波器选取、异常分离阈值确定、最终调平效果检验等实际问题进行讨论,这对理解微调平并应用于实际资料处理均具有意义。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号