基于地质统计方法与DEM的地震灾情空间插值研究

历次破坏性地震的震害调查和强震观测资料显示, 地形地貌对地震灾害有着显著的影响. 地震发生后, 为了能够及时、 准确地为地震救灾指挥提供灾情分布信息, 该文借鉴了地质统计学方法, 利用灾情速报人员上报的地震现场离散点灾情短信对灾区进行灾情空间模拟的同时, 将数学高程模型(DEM)中所包含的高程、 坡度等地形地貌信息作为影响因素引入协克里金(Co-Kriging)插值; 并以汶川M_S8.0地震灾情短信数据为例, 分别对确定性插值、 地质统计学插值结果与有无考虑坡度因素的地质统计学插值结果进行了交叉检验. 结果表明, 考虑坡度影响因素的协克里金插值在合适的模型和参数下取得了最优的灾情模拟效果. 该方法为地震应急期间进行较高精度的灾情模拟提供了一种新的可行思路.     

位场曲化平的插值-迭代法

将起伏曲面B上的位场向下延拓至曲面最低点的平面A的插值-迭代法步骤是：1）将曲面B上的场值放置在水平面A上具有相同水平坐标的点上,作为A上的初值;2）用若干水平面切割B,从A的初值,用快速傅里叶变换法(FFT)向上延拓出这些平面的场值,用插值的方法从这些平面的场值计算曲面B的场值;3）根据B上的实测值与计算值的差值,对A上的值进行加权改正;4）重复步骤2）和3）,直到B上的差值小到可以忽略.这种插值-迭代法具有高的计算速度,比通常的FFT法延拓得更深,可以超过10倍点距.文中给出计算实例.     

三维地质空间插值重构方法的对比研究

钻孔数据是工程地质、水文地质、钻井测井和地震勘探等各项地质工作开展的基础资料.钻孔数据能够提供地质(地层)构造和矿产资源的分布信息,完整准确地表达复杂地质现象的边界条件及地质体内含的各种地质构造,直观地再现地质单元的空间展布及相互关系,最大限度地提高地质分析的直观性和准确性.     

基于克里格方法的相空间重构统计特征研究与应用

提高矿藏储量的预测精度是资源评价的难点之一.本文利用地质数据的相关性和变异性,建立了相空间重构的克里格方法,并结合限制性克里格法,给出了带有限制特征的计算理论,并设计出在特定铜矿中的储量计算流程,用阿舍勒铜矿体作为实例模拟了储量估计.数值计算结果表明,相空间重构克里格储量计算方法能够较为理想的预测储量.     

任意曲线上二维重磁位场转换的B样条函数法

重磁位场转换在重磁资料解释中是十分重要和必要的。随着物探仪器的改进和资料解释精度要求的提高,人们对位场转换方法的研究越来越深入。本文提出一种原理简单、计算精度高的新方法--B样条函数法,应用其插值、求导、求积的性质,以解决曲线上二维重磁位场转换的两类问题。一是Z_α分量向H_α分量的转换及向上延拓,二是水平和垂向导数的换算及磁位的计算。 样条函数为分段多项式函数。B样条函数是用σ函数逼近的一类样条函数,其节点     

基于储层骨架的多点地质统计学方法

建立弯曲、连续的河道形态是河流相储层建模的重点和难点.在分析了不同随机建模方法优缺点基础上,提出了一种针对河流储层新的建模方法,即基于储层骨架的多点地质统计学方法(SMPS).它综合利用了基于目标建模的优点以及多点地质统计学方法的长处.根据基于目标方法建模思路,预测河道骨架(即河道中线)分布,提出采用搜索窗预测河道骨架.随后,利用河道骨架对多点地质统计学中的数据事件选择进行约束,从而有效提高数据事件选择的合理性,达到更准确建立河道地质模型的目的.概念模型和实际储层建模对比研究表明,SMPS相比较于已有的序贯指示建模、Snesim和Simpat,确实能够更有效的建立起河道储层地质模型.这一新方法的提出对推动储层随机建模方法研究具有重要理论意义,对油田生产建立高精度储层地质模型也具有重要现实意义.     

基于Kriging方法的三维地质可视化系统构建

针对传统三维地质建模方法存在的诸多问题,对基于Kriging方法的三维地质体构建进行研究。选取某地震工程项目钻孔数据为主要研究对象,在Arc GIS平台上通过地统计分析,实现了复杂三维地质体的快速、准确建立。在此基础上,采用Arc Engine作为二次开发组件,完成该地质体三维展示系统的开发,实现三维展示、专题图出图、三维分析等功能。     

基于含砂率的相控地质统计学反演储层预测技术研究

陆相湖盆三角洲沉积储层横向变化快、纵向薄砂泥互层,储层预测难度大.针对该问题,本文以乌石油田群流沙港组二段中层序为背景,提出了一种基于含砂率的相控地质统计学反演储层预测技术.通过单井统计,建立含砂率与沉积微相的关系;再利用平面含砂率横向约束建立三维含砂率初始模型;进而运用贝叶斯理论将三维含砂率初始模型作为约束条件融入到...     

基于InSAR与GPS观测的汶川同震垂直形变场的获取

发震断层的形变是断层活动的重要参数之一,对认识断层性质、震源机理有重要作用。文中以逆冲性质为主的汶川地震为例,采用符合地表水平形变特征的Biharmonic样条插值对GPS水平形变矢量插值,然后再分解为EW和SN向分量。利用可靠的GPS观测值对InSAR参考点进行校正,统一两者的坐标系。通过对汶川地震视线向形变场剖面与GPS对比分析发现,断层上盘GPS与InSAR观测参考点相差9.93cm,而下盘则为-11.49cm。在此研究基础上,通过GPS水平形变场与InSAR视线向形变场联合解算,获取了汶川地震垂直连续形变场。结果表明,断层两侧垂直形变衰减较快,横跨断裂带形变量30cm的宽度不超过50km;沿发震断层附近垂直形变高值区分布不均匀,主要集中分布在发震断裂的汶川县城至都江堰段、茶坪—北川—南坝段和青川段。这3段各有特色,南段断层两侧垂直形变极不对称,主要以上盘剧烈抬升为主,最大抬升区域在映秀镇至连山坪一带,抬升量达到5.5m。中段表现为较强的反对称性,断层一侧抬升另一侧沉降。该段上盘最大抬升区域在茶坪东侧,抬升量为255cm,下盘最大沉降量在永庆,沉降量为-215cm。北端垂直形变相对较小,主要分布在青川北侧,呈对称分布,在发震断层最北端,最大抬升量为120cm。     

基于地震监测台网资料近实时插值计算震动图的初步研究

本文首先介绍了国内外有关震动图研究方面的动态。针对国内大部分地区的地震监测台网密度还比较低的现状,在已知震源位置和发震时刻的前提下,作者提出了利用实际地震台网观测资料近实时插值计算网格点处地震图的方法,在插值计算时充分考虑地震波在时间上的传播特征及在空间上的衰减规律。随后,对该方法的适用性进行了验证。以0．05°（约5Km）为尺度先对整个福建地区进行网格化划分,并以顺昌MIA．9级地震为例采用以上方法计算得到所有网格点处的地震时程。随后以1s为间隔,分别得到该时间间隔内的最大值等值线图,最后将所有的等值线图连续播放,得到了一种新形式的“ShakeMap”。作者还给出了顺昌地震最终的峰值地震动（PGV、PGA）分布,并分别对实际记录最大值及插值记录最大值随震中距的衰减关系进行了对比分析,结果充分证实了本文结果的可靠性。     

基于实测值分段线性内插模型的InSAR视线向同震位错分解-以青藏高原昆仑山MemＳsub8.1地震为例x,auto,auto,415px);}stylehf=httprztl.com business week launch

利用现场GPS定位的实测值，在尝试过最小二乘拟合并取得初步成果的基础上，进而采用更符合破裂带形态的线性立方插值函数，在破裂带主断面上，建立起一种理论与实测相结合的InSAR视线向(LOS)变化量的分解方程，得到了InSAR视线向位移分解的具有唯一性的解析解. 本文的解析法吸纳了GPS定点实测值的精度优势，利用InSAR全天候、准实时获取连续形变场的技术特点，通过数学近似，最终获得了主破裂带上连续变化的水平位错及垂直位错同震形变曲线.      

海上高分辨率地震实例研究

在莺歌海盆地油气勘探中,应用高分辨率地震资料评价ＤＦ１一１气田,探明了储集层顶面的水道冲蚀,给出了正确的气田地质模式;并应用高分辨率地震资料做气藏描述,比较准确地估算了储集层的厚度与孔隙度,减少了评价井数．应用高分辨率地震的亮点加平点技术,在原认为是空构造的ＬＤ８一１构造发现了高产天然气。理论和实践都表明,若要得到高分辨率的地震资料,提高地震子波的频率是很重要的,但保护低频比提高高频更重要。高分辨率地震波阻抗剖面比地震偏移剖面更直接、更清晰地反映地下地质情况,直接用波阻抗剖面做地震解释和油藏描述可得到更好的效果。     

油田放射性异常成因研究

为探索地表放射性异常与地下油气藏的关系，在４种类型的９个油气藏上收集了３１０口井的测井数据，同时在５个油气藏上进行了地面放射性测量．上述资料证实，异常从油层到地表具有延续性，在同一油层上方放射性异常随测点与油层距离加大而减弱，复合式油气藏各油层异常独立存在，不同类型油藏异常总体形态相似，异常幅值与油藏圈闭条件有关．这些为油气物质垂向运移及放射性异常深部形成机制提供了直接证据．在此基础上，综合前人研究成果，提出了油田放射性异常物质来源主要是生油岩的观点和油气盆地内流体微运移循环的设想，建立了油田放射性异常综合成因模式．     

多层球体对电磁场散射响应的求解方法

为建立多层球体对电磁场散射响应的求解方法,将球外入射电磁场分为TM和TE极化波两部分,通过求解Debye位方程可得到多层球体各层中两种极化场的通解形式,利用边界条件建立了TM波（或TE波）球面波阻抗（或导纳）的递推关系,借助于阻抗（或寻纳）导出了球体对TM波（或TE波）所产生散射场的一般表达式．文中具体求解了球体对平面波径向和横向偶极子场的散射响应．     

多层球体对电磁场散射响应的求解方法

为建立多层球体对电磁场散射响应的求解方法,将球外入射电磁场分为TM和TE极化波两部分,通过求解Debye位方程可得到多层球体各层中两种极化场的通解形式,利用边界条件建立了TM波（或TE波）球面波阻抗（或导纳）的递推关系,借助于阻抗（或寻纳）导出了球体对TM波（或TE波）所产生散射场的一般表达式．文中具体求解了球体对平面波径向和横向偶极子场的散射响应．     

COMPARISON STUDY OF TWO KINDS OF CODES TO MEASURE FAULT-OFFSETS BASED ON MATLAB: A CASE STUDY ON EASTERN ALTYN TAGH FAULT

Geomorphic offsets displaced by coseismic surface rupture can be analyzed to identify earthquake recurrence behavior. Therefore, obtaining a sufficient and precise along-fault offset dataset is vital to identify long-term earthquake recurrence behavior. Furthermore, knowledge of along-fault slip distribution during a single-earthquake or multi-earthquakes is important for other reasons, including a better understanding of the relationship between earthquake size and coseismic displacements, fault kinematics and fault mechanics. A recent flourish of offsets-measuring software and high-resolution topographic data together offer an unprecedented opportunity to measure high-density fault offsets. Here, we introduce and compare two kinds of most popular software, LaDiCaoz and 3D_Fault_Offsets. We describe the workflow and principle of the two codes by taking a fault-offset example on the eastern Altyn Tagh Fault. LaDiCaoz iterates over the channel morphology and position parameters and determines the summed absolute elevation difference Σ[Δ(elevation)] between both transverse profiles. The optimal horizontal offset is defined by the parameter combination that results in the least mismatch between two profiles. Compared with LaDiCaoz, the principle of 3D_Fault_Offsets is more complicated by measuring the offset in three dimensions. It mathematically identifies and represents nine of the most prominent geometric characteristics of common sublinear markers along faults in three dimensions, such as the streambed(minimum elevation), top, free face and base of channel banks or scarps(minimum Laplacian, maximum gradient, and maximum Laplacian), and ridges(maximum elevation). By calculating best fit lines through the nine point clouds on either side of the fault, the code computes the lateral and vertical offsets between the piercing points of these lines onto the fault plane, providing nine lateral and nine vertical offset measures per marker. Through a Monte Carlo approach, the code calculates the total uncertainty on each offset. Although both 3D_Fault_Offsets and LaDiCaoz are developed based on the Matlab platform, there are significant differences in principles, linear marker, software interface, repeatability, input-file types, degree of automation, adaptability, output file types, etc. In this part, we compare and summarize their features, advantages, and disadvantages. Finally, we calculate the correlation of two groups of fault-offset data derived from the two methods along the eastern ATF. By doing this, we try to explore if the two methods can be crosschecked and to study how sinuosity of the linear geomorphic markers affect the measuring results. By discussing and comparing the accuracy of the two measuring methods, we consider that LaDiCaoz is better than 3D_Fault_Offsets in accuracy aspect. In our opinion, there exist some disadvantages in the both software, and higher automation and introduction of artificial intelligence will be the future development direction.     

A CASE STUDY OF SEDIMENT IMPACT ON THE ECOSYSTEM

IINTRODUCTIONInthepast,thesedimenttransportisconsideredinmostengineeringproblemsbecauseofitssignificanteffectonreservoir,navigation,floo'dingcontrolandstabilityofcoastallines(Xuetal.,1992;Heetal.,1997).Inrecentyears,impactsofsedimentontheecologicalhabitatshavebeennotedbymanyenvironmentalists(Klaus,1997).Itseemsthattheconceptofsedimenthazardsorsedimentimpactsmustbeextendedandtheenvironmentalorecologicalimpactmustbeincludedinanycomprehensiveassessment.Forexample,thelarge-scaleconstructiono…     

2008年汶川8.0级地震地表位移场的模拟——映秀-北川断裂逆冲兼右旋走滑错动形成的地表位移场

基于Yoshimitsu Okada及Steketee的断裂位错模型,从理论上反演了汶川地震时龙门山中央主断裂中段附近区域(距离断裂30km)内由断裂逆冲兼右旋错动形成的地表位移场的空间演化,包括地表垂直和水平位移场的基本特征。模拟表明了断裂错动过程中近断裂区域地表位移场的空间演化情况,而该空间演化情况通过野外科考和有限的GPS测站数据是无法细致描述的。文中将模拟结果与野外地表破裂带的科考成果进行比较,发现两者在变化趋势上近似表现一致;且位移场在离开断裂出露处迅速衰减的特征与现有研究成果一致;同时模拟结果也表明位移场的衰减速度在下盘强于上盘。模拟结果最终表明:断裂错动形成的地表垂直位移存在较大的空间不均匀性且较大的幅值主要集中于断裂的两端,即映秀和北川附近,垂直和水平位移场在断裂端部的变化都强于断裂的中部;地表垂直位移沿走向的变化在断裂上盘强于下盘;水平位移除了在断裂两端变化比较剧烈外,在空间上分布比较均匀;总体上,除断裂端部外位移场的幅值在上盘大于下盘。     

重力异常视深度滤波及应用

为了解观测面以下位场场源体产生的异常,设计了一种波数域的滤波器,该滤波器可在无需岩石物性资料的前提下用于位场数据的反演．对位场资料用解析法和差分法上、下延拓,以消除浅表干扰源的影响,并最后下延至观测面以下,根据下延后场的极值确定场源的质心深度．理论模型的试算和实际资料的处理均获得了令人满意的效果．