深地震反射剖面构造信息识别研究

随着我国深部探测计划的实施,近年来在中国大陆地区实施了多条深地震反射剖面,获得了大量揭示地壳上地幔精细结构的宝贵资料. 为了更充分地挖掘这些资料中所包含的深部构造信息,本文在研究和总结前人工作的基础上,提出一种快速识别深地震反射剖面构造格架的方法. 该方法通过数据预处理、强振幅提取、中值滤波、对象识别、连续性计算和连续性滤波实现对深地震反射剖面构造格架识别. 同时,该方法还可通过对象倾角计算对复杂区域进行属性分析.     

重力资料识别鄂尔多斯盆地线性构造方法研究

目的利用重力异常数据微弱信息挖掘技术识别线性构造的方法,并介绍其实现过程.方法采用水平梯度法和二次垂直导数等方法提取微弱信息,利用数字图象显示技术成图.结果识别出鄂尔多斯盆地线性构造多条,分析对比鄂尔多斯盆地区域地质构造资料及重力解释成果,具有较好的吻合性.为深入研究该地区线性构造补充了新的证据.结论该方法提取了区域重力异常资料中的微弱信息,弥补用等值线图示构造信息的不足,有助于线性构造的分析.水平梯度法(最大模)和二次导数处理方法与数字图象成图相结合取得良好的解释效果.     

湖北襄樊—福建罗源的大地电磁测量与地壳上地幔电性特征研究

湖北襄樊—福建罗源一千多公里的剖面上作了23个大地电磁测深点。该剖面穿过不同的地质构造单元。 大地电磁测深的结果清楚地反映了剖面内地壳和上地幔的电性结构特征。例如构造单元的形态、边界、断层位置、地壳和上地幔的高导层发育情况等。结果对于了解该区的深部构造背景、探讨深浅构造之间的关系提供了有益信息     

基于数字高程模型的南天山地貌特征研究

基于STRM-DEM((90 m)数据,利用ArcGIS空间分析模块,通过人工判定法、最大高差法以及最大高差与面积比法联合对比,获取研究区最佳分析窗口为31×31,面积约为7. 78 km~2,根据获取的最佳窗口提取出西南天山地区地形高程(最大高程、最小高程、平均高程)和地形起伏度。分析认为:西南天山山势陡峻,总体呈中东部高、西部低的地貌特征,且以高地形起伏度为主,沿西南天山走向方向,高山区地形起伏度自东向西逐渐降低;高海拔地区呈北东—南西向展布规律,反映了该地区遭受了南北向强烈挤压隆升作用;并通过条带状地形剖面分析研究区夷平面地貌特征,认为西南天山地区发育4级层状地貌面,反映了该地区的多次构造运动。     

高起伏区无人机摄影测量精度评估?以西秦岭光盖山-迭山断裂为例

经过近10年的迅速发展,无人机摄影测量已成为活动构造研究的常用方法之一。但对于无人机摄影测量的精度评估,尤其是高起伏地区的精度评估存在不足。为此,选择白龙江北岸光盖山-迭山断裂沿线的黑峪寺、化马村,开展无人机摄影测量,并构建正射影像（DOM）和数字地表模型（DSM）,配合差分GPS测绘进行校正和精度验证。通过对比实测控制点和图像提取点分析点精度,通过对比实测剖面与提取剖面分析剖面精度。研究结果表明,未经控制点校正的图像提取点与实测点存在较大误差,水平误差为5～8 m,垂直误差为几十米至上百米,但通过少数控制点校正后,点精度可达20 cm以内;6条实测剖面与提取剖面（提取自控制点校正后的图像）平均垂直精度总体为分米级,即0.16～0.65 m,标准差为0.13～0.69 m,略低于低起伏区的精度,对于测量条件恶劣的高起伏区,该精度是可接受的;异常高的垂直误差常出现在地形突变、低矮植被密集、行走困难等测量条件不理想位置。图像控制点中心点的准确识别、提取剖面线的修正准确性等因素也会影响精度评估的可靠性。     

云南南部地区深部电性结构特征研究

在云南南部地区布设了一条孟连-罗平的北东向大地电磁测深剖面,以开展该地区的深部电性结构探测和孕震环境探查.沿该剖面进行了114个大地电磁测深点的观测,经过对观测资料的远参考Robust处理、定性分析和二维反演,得到了沿该剖面地壳、上地幔电性结构模型,从模型的电性结构特征进一步探讨了剖面穿过的3个地震区的深部地震孕育环境.研究结果表明：沿剖面的地壳上地幔电性结构反映出与区域地质构造资料基本一致的构造特征;该区的三个强震带地球深部都存在壳内低阻体,地震发生在电阻率梯度带上;断裂带的两侧块体介质的电阻率差异是强震活动带重要的深部背景.     

龙门山中北段流域地貌特征及其构造意义

龙门山中北段位于青藏高原东缘,该区作为高原向东扩展的前缘部位,其地形与河流水系的演化记录了高原隆升与挤压扩展及其气候环境效应的各种信息。龙门山中北段构造活动有明显差异,从中段逆冲为主转化为北段的走滑为主,本文采用定量化地貌参数从构造地貌的角度揭示了区域构造活动的差异。龙门山中北段地貌因子（坡度、地形起伏度和条带状剖面）的阶梯状分布特点,显示了高原扩展的逆冲推覆特征,在中央断裂处构造抬升作用最强,同时显示出了南北向构造活动减弱的趋势,由中段的逆冲转换为北段逆冲兼走滑的形式。北川-映秀断裂两侧流域的HI值也显示了断裂上盘高、下盘低、沿走向减弱的趋势。综合分析认为,本区构造活动是地貌演化的主控因素,龙门山中北段地形存在差异,北川-映秀断裂两侧的小流域地貌指数分析显示,构造抬升活动自南向北减弱,中段以逆冲为主,北段为逆冲兼走滑。     

若尔盖盆地-西秦岭造山带结合部位深反射资料的静校正和去噪技术

深地震反射剖面技术是探测岩石圈精细结构的有效手段.通常情况下工区地质情况复杂,尤其在盆山结合部位,地表地形起伏大,地下构造复杂,其深地震反射资料具有低信噪比、干扰强、构造复杂等特点,给后续处理和解释造成很大困难,因此获得真实的叠加剖面是地质解释的前提和基础.复杂地区低信噪比深地震反射资料处理的关键是做好静校正和去噪工作.本文以若尔盖盆地－西秦岭造山带接合部位深地震反射资料作为例,通过方法试验和参数测试,找到适合该工区的静校正方法和去噪技术,得到较好的处理结果,为揭示若尔盖盆地－西秦岭造山带结合部位的细结构提供了可靠的依据.     

一种中浅地层剖面高精度成像的处理方法

在海洋工程物探中浅地层剖面调查中,调查船的上下起伏会对地层剖面的采集造成影响,即便是在姿态传感器的补偿下,仍然会存在海底及地层波状起伏的假成像现象.载体搭载的地层剖面系统记录的多为相对于载体的数据,而不是相对于海面的深度.在实际工程作业中,与地层剖面同时采集的多波束数据通常包含有高精度的地形数据,本文研究的目的是通过采集不同水深、不同区块、不同类型的工程物探资料,将测深和剖面探测两种调查方式结合,研究采用多波束地形数据来校正地层剖面的海底,消除涌浪、姿态、位置、滤波等参数引起的成像偏差,精准呈现真实的海底地形,实现高精度的剖面成像,有效提升工程物探资料研究评价的准确性和可靠性,为以后在类似区域进行海底浅层灾害地质研究提供参考和借鉴.     

陕渝黔桂1800km超长探测剖面重力异常场特征及深部地壳结构探榷

通过对跨越华北克拉通、秦岭造山带、扬子克拉通几大地质构造单元的,北起陕西榆林经秦岭过重庆鱼泉、贵州贵阳并向南到广西凭祥全长1810km超长重力探测剖面的数据进行处理分析和解释,构建了沿剖面的二维地壳密度结构模型,并详细分析了沿剖面壳内各界面与Moho界面展布的深部结构和构造特征、构划出了该剖面的深部断裂分布,探讨了剖面辖区跨越的克拉通、造山带、接触带或耦合带等一系列的区域构造的差异,同时对其可能的地质构造含义进行初步了解释,以期能对深化认识该剖面跨越地区的地壳结构、构造单元划分及动力学等研究,提供相关重力场的依据.     

Plumes,subaxial pipe flow,and topography along the Mid-Oceanic Ridge

If plate thickness depends on crustal age, the region of extensive partial melting below the spreading axis will be wider around fast-spreading ridges. The melt region creates a subaxial conduit channeling partial melts away from ridge-centered hot spots. The channel is here modeled by an elliptical pipe of semiminor (vertical) axis 2 × 10⁶ cm (20 km) and semimajor (horizontal) axis KS, where S is spreading half-rate (cgs) and K is a constant of magnitude 10¹⁴ to 10¹⁵ seconds. This simple analytical model is used to explain the observation that maximum hot spot elevations on the Mid-Oceanic Ridge fall dramatically with increasing spreading rate (there are no Icelands or Afars on the East Pacific Rise!). A hot spot under a fast-spreading ridge has a broad pipe in which to discharge its partial melts; hence, only a slight topographic gradient and a low elevation is needed to discharge the mass flux rising out of the deeper mantle at the hot spot center. A second factor is that partial melts are “used up” faster in the accretion process on fast-spreading ridges. In the simple analytical model, both factors operating together explain the rapid fall of hot spot heights with increasing spreading half-rate. This result indirectly helps confirm the idea of horizontal pipe flow below the Mid-Oceanic Ridge.A theoretical topographic profile through a hot spot on the Mid-Oceanic Ridge is derived from the assumption that the pressure — i.e., topographic — gradient at a distance x from the hot spot is sufficient to supply all the accreting lithosphere downstream of x, out to x_n, the limit of topographic hot spot influence. The predicted profile is quadratic in x and concave upward, and resembles several observed profiles where neighboring hot spots are not so close as to confuse the profiles. Some observed profiles are more nearly linear or even convex upward. This could be explained, for example, by downstream increases in viscosity or decreases in pipe dimensions.A hot spot on a ridge spreading at much less than 1 cm/yr half-rate would produce an enormous elevation of the ridge axis, according to our model, because the pipe would be very narrow. Such a large topographic high would create a large gravity potential which would cause the plates to move apart faster, thereby widening the pipe, and reducing the topographic high. The system of ridges and hot spots may thus be self-regulating with respect to plate speeds; this could explain why spreading half-rates on the Mid-Oceanic Ridge are in many areas as low as 1.0 cm/yr but very rarely as low as 0.5 cm/yr.     

APPLICATION OF DEM GENERATION TECHNOLOGY FROM HIGH RESOLUTION SATELLITE IMAGE IN QUANTITATIVE ACTIVE TECTONICS STUDY: A CASE STUDY OF FAULT SCARPS IN THE SOUTHERN MARGIN OF KUMISHI BASIN

Traditional method to generate Digital Elevation Model (DEM)through topographic map and topographic measurement has weak points such as low efficiency, long operating time and small range. The emergence of DEM-generation technology from high resolution satellite image provides a new method for rapid acquisition of large terrain and geomorphic data, which greatly improves the efficiency of data acquisition. This method costs lower compared with LiDAR (Light Detection and Ranging), has large coverage compared with SfM (Structure from Motion). However, there is still lack of report on whether the accuracy of DEM generated from stereo-imagery satisfies the quantitative research of active tectonics. This research is based on LPS (Leica Photogrammetry Suit)software platform, using Worldview-2 panchromatic stereo-imagery as data source, selecting Kumishi Basin in eastern Tianshan Mountains with little vegetation as study area. We generated 0.5m resolution DEM of 5-km swath along the newly discovered rupture zone at the south of Kumishi Basin, measured the height of fault scarps on different levels of alluvial fans based on the DEM, then compared with the scarp height measured by differential GPS survey in the field to analyze the accuracy of the extracted DEM. The results show that the elevation difference between the topographic profiles derived from the extracted DEM and surveyed by differential GPS ranges from -2.82 to 4.87m. The shape of the fault scarp can be finely depicted and the deviation is 0.30m after elevation correction. The accuracy of measuring the height of fault scarps can reach 0.22m, which meets the need of high-precision quantitative research of active tectonics. It provides great convenience for rapidly obtaining fine geometry, profiles morphology, vertical dislocations of fault and important reference for sites selection for trench excavation, slip rate, and samples. This method has broad prospects in the study of active tectonics.     

On the exact shape of the horizontal profile of a topographically rectified tidal flow

Abstract

Starting from the nonlinear shallow water equations of a homogeneous rotating fluid we derive the equation describing the evolution of vorticity by a fluctuating bottom topography of small amplitude, using a multiple scale expansion in a small parameter, which is the topographic length scale relative to the tidal wave length. The exact response functions of residual vorticity for a sinusoidal bottom topography are compared with those obtained by a primitive perturbation series and by harmonic truncation, showing the former to be invalid for small topographic length scales and the latter to be only a fair approximation for vorticity produced by planetary vortex stretching. In deriving the exact shape of the horizontal residual velocity profile at a step-like break in the bottom topography, it is shown that the Lagrangian profile only exists in a strip having the width of the amplitude of the tidal excursion at both sides of the break, and that it vanishes outside that interval. Moreover, in the limit of small amplitude topography at least, it vanishes altogether for the generation mechanism by means of planetary vortex stretching. The Eulerian profile is shown to extend over twice the interval of the Lagrangian profile both for production by vortex stretching and by differential bottom friction. These finite intervals over which the residual velocity profiles exist for a step-like topography are not reproduced by harmonic truncation of the basic equation. This method gives exponentially decaying profiles, indicating spurious horizontal diffusion of vorticity. In terms of orders of magnitude, the method of harmonic truncation is reliable for residual velocity produced by vortex stretching but it overestimates the residual velocity produced by differential bottom friction by a factor 2.     

Parallels plane projection and its geometric features

A new equivalent map projection called the parallels plane projection is proposed in this paper. The transverse axis of the parallels plane projection is the expansion of the equator and its vertical axis equals half the length of the central meridian. On the parallels plane projection, meridians are projected as sine curves and parallels are a series of straight, parallel lines. No distortion of length occurs along the central meridian or on any parallels of this projection. Angular distortion and the proportion of length along meridians (except the central meridian) introduced by the projection transformation increase with increasing longitude and latitude. A potential application of the parallels plane projection is that it can provide an efficient projection transformation for global discrete grid systems.     

Topographic stress and rock fracture: a two‐dimensional numerical model for arbitrary topography and preliminary comparison with borehole observations

Theoretical calculations indicate that elastic stresses induced by surface topography may be large enough in some landscapes to fracture rocks, which in turn could influence slope stability, erosion rates, and bedrock hydrologic properties. These calculations typically have involved idealized topographic profiles, with few direct comparisons of predicted topographic stresses and observed fractures at specific field sites. We use a numerical model to calculate the stresses induced by measured topographic profiles and compare the calculated stress field with fractures observed in shallow boreholes. The model uses a boundary element method to calculate the stress distribution beneath an arbitrary topographic profile in the presence of ambient tectonic stress. When applied to a topographic profile across the Susquehanna Shale Hills Critical Zone Observatory in central Pennsylvania, the model predicts where shear fractures would occur based on a Mohr–Coulomb criterion, with considerable differences in profiles of stresses with depth beneath ridgetops and valley floors. We calculate the minimum cohesion required to prevent shear failure, C_min, as a proxy for the potential for fracturing or reactivation of existing fractures. We compare depth profiles of C_min with structural analyses of image logs from four boreholes located on the valley floor, and find that fracture abundance declines sharply with depth in the uppermost 15 m of the bedrock, consistent with the modeled profile of C_min. In contrast, C_min increases with depth at comparable depths below ridgetops, suggesting that ridgetop fracture abundance patterns may differ if topographic stresses are indeed important. Thus, the present results are consistent with the hypothesis that topography can influence subsurface rock fracture patterns and provide a basis for further observational tests. Copyright © 2014 John Wiley & Sons, Ltd.     

Discriminating between the roles of late Pleistocene palaeodischarge and geological‐topographic inheritance in fluvial longitudinal profile and channel development

This paper investigates the influences of palaeohydrology and geological‐topographic inheritance in shaping the channel of the lower River Suir, southeast Ireland. Results of acoustic surveys of the lower River Suir and Waterford Harbour reveal two scales of pseudo‐cyclic river bedforms. Longitudinal elevation profiles of the geological topography (undulating bedrock and till‐mantled bedrock) bounding the present floodplain swath reveal pseudo‐cyclicity in that terrain too. Spectral and statistical analyses are used to quantify the cyclicity of the long profile and geological‐topographic series. These methods show that the dominant cyclicity of the long profile reflects autocorrelation more than inheritance of cyclicity from the bounding geological topography. The cyclicity of the long profile mainly reflects a hydraulic control on pool‐spacing, although some cyclicity probably has been inherited from the geological‐topography. Channel‐forming palaeodischarge is estimated based on the dominant pool‐spacing revealed by spectral analysis, validated using relationships between meander wavelength, channel cross‐sectional geometry and hydraulically‐informed discharge reconstruction. The palaeodischarge estimates are in close agreement and are two orders of magnitude greater than present flood maxima. Significantly, these palaeodischarge estimates also agree closely with palaeodischarge calculated for the submerged Pleistocene palaeochannel that extends across the near‐shore continental shelf from Waterford Harbour. The pool‐sequence of the lower Suir and the submerged palaeochannel represent a former land‐system that was active during a period of low relative sea level during the last glacial. More broadly, the paper offers insights into the landscape evolution of formerly glaciated regions that experienced very wide discharge variability during and after the transition from glacial to interglacial regimes, in a context of complex relative sea level change. Copyright © 2017 John Wiley & Sons, Ltd.     

Parallels plane projection and its geometric features

A new equivalent map projection called the parallels plane projection is proposed in this paper. The transverse axis of the parallels plane projection is the expansion of the equator and its vertical axis equals half the length of the central meridian. On the parallels plane projection, meridians are projected as sine curves and parallels are a series of straight, parallel lines. No distortion of length occurs along the central meridian or on any parallels of this projection. Angular distortion and the proportion of length along meridians (except the central meridian) introduced by the projection transformation increase with increasing longitude and latitude. A potential application of the parallels plane projection is that it can provide an efficient projection transformation for global discrete grid systems.

     

3-D shell analysis of cylindrical underground structures under seismic shear (S) wave action

The 3-D shell theory is employed in order to provide a new perspective to earthquake-induced strains in long cylindrical underground structures, when soil-structure interaction can be ignored. In this way, it is possible to derive analytical expressions for the distribution along the cross-section of axial, hoop and shear strains and also proceed to their consistent superposition in order to obtain the corresponding principal and von Mises strains. The resulting analytical solutions are verified against the results of 3-D dynamic FEM analyses. Seismic design strains are consequently established after optimization of the analytical solutions against the random angles which define the direction of wave propagation relative to the longitudinal structure axis, the direction of particle motion and the location on the structure cross-section. The basic approach is demonstrated herein for harmonic shear (S) waves with plane front, propagating in a homogeneous half-space or in a two layer profile, where soft soil overlays the bedrock.     

Mohr圆分析方法在川西-藏东MT剖面资料解释中的应用

用Mohr圆分析法对川西—藏东全测区 76个大地电磁测深点的资料进行了分析 ,结果表明 :测区以 2D区域构造为主要特征 ,横向各向异性不太强 ;区域构造走向大致呈近SN向 ,最大偏离不超过± 30°。文中介绍了 4个测点的分析结果     

On the effect of topography on surface wave propagation in the ambient noise frequency range

Due to the increasing popularity of analyzing empirical Green’s functions obtained from ambient seismic noise, more and more regional tomographical studies based on short-period surface waves are published. Results could potentially be biased in mountainous regions where topography is not small compared to the wavelength and penetration depth of the considered waves. We investigate the effect of topography on the propagation of short-period Rayleigh waves empirically by means of synthetic data using a spectral element code and a 3-D model with real topography. We show that topography along a profile through the studied area can result in an underestimation of phase velocities of up to about 0.7% at the shortest investigated period (3 s). Contrary to the expectation that this bias results from the increased surface distance along topography, we find that this error can be estimated by local topographic contrasts in the vicinity of the receiver alone. We discuss and generalize our results by considering topographic profiles through other mountain ranges and find that southern Norway is a good proxy to assess the topography effect. Nevertheless, topographic bias on phase velocity measurements is in general not large enough to significantly affect recovered velocity variations in the ambient noise frequency range.