首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 62 毫秒
1.
针对重力与地震联合反演存在的问题,结合已有的研究成果,本文研究实现了速度和密度随机分布共网格单元模型的建模技术,以适应密度和速度剧烈变化的复杂模型及联合反演的计算要求.重力正演利用了该网格的二度半体模型,并进一步改进了地震走时的二维射线追踪计算方法,以适用于速度随机分布的网格介质.结合改进的模拟退火算法,实现了这种共网格条件下的重力与地震资料的同步联合反演.模型试验证明了重力与地震联合反演可以准确确定复杂物性界面的密度和速度结构,适用于物性界面不完全一致和物性变化剧烈的复杂模型,并且联合反演结果要优于单独的重力反演.带先验信息约束下的实际资料的联合反演,进一步证明了该方法的适用性和效果,可提高反演精度并减少多解性.  相似文献   

2.
三维连续密度分布的重力计算及应用   总被引:2,自引:0,他引:2       下载免费PDF全文
楼海  王椿镛 《地震学报》1999,21(3):297-304
给出了基于矩形网格模型的三维连续密度分布的重力解析公式.这种密度模型与一类地震层析算法采用的速度模型有相同的参数化形式.因此,地震层析成像的速度模型可以很方便地通过适当的速度密度函数转换为密度模型,从而进行地震 重力资料的综合解释.一致性的速度密度模型还有利于地震重力资料的联合反演.试验表明,在相同计算精度下,变密度分布情形的线性密度算法比常密度算法要快.利用这个算法对大别山东段的三维地震层析结果进行了地震重力的综合解释.   相似文献   

3.
地球物理反演是探索地下结构的最佳途径之一.地震波可以穿透到地球深部进行直接采样,是探测地球深部的主要方法.重力是结构体密度分布与地表观测点之间距离的体积积分效应,重力异常随着源深度的增加衰减很快,其对浅部结构的灵敏度明显优于地震数据.地震和重力联合反演能够相互补充和约束,提高空间分辨率,使反演结果更加稳定可靠.本文首先介绍了联系地震和重力数据的速度-密度经验关系,随后分别介绍了重震联合反演的3种常用方法—顺序反演、同步反演和交叉梯度反演,简要阐述了各种方法在国内外的应用情况.分析认为顺序反演将两类数据分开独立进行计算,原理简单,易于操作实现.但是该方法依赖于先验模型和速度-密度经验关系,可能存在分辨率较低区域(如模型边界)的误差放大效应.同步反演采用将地震和重力数据放在同一方程组中同时反演的策略,减弱了单一数据先验模型对结果的影响,但两种数据的同时运用势必引入数据权重分配问题.交叉梯度寻求不同物理参数模型在结构上的相似性,对潜在的岩石物性关系做了最少的假设,一定程度上降低了反演的非唯一性,但强制性地匹配模型的结构不一定完全符合地下介质的物性分布.因此使用交叉梯度方法反演时应注意模型的推导需要遵循客观标准,以控制模型的结构相似性和数据拟合度.最后指出重震联合反演中的速度-密度经验关系和数据的权重分配仍是值得探究的问题.  相似文献   

4.
基于变密度模型的位场界面反演   总被引:2,自引:1,他引:2       下载免费PDF全文
在波数域中用重力反演莫霍面深度时通常假定壳幔密度差为一常数,但这只是一种近似的密度模型,本文采用了密度随深度呈指数变化的变密度模型来反演莫霍面深度,给出了利用指数密度模型在波数域中计算重力异常的正演公式及界面深度的反演公式.利用指数密度模型及重力资料反演了青藏高原莫霍面的深度,分析了莫霍面的特征.结果表明,青藏高原莫霍面呈现出边缘浅、中部深的特点,边缘变化快、梯度大,中间变化梯度趋缓.中心地带的羌塘地体莫霍面深度达74 km,向四周慢慢变浅至67 km左右,边缘地区突然变浅至50km左右.通过常密度模型、变密度模型及地震反演得到的莫霍面的比较,证实变密度模型更适合于莫霍面结构的反演.  相似文献   

5.
储层重力密度反演后验约束正则化方法   总被引:2,自引:1,他引:1       下载免费PDF全文
本文针对蒸汽辅助重力泄油(SAGD)生产中开发监测问题,发展了综合应用地震及重力数据反演储层密度的联合反演算法.通过测井数据建立纵波阻抗与密度的直接关系,并推导出这种关系下重力与纵波阻抗数据联合反演的计算方法,从而计算出蒸汽腔体密度分布规律.文中应用密度反演后验约束正则化方法,采用Tikhonov正则化模型,通过波阻抗数据作为约束进行联合反演,在算法上提高了稳定性,同时得到较高的反演精度.文中对SAGD生产中的理论模型进行了方法试算,并分析了算法的误差,最终应用于SAGD生产的实际数据中,通过最终反演结果分析,该方法取得了很好的应用效果.  相似文献   

6.
二度半长方体组合模型的重力模拟退火反演   总被引:6,自引:4,他引:6       下载免费PDF全文
针对重力异常反演的现状,采用矩形单元体网格剖分的二度半体组合模型,模拟任意的多层密度界面或形体,实现了重力正演.该建模方法适应性强,适合于实际工作中常见的密度横向变化的复杂模型.使用改进的全局寻优的快速模拟退火算法,对重力异常进行反演,结合这种灵活的密度建模方法,反演过程中只需要反演各矩形单元体的密度参数,即可同时得到地质体的界面或形态以及密度值分布.通过复杂界面和形体的模型试验,证明了方法的效果.在江苏下扬子地质构造复杂区的实际应用说明,在地震等先验信息约束下,该重力反演方法提高了反演精度并减少了多解性,可有效解决古生界泥盆系-志留系目的层分布等地质问题.  相似文献   

7.
利用福建中部和东南沿海地区的重力布格异常资料,通过地震与重力3维联合反演所分离的不同层面间的重力异常信息,采用基于单层位理论和球近似条件下的简便公式,解算并分析研究不同过层间的密度变化分布图像及其特征,结果表明:不国圈层的密度分布特征与地表地质构造密切相关。  相似文献   

8.
为深入理解研究区的构造运动规律,进一步探求地壳物质可能的运移模式,本文基于Crust1.0地壳结构模型和最新莫霍面深度模型对自由空气重力异常数据进行联合约束,采用三维重力反演技术得到了郯庐断裂带中段的三维密度结构。反演结果显示,研究区的背景场稳定,断裂构造发育部位表现为重力密度低值,特别是营潍断裂带正处于大范围的重力密度低值区并纵贯整个研究区,由此可推断重力密度低值区可指示断裂构造的存在。   相似文献   

9.
根据重力数据反演结果的数据类型不同,可将重力反演概括为两类:一类是应用线性或者非线性算法,直接计算成像空间的密度分布;第二类是计算得到一种所谓等效密度,用以勾画重力异常场源分布,称为概率成像方法.相对于第一类方法,第二种方法计算更为简单,反演过程也更为稳定.但在一些应用中,需要获取地下介质的实际密度值,进而通过岩性层析等来确定岩性分布,此时,等效密度反演方法表现出了其局限性.本文从概率密度方法出发,应用一种迭代计算,可完成将概率密度值到实际密度值的转化,同时,增加了反演结果的分辨率.本文应用两个理论模型证明了该方法的有效性.  相似文献   

10.
蒙古及周边地区重力异常和地壳不均匀体分布   总被引:4,自引:1,他引:3       下载免费PDF全文
基于全球EGM2008自由空气重力异常模型,本文计算了蒙古及周边地区的布格重力异常和AiryHeiskanen均衡重力异常.在此基础上,本文采用Crust 1.0地壳模型为参考,通过重力正演方法,对蒙古及周边地区不同深度地壳密度结构模型的重力异常进行了计算,并对得到的正演布格重力异常与实际重力异常进行了对比和分析.研究结果表明:蒙古西部杭爱山地区与阿尔泰山地区的构造变形差异性明显,现今均衡重力异常中杭爱山周边没有明显的均衡异常高值区,而阿尔泰山地区西南方向存在均衡重力异常高值分布,分析与新构造运动密切相关;Crust 1.0模型给出的壳幔横向密度不均匀体分布对于计算Moho面起伏引起的重力异常作用明显;Crust 1.0给出的地壳内界面变形可以反映深大活动断裂的深部构造变形.研究结果对于认识蒙古东西部构造特征差异,以及现今西部活动断裂的地球物理场特征具有参考意义,也可以为进一步应用Crust 1.0模型为参考开展三维密度结构反演提供一定帮助.  相似文献   

11.
深内部地球结构对内核平动振荡本征周期的影响   总被引:2,自引:2,他引:0       下载免费PDF全文
地球固态内核的平动振荡是地球的基本简正模之一,又称Slichter模,其本征周期大约为几个小时,与地球内部结构密切相关.为了研究影响内核平动振荡的本征周期与内部结构的依赖关系,本文利用球对称、非自转、弹性和各向同性地球模型(SNREI),通过自由振荡运动方程的数值积分,以地球模型PREM为基础,理论上系统研究了地球内部介质(包括密度、地震波速等)分布异常对Slichter模本征周期的影响.数值结果表明,Slichter模周期随着内外核边界(ICB)密度差的增加以类似于双曲线的特征显著减小,当ICB密度差从597 kg·m-3减小到200 kg·m-3时,周期增大66.44%,当ICB密度差从597 kg·m-3增大到1000 kg·m-3时,周期减小21.48%;Slichter模周期随着核幔边界(CMB)密度差的增大而缓慢增大;相对于PREM,地球模型1066A在ICB和CMB的密度差分别相差45.321%和1.132%,内部地震波速度和密度梯度也存在差异,但是,当密度差减小到1066A模型提供的数值时,得到的Slichter模周期与基于1066A获得的结果(4.599 h)非常接近,差异分别只有3.762%和0.037%;表明Slichter模本征周期与地球内部介质的精细结构关系不大,而对ICB的密度差非常敏感.内、外核P波波速分布异常对Slichter模周期的影响基本相当,当内核和外核P波波速均增加5%时,Slichter周期分别减小1.02%和1.69%,P波波速分别减小5%时,Slichter模周期分别增加1.27%和1.847%,内核S波波速分布异常比P波波速分布异常对Slichter模周期的影响小1个量级;与地核相比,地幔中的地震波速异常对Slichter模本征周期的影响小1~2个量级;表明地核中地震波速异常对Slichter模周期的影响很小,目前有关Slichter模周期理论计算的差异主要来自于所采用的地球模型中内核边界的密度差的差异,本文结果可以为Slichter模的研究、探测及其对地球深内部结构的约束提供理论依据.  相似文献   

12.
The World Ocean Atlas 1998 is used to determine the global field of the meridional density ratio R hy =T/S, where temperature and salinity changes T and S are evaluated along meridians, in and below the mixed layer. The focus of the analysis is the identification of regions where the R hy field matches the values R =2 sometimes suggested as the commonly perceived state of the ocean and R =1, the condition of density compensation. Results are presented through fields of the meridional Turner angle Tu hy =arctan(R hy ) and through histograms of Tu hy for the Pacific, Atlantic and Indian Oceans at the ocean surface and at 300 m depth. At the 300-m depth level, which in the subtropics is representative of conditions in the permanent thermocline, the most frequently encountered values of the meridional density ratio are R hy =3.2 in the North and South Pacific, R hy =2.0 in the South Atlantic and Indian and R hy =1.6 in the North Atlantic Ocean. Conditions in the mixed layer are more variable and show seasonal differences, but R hy =2.0 occurs prominently in all ocean regions during winter and in all regions but the Atlantic during summer. Summer values for the Atlantic Ocean are R hy =3.2 in the Northern Hemisphere and R hy =2.4 in the Southern Hemisphere. Detailed analysis of R hy across the Subtropical Front (STF) confirms the most frequently observed values but shows zonal variation along the front in some oceans. Nearly complete density compensation (R hy =1) in the mixed layer is encountered in the STF of the eastern North Pacific, the eastern South Pacific and the eastern Indian Ocean. The eastern Indian Ocean south of Australia is also the only region where complete density compensation in the STF occurs below the mixed layer.Responsible Editor: Neville Smith AcknowledgementWe thank Dan Rudnick for helpful comments and discussion during the preparation of this paper.  相似文献   

13.
We present new experimental results on the mechanisms through which steady two-dimensional density currents lead to the formation of a stratification in a closed basin. A motivation for this work is to test the underlying assumptions in a diffusive “filling box” model that describes the oceanic thermohaline circulation (Hughes, G.O. and Griffiths, R.W., A simple convective model of the global overturning circulation, including effects of entrainment into sinking regions, Ocean Modeling, 2005, submitted.). In particular, they hypothesized that a non-uniform upwelling velocity is due to weak along-slope entrainment in density currents associated with a large horizontal entrainment ratio of E eq ?~?0.1. We experimentally measure the relationship between the along-slope entrainment ratio, E, of a density current to the horizontal entrainment ratio, E eq, of an equivalent vertical plume. The along-slope entrainment ratios show the same quantitative decrease with slope as observed by Ellison and Turner (, 6, 423–448.), whereas the horizontal entrainment ratio E eq appears to asymptote to a value of E eq?=?0.08 at low slopes. Using the measured values of E eq we show that two-dimensional density currents drive circulations that are in good agreement with the two-dimensional filling box model of Baines and Turner (Baines, W.D. and Turner, J.S., Turbulent buoyant convection from a source in a confined region, J. Fluid. Mech., 1969, 37, 51–80.). We find that the vertical velocities of density fronts collapse onto their theoretical prediction that U =-2?2/3 B 1/3 E eq 2/3 (H/R) ζ, where U is the velocity, H the depth, B the buoyancy flux, R the basin width, E eq the horizontal entrainment ratio and?ζ?= z/H the dimensionless depth. The density profiles are well fitted with?Δ?= 2?1/3 B 2/3 E eq ?2/3 H -1 [ln(ζ )?+?τ ], where?τ?is the dimensionless time. Finally, we provide a simple example of a diffusive filling box model, where we show how the density stratification of the deep Caribbean waters (below 1850?m depth) can be described by a balance between a steady two-dimensional entraining density current and vertical diffusion in a triangular basin.  相似文献   

14.
本文描述了X 线CT 机的几个主要技术特性,并根据这些特性,计论如何去选择一台价格低、图像质量高的CT 扫描机。  相似文献   

15.
采用动力学方程,求解了定态情况下磁层中上行离子沿磁力线的分布函数,针对不同的K指数,分别对北半球极光带区起源的上行离子O、H和He在子午面内沿磁力线的密度分布及其特性进行了研究.结果表明:1.沿磁力线向外,上行离子密度在近地空间呈急剧下降趋势,在远地空间呈缓慢下降趋势;2.重离子或初始能量较小的离子,其密度沿磁力线向外下降较快;3.K指数越大,离子进入磁层的空间范围越大;4.离子的投掷角分布对密度分布的影响甚小;5.离子密度沿磁力线的下降程度随到X轴距离的增大而呈增大趋势;6.在典型参数条件下,求得上行离子O在等离子体片边界附近的密度为10-3-10-2个cm-3,这与观测结果相一致.  相似文献   

16.
地震层析地球内部密度分布横向不均匀研究   总被引:2,自引:1,他引:2  
本文使用Woodhouse和Dziewonski等分别用S波和P波层析分析提供的,以670km深度为分界的上地幔和下地幔地震波传播所反映出的横向不均匀的球谐函数系数,假定地幔内部密度与正常水平的偏离分别与P波和S波的速度Vp和Vs成比例(下地幔δρ=γδVp,上地幔δρ=σδVS2),并用一组新的球谐函数系数表示全地幔密度分布模型.计算结果表明:1.比例系数γ=0.199,σ=3.13×10-5时,该全地幔密度模型可以70%恢复Dziewonski提供的下地幔密度模型,50%恢复Woodhouse提供的上地幔密度模型.合成密度分布基本上反映了分别由P波和S波速度异常所提供的地幔密度分布特征;2.所有穿过太平洋中部的剖面均明显地显示出一个低密度异常从地表一直延伸到核-幔边界;3.所有大陆下部均存在一个高密度异常区,而且深入到1200-1300km深度处;4.由密度异常的正、负交替基本上将部分地区地幔分成为3层,其深度依次为:第1层25-1300km左右,第2...  相似文献   

17.
Abstract

The paper explores some of the many facets of the problem of the generation of magnetic fields in convective zones of declining vigor and/or thickness. The ultimate goal of such work is the explanation of the magnetic fields observed in A-stars. The present inquiry is restricted to kinematical dynamos, to show some of the many possibilities, depending on the assumed conditions of decline of the convection. The examples serve to illustrate in what quantitative detail it will be necessary to describe the convection in order to extract any firm conclusions concerning specific stars.

The first illustrative example treats the basic problem of diffusion from a layer of declining thickness. The second adds a buoyant rise to the field in the layer. The third treats plane dynamo waves in a region with declining eddy diffusivity, dynamo coefficient, and large-scale shear. The dynamo number may increase or decrease with declining convection, with an increase expected if the large-scale shear does not decline as rapidly as the eddy diffusivity. It is shown that one of the components of the field may increase without bound even in the case that the dynamo number declines to zero.  相似文献   

18.
A densely spaced gravity survey across the San andreas fault zone was conducted near Bear Valley, about 180 km south of San Francisco, along a cross-section where a detailed seismic reflection profile was previously made byMcEvilly (1981). WithFeng andMcEvilly's velocity structure (1983) of the fault zone at this cross-section as a constraint, the density structure of the fault zone is obtained through inversion of the gravity data by a method used byParker (1973) andOldenburg (1974). Although the resulting density picture cannot be unique, it is better constrained and contains more detailed information about the structure of the fault than was previously possible. The most striking feature of the resulting density structure is a deeply seated tongue of low-density material within the fault zone, probably representing a wedge of fault gouge between the two moving plates, which projects from the surface to the base of the seismogenic zone. From reasonable assumptions concerning the density of the solid grains and the state of saturation of the fault zone the average porosity of this low-density fault gouge is estimated as about 12%. Stress-induced cracks are not expected to create so much porosity under the pressures in the deep fault zone. Large-scaled removal of fault-zone material by hydrothermal alteration, dissolution, and subsequent fluid transport may have occurred to produce this pronounced density deficiency. In addition, a broad, funnel-shaped belt of low density appears about the upper part of the fault zone, which probably represents a belt of extensively shattered wall rocks.  相似文献   

19.
Numerical simulations of variable-density flow and solute transport have been conducted to investigate dense plume migration for various configurations of 2D fracture networks. For orthogonal fractures, simulations demonstrate that dispersive mixing in fractures with small aperture does not stabilize vertical plume migration in fractures with large aperture. Simulations in non-orthogonal 2D fracture networks indicate that convection cells form and that they overlap both the porous matrix and fractures. Thus, transport rates in convection cells depend on matrix and fracture flow properties. A series of simulations in statistically equivalent networks of fractures with irregular orientation show that the migration of a dense plume is highly sensitive to the geometry of the network. If fractures in a random network are connected equidistantly to the solute source, few equidistantly distributed fractures favor density-driven transport. On the other hand, numerous fractures have a stabilizing effect, especially if diffusive transport rates are high. A sensitivity analysis for a network with few equidistantly distributed fractures shows that low fracture aperture, low matrix permeability and high matrix porosity impede density-driven transport because these parameters reduce groundwater flow velocities in both the matrix and the fractures. Enhanced molecular diffusion slows down density-driven transport because it favors solute diffusion from the fractures into the low-permeability porous matrix where groundwater velocities are smaller. For the configurations tested, variable-density flow and solute transport are most sensitive to the permeability and porosity of the matrix, which are properties that can be determined more accurately than the geometry and hydraulic properties of the fracture network, which have a smaller impact on density-driven transport.  相似文献   

20.
Sedimentation and density are important factors influencing the growth of wetland plants. The aim of this study was to elucidate the role of plant density in acclimation to burial stress by investigating the growth, biomass allocation, and carbohydrate contents of Polygonum hydropiper L. var. flaccidum, one of the dominant species in the Dongting Lake wetland, China. Experimental treatments, conducted in a greenhouse, combined three plant densities (16, 144, and 400 plants m−2) with two burial depths (0 and 15 cm) in a factorial design. Greater burial depth and higher plant density had negative effects on biomass accumulation and plant height with some exceptions. Plant growth, although lower at higher densities in general, was unaffected by burial depth at medium and high densities, indicating that the negative effect of burial stress was ameliorated by increased plant density. Deeper burial led to higher stem mass fraction only in the low-density treatment. Belowground mass fraction increased significantly with increasing density in both burial treatments. Moreover, higher density and deeper burial led to higher soluble sugar content but lower starch content. These data indicate that higher density facilitates acclimation of P. hydropiper to burial stress through increased soluble sugar content.  相似文献   

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

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