羌塘盆地中部地区地壳S波速度结构及构造意义

羌塘盆地是我国最大的海相盆地.本文根据在羌塘盆地内布设的27个宽频带地震观测台站记录的远震波形数据,利用非线性复谱比反演算法得到各台站下方100 km深度范围内S波速度结构.结果表明.羌塘地区Moho深度较为平缓,平均深度为61 km;北羌塘地壳内低速层广泛分布;北羌塘具有两个较大的沉积盆地,龙尾错和白滩湖坳陷,沉积厚度分别有10 km和15 km.尽管北羌塘下地壳受到强烈的新生代火山岩作用改造,但是这种深部岩浆热作用会加速烃源岩中有机质的热演化历程.北羌塘两个盆地具有很好的油气前景.与北羌塘低速层分布相比,南羌塘下低速层更深,可能与班公怒江洋于中生代的俯冲消减及拉萨地体北向俯冲有关.     

基于Theta法改进的均衡滤波器在位场边界识别中的应用——以庐枞矿集区为例

利用观测数据确定地质体的边界位置是位场数据解译中的一项重要工作,传统的边界识别滤波器通常不能均衡深、浅部的地质体边界,近些年相关研究开始致力于发展均衡边界滤波器.本文基于Theta图法定义了新的边界识别滤波器,详述了滤波器波数域及空间域的主要计算公式,通过模型验证,该滤波器显著压制了Theta图法对深部地质体边界的放大作用,较好地平衡了深部和浅部边界.通过与传统的边界识别滤波器对比,本文定义的滤波器能够清晰且更加收敛地圈定出地质体的水平边界位置.以长江中下游成矿带庐枞矿集区为例,开展了1∶5万重磁数据的处理分析,并结合物性资料进行了讨论,结果表明:重磁数据的检测结果精确刻画了郯庐断裂带的位置;庐枞盆地的磁力数据检测边界整体与盆地的地质边缘一致,明确了边界断裂在深部倾向盆地内部;识别出庐枞盆地外围一系列环形边界,这些边界封闭区域与最新勘探发现的深部岩体及铁铜矿化体相对应,对于指导区域深部找矿工作有着重要的意义.     

吐鲁番-哈密盆地中新生代构造变形的数值模拟分析

根据地质、地震剖面、重力等资料及野外实地考察结果，分析了新疆吐哈盆地内各块体的力学性质和运动方式，在此基础上提出了吐哈盆地的发育模式。根据地质模型，并充分考虑了地质时间因素，建立了包括大区域平面模型、南北向剖面模型和台北凹陷三维模型等计算模型，用有限元数值计算方法模拟了古构造应力场和构造变形场。结果表明：吐哈盆地自侏罗纪以来一直受到近南北向的挤压，但盆地各段边界上的作用方向略有变化。模拟结果还显示出盆地的不对称形态以及主应力方向随深度的变化     

联合星载GPS和KBRR星间测速数据反演GRACE时变重力场模型

高精度GRACE卫星时变重力场反演一直是卫星重力测量中的难题.为了恢复高精度的时变地球重力场模型,本文联合GRACE卫星的星载GPS和KBR星间测速观测数据,在对GRACE卫星进行精密定轨的同时,解算出60阶月平均地球重力场模型.通过对GRACE卫星的定轨精度、星载GPS相位和KBR星间测速数据的拟合残差以及时变地球重力场模型解算精度等分析,表明:(1)与美国宇航局喷气推进实验室(JPL)发布的约化动力学精密轨道相比,本文确定GRACE卫星轨道三维位置误差小于5 cm.(2)星载GPS相位数据拟合残差为5~8 mm,KBR星间测速数据拟合残差为0.18~0.30μm·s~(-1).(3)解算的月平均重力场模型与美国德克萨斯大学空间研究中心(CSR)、德国地学研究中心(GFZ)和JPL发布的RL05模型精度接近,时变信号在全球范围内具有很好的空间分布一致性.通过计算亚马逊流域和长江流域的水储量变化,本文与上述三个机构的计算结果无明显差异,且相关系数均达0.9以上.可见,本文建立的卫星轨道与重力场同解算法具有反演高精度GRACE时变重力场能力,为我国卫星重力场反演提供了重要的技术支持.     

由多个地震震源机制解求川滇地区平均应力场方向

基于震源断层面解的空间取向和断层滑动方向,写出相应力轴张量在地理坐标系中的表达式,进而给出计算平均力轴张量及主值的方法,即通过求解相应的本征方程得到. 对使用多个震源机制解的T,B,P轴参数计算平均应力场的方法,以用滑动方向拟合法反演富蕴、唐山地区平均应力场数据进行验证,二者结果一致. 选择具有地震构造意义的地块或地震带内大量地震的震源机制解研究区域平均应力场. 根据川滇13个地震带（区）的256次中强地震的震源机制解,给出了各带（区）应力张量的定量分析结果. 该方法算法简便,使由大量地震震源机制解资料分析构造应力场的方法走向定量化.      

GSFC（9/80）地磁模型的评价和中国地磁台站异常值的确定

本文利用1960-1980年中国地磁台站资料检验了GSFC（9/80）地磁模型,结果表明该模型是成功的。据此通过地磁台站1960-1980年的观测值与GSFC（9/80）模型计算值的比较,确定了中国地磁台站的异常值。本文还建议在基本磁场的研究中从台站观测年均值中减去较精确测定的台站磁异常值,以提高资料的精确度。     

东海陆架盆地及其周边海域地质、地球物理场特征

分东海陆架盆地及其周边海域的重、磁场特征、对研究区域的地球物理场资料进行了小波变换与多尺度分布与计算，计算了研究了海域重，磁资料的1－4阶小波变换逼近、细节、根据分析与计算可知，在东海陆架盆地及其周边域，重、磁异常的主要特征没有内陆盆地那样明显，但盆地内的绝大部分地区均位于异常相对缓变区，在盆地内部有时有局部异常的圈闭及剧变区，但总的特征为缓变的，而周边区域多为相对剧变区，但也存在局部区域的缓变, 、磁场特征表明，盆地基底相对盆地周边区域主要由密度较大，磁性中的岩石组成，因此盆地中部以高正磁异常为主，布格重力异常 对周边地区高，东海陆架盆地及其周边海域的重力场的主要走向为NNE向，局部重力异常较为平缓，磁异常走向没有重力异常明显，但主要也为NNE向，重、磁场的分布特征表明，东海陆架盆地及其邻域的地质构造为大陆架构造的延伸。     

关中地区垂直形变场和应力场的特征

本文利用关中盆地1972—1986年的4期大面积垂直形变资料,采用经典静态平差方法,求出该盆地相对于1972年的垂直位移量,并作出位移等值线图。据此讨论了近年来该盆地的垂直形变场的变化特征。结果表明,该地区现今的地壳垂直形变基本上继承了第四纪时期的特征,但在盆地西部有差异。本文将关中盆地简化为一个平面力学模型,采用有限元法反演关中盆地1976年、1980年和1986年相对于1972年的应力值,研究了该盆地的应力场分布特征。     

大地电磁资料偏移成像

在研究Zhdanov的二维电磁场偏移成像的基础上, 针对它在某些方面的不足, 研究并提出了改进算法. 首先, 用一维反演结果确定初始背景电阻率, 然后在延拓计算过程中, 利用迭代算法进行逐层修正、 延拓计算, 从而获得较准确的电阻率; 第二, 研究了有限差分方程的改进算法, 根据大地电磁场变化的特点, 设计了纵向以等比级数变网格的迭代算法. 通过以上两方面技术的改进, 经理论模型和实际资料试算验证, 取得了较好的结果.      

南海北部陆架西区盆地地质、地球物理场特征及其深部结构

分析了南海北部陆架西区盆地的地质、地球物理场特征，计算了研究海域重、磁资料的1阶小波细节、4阶小波逼近变换。根据分析与计算可知，研究区的布格重力异常以北西低的负值，东南高的正值为特征。在东部及东南部异常等值线走向为北东；西部异常等值线以北西走向为特征；西北地区异常以北东东、北东走向的局部等值线圈闭为特征。磁场的展布十分复杂，按磁异常的变化程度可分为三个变化区，即磁异常平静区、剧变区及缓变区。磁异常的平静区位于研究区的西部，即莺歌海盆地所在位置，这一带磁异常等值线极为稀疏，异常值为负背景异常。剧变区位于海南岛，该地区的磁异常变化极为剧烈，异常特征以局部小圈闭为特征，等值线分布密集。磁异常的平缓区位于平静区及剧变区之外的其它地区。琼东南盆地、北部湾盆地的磁异常具有此特征。根据重、磁场资料以及南海北部盆地钻井取样的测试结果、同时参考穿越南海地学断面的结果，对研究区的地壳结构进行了反演计算，计算表明南海陆架盆地区域地壳结构较为复杂，研究区的地壳厚度在22－33km之间，总的趋势由陆向洋地壳厚度逐渐减薄，反映出该区域地壳具有陆壳、拉伸陆壳、过渡壳的性质，同时存在有上地幔隆起区及凹陷区。磁性底界面厚度在17－24km之间变化，其中在莺歌海盆地较深，在海南岛地区磁性界面较浅。     

Joint parameter estimation from magnetic resonance and vertical electric soundings using a multi‐objective genetic algorithm

Magnetic resonance sounding (MRS) has increasingly become an important method in hydrogeophysics because it allows for estimations of essential hydraulic properties such as porosity and hydraulic conductivity. A resistivity model is required for magnetic resonance sounding modelling and inversion. Therefore, joint interpretation or inversion is favourable to reduce the ambiguities that arise in separate magnetic resonance sounding and vertical electrical sounding (VES) inversions. A new method is suggested for the joint inversion of magnetic resonance sounding and vertical electrical sounding data. A one‐dimensional blocky model with varying layer thicknesses is used for the subsurface discretization. Instead of conventional derivative‐based inversion schemes that are strongly dependent on initial models, a global multi‐objective optimization scheme (a genetic algorithm [GA] in this case) is preferred to examine a set of possible solutions in a predefined search space. Multi‐objective joint optimization avoids the domination of one objective over the other without applying a weighting scheme. The outcome is a group of non‐dominated optimal solutions referred to as the Pareto‐optimal set. Tests conducted using synthetic data show that the multi‐objective joint optimization approximates the joint model parameters within the experimental error level and illustrates the range of trade‐off solutions, which is useful for understanding the consistency and conflicts between two models and objectives. Overall, the Levenberg‐Marquardt inversion of field data measured during a survey on a North Sea island presents similar solutions. However, the multi‐objective genetic algorithm method presents an efficient method for exploring the search space by producing a set of non‐dominated solutions. Borehole data were used to provide a verification of the inversion outcomes and indicate that the suggested genetic algorithm method is complementary for derivative‐based inversions.     

Modeling of the Bosphorus exchange flow dynamics

The fundamental hydrodynamic behavior of the Bosphorus Strait is investigated through a numerical modeling study using alternative configurations of idealized or realistic geometry. Strait geometry and basin stratification conditions allow for hydraulic controls and are ideally suited to support the maximal-exchange regime, which determines the rate of exchange of waters originating from the adjacent Black and Mediterranean Seas for a given net transport. Steady-state hydraulic controls are demonstrated by densimetric Froude number calculations under layered flow approximations when corrections are applied to account for high velocity shears typically observed in the Bosphorus. Analyses of the model results reveal many observed features of the strait, including critical transitions at hydraulic controls and dissipation by turbulence and hydraulic jumps. It is found that the solution depends on initialization, especially with respect to the basin initial conditions. Significant differences between the controlled maximal-exchange and drowned solutions suggest that a detailed modeling implementation involving coupling with adjacent basins needs to take full account of the Bosphorus Strait in terms of the physical processes to be resolved.     

内磁层电流体系的地磁效应

对Ⅱ区场向电流及其伴随的部分环电流和电离层电流组成的内磁层三维电流体系（PRFI电流系）的磁场效应进行了数值计算.这一三维电流体系在中低纬度地面产生的磁场呈现出特殊的纬度分布:X分量几乎不随纬度变化,Y分量随纬度增高近似呈线性变化.这些特征明显不同于对称环电流的磁场分布特征（X∝ cosφ,φ是纬度,Y=0）,也不同于DP2、S_q、L等电流体系的磁场分布特征.利用这一特征我们可以从地磁台子午链观测到的磁场扰动中分离出PRFI电流系的贡献.用1989年3月磁暴的实例检验了上述模型,观测结果与理论结果符合得很好.分析结果还表明,最大的D_st指数并不一定对应着最强的对称环电流.     

四川盆地综合地质、地球物理研究

对四川盆地的地质、地球物理场特征进行了分析与探讨，对重力、航磁数据进行了延拓、求导，并反演计算了该盆地莫霍界面与磁性项、底界面，为认识盆地的断裂系统、地壳深部结构提供了基础，在此基础上收集了四川盆地的人工地震与测井资料，对该盆地进行了盆地沉降史的动态模拟，模拟表明，四川盆地在65Ma以后，属于稳定沉积阶段。     

GEOMORPHIC FEATURES OF THE SHULE RIVER DRAINAGE BASIN IN QILIANSHAN AND ITS INSIGHT INTO TECTONIC IMPLICATIONS

Because of the strong uplift of the Qilian Shan since late Cenozoic,the drainage basins that are derived from the mountains have undergone strong tectonic deformation.So the typical geomorphology characteristics of these drainage basins may indicate the strong tectonic movement in the region.For example,the Shule River drainage basin,which originates from the western part of the Qilian Shan owns unique geomorphology characteristics which may indicate the neotectonic movement. Stream networks of the Shule drainage basin extracted from the DEM data based on GIS spatial analysis technology are graded into five levels using Strahler classification method.Four sub-catchments,numbered 1,2,3 and 4 are chosen for detailed analysis.Furthermore,the four sub-catchments,the hypsometric integral curves,Hack profiles,SL index and average slope of the Shule drainage basin are determined by GIS tools.In addition,we analyzed the slope spectrum of the Shule drainage basin. The average elevation of the Shule drainage basin is very high,however,the slope of the drainage basin is very low,the gentle slope occupies so large area proportion that the slope spectrum shows a unimodal pattern and a peak value is in low slope region (0°~5°),so tectonic movement has a strong influence on the drainage basin.Under the intensive impact of the tectonic movement of the active fault and regional uplift,the hypsometric integral curve is sigmoid,revealing that the Shule drainage basin is in the mature stage.The Hack profile is on a convex,the longitudinal profile is best fitted by linear fitting and the abnormal data of the SL index of the Shule River has a good fit with the section through which the active fault traverses,that means the tectonic movement of the active fault has strong influence on the river's SL index.It is worth noting that lithologic factors also have great impact on the river geomorphology in some sections. According to the above analysis,we recognize that in the interior of active orogen,the evolution of river geomorphology usually is influenced by tectonic movement and reveals the regional neotectonics in turn.     

Tomographic imaging of P- and S-wave velocity structure beneath Costa Rica

45287 P-wave and 26813 S-wave arrival times from the data base of the Costa Rica network have been tomographically inverted to image the structure beneath Costa Rica. A regularized recursive least squares inverse method was used to produce the high resolution and minimum variance model parameter estimates. The first arrival times are calculated using a finite difference technique, which allows for flexible parameterization of the velocity model and easy inclusion of topography and source-receiver geometry. The P wave velocity structure and hypocenters are determined simultaneously, while the S wave velocity structure is determined using the relocated seismicity and an initial model derived from the P wave model assuming an average P to S wave velocity ratio of 1.78. The most prominent features in the inverted model are a low velocity structure under the volcanic chain in the center of the country, which is related to the hot material connected with the active volcanoes; and a high velocity zone in the mantle, which is related to the Cocos plate subducted under Costa Rica.     

利用黄土剖面的古土壤年龄研究毛毛山断裂的滑动速率

利用１４Ｃ、热释光（ＴＬ）样品年代及扩散方程计算结果，结合区域黄土剖面中古土壤年龄，对毛毛山地区晚第四纪各级地貌年龄进行了对比研究。根据毛毛山活动断裂水平位移和垂直位移分布明显的分组特征，求得毛毛山断裂带不同段落不同时段的平均滑动速率。大约自中更新世晚期以来，毛毛山断裂走滑段的平均水平滑动速率为２．３～３．９ｍｍ／ａ，垂直滑动速率为０．０７～０．１９ｍｍ／ａ；天祝盆地倾滑段垂直滑动速率为０．１１～０．８６ｍｍ／ａ。沿断裂带滑动速率具明显的非均匀性特点，表现为自东向西水平位移具累积滑动亏损特征，垂直位移则具补偿性     

3D acoustic waveform inversion in the Laplace domain using an iterative solver

In this paper we propose a 3D acoustic full waveform inversion algorithm in the Laplace domain. The partial differential equation for the 3D acoustic wave equation in the Laplace domain is reformulated as a linear system of algebraic equations using the finite element method and the resulting linear system is solved by a preconditioned conjugate gradient method. The numerical solutions obtained by our modelling algorithm are verified through a comparison with the corresponding analytical solutions and the appropriate dispersion analysis. In the Laplace‐domain waveform inversion, the logarithm of the Laplace transformed wavefields mainly contains long‐wavelength information about the underlying velocity model. As a result, the algorithm smoothes a small‐scale structure but roughly identifies large‐scale features within a certain depth determined by the range of offsets and Laplace damping constants employed. Our algorithm thus provides a useful complementary process to time‐ or frequency‐domain waveform inversion, which cannot recover a large‐scale structure when low‐frequency signals are weak or absent. The algorithm is demonstrated on a synthetic example: the SEG/EAGE 3D salt‐dome model. The numerical test is limited to a Laplace‐domain synthetic data set for the inversion. In order to verify the usefulness of the inverted velocity model, we perform the 3D reverse time migration. The migration results show that our inversion results can be used as an initial model for the subsequent high‐resolution waveform inversion. Further studies are needed to perform the inversion using time‐domain synthetic data with noise or real data, thereby investigating robustness to noise.     

小震震源参数与震源机制的联合测定

本文采用一个非均匀破裂圆盘位错模式,考虑到介质-仪器系统的影响,计算了它在不同幅射方向的P波初动位移与半周期,以及几组用于测定震源参数的理论曲线。在采用1/4周期测定介质的Q值和球面拟合震源机制解的基础上,给出了一个小震震源参数与震源机制的联合测定方法。该方法的实质是,将震源断层面作为被测参数的一部分,并同时考虑它对其它被测参数的影响,从而获得单个小震震源参数与震源机制的相应解。最后,根据新丰江水库区四个台的小震资料作了实测与分析,所得结果为:介质的平均Q值约450,地震矩1019——1020达因厘米,平均应力降约35巴;其震源机制解亦反映了本区构造应力场的主要特征。      

中国岩石层底面地幔对流状态的探讨

本文应用卫星测量地球重力场的球谐函数系数和地幔流体力学方程,计算了中国地区岩石层下的地幔对流状态。计算中选用了三种不同的模式:模式Ⅰ,完全场（2-30阶）;模式Ⅱ,高阶场（13-25阶）;模式ⅢI,低阶场（4-18阶）。其目的是探求一种与我国地质和地球物理实测资料相符合的最佳模式。考虑到岩石层厚度与所假定的地幔对流区厚度相比是很薄的,因而可以把地幔流对岩石层底面产生的拖曳力近似地看作岩石层内构造力。将上述三个模式分别与中国及邻区的131个地震（1961-1977,M≥5）震源机制解和近代地质构造形变等资料进行比较。结果表明,模式Ⅰ与我国大多数地区的实测资料符合得最好,该模式所显示的地幔对流特征与相应各地区近代构造运动特点相符。模式Ⅰ为探讨近代构造运动的力源问题,提供了新的线索。