基于局部采样MCMC方法的时移探地雷达反演

马尔科夫链蒙特卡洛方法(MCMC)是一种启发式的全局寻优算法,可以用来解决概率反演的问题.基于MCMC方法的反演不依赖于准确的初始模型,可以引入任意复杂的先验信息,通过对先验概率密度函数的采样来获得大量的后验概率分布样本,在寻找最优解的过程中可以跳出局部最优得到全局最优解.MCMC方法由于计算量巨大,应用难度较高,在地...     

各向异性扩散在电磁测深曲线平滑中的应用

电磁类地球物理方法由于极易受到各类噪声的干扰,使得估计的视电阻率曲线或相位曲线发生畸变,严重地影响了反演解译精度,如何对这类观测曲线进行合理有效地平滑,是目前数据处理中的重点,也是难点.本文把各向异性扩散（anisotropic diffusion AD）引入到曲线平滑中,提出了各向异性扩散的曲线平滑方法,以平滑点梯度值的降函数作为扩散速度,在梯度较大的位置予以较小的扩散速度以保护特征,在梯度较小的位置予以较大的扩散速度平滑噪声引起的扰动.同时为了减小“飞点”(outliers)对平滑的干扰,采用局部单调扩散进行预处理.实验结果显示,本文提出的平滑方法有效地平滑了噪声与“飞点”的干扰,恢复了曲线的基本形态,有效地保持了曲线基本特征.     

Studies on turbulent diffusion processes and evaluation of diffusivity values from hydrodynamic observations in Corpus Christi Bay

The physical process of dispersion which can be attributed to turbulence (turbulent diffusion) or shear (shear-augmented diffusion) within the flow field is very important as it ultimately governs the distribution of constituents of interest within the environment. A series of diffusion experiments were conducted in Corpus Christi Bay, TX with the purpose of characterizing turbulent diffusion through dispersion coefficients or turbulent diffusivity, K_i (i=x, y, z) dependent on the degree of randomness or turbulence intensity, I.Measured with a boat-mounted acoustic doppler current profiler (ADCP), the Eulerian velocity time-series of fluid particles in random motion, u_i was used in the evaluation of the Eulerian time-scale of turbulence, T^E based on the velocity correlation function, R^E with T^E being related to the Lagrangian time-scale T^L through a scaling parameter, β(=T^L/T^E). Surface currents were obtained with high frequency (HF) Radar equipment deployed over the study area from which the horizontal velocity gradients were determined.Within the spatial scale of the experiment (1000 m), the observed low horizontal gradients (10⁻⁴ s⁻¹) allowed for the generation of velocity time-series from an ADCP mounted on a moving platform. A numerical scheme for evaluating turbulent diffusivity values was developed on the basis of Eulerian current measurements and calibrated through the statistics of an evolving dye patch for the scaling parameter β which in this scheme was found to be in the range 1–3.     

Evaluation and application of dialysis porewater samplers for microbiological studies at sediment-water interfaces

The equilibrium diffusion technique has become a valuable tool for ecological and biogeochemical studies in aquatic environments. In sediment ecosystems, changes in concentration of microbial metabolites with increasing depth can be determined dependably and reproducibly with this technique. Since the permeation characteristics of the membranes employed are crucial, selecting the proper membrane requires knowledge about its behavior under conditions which prevail in the natural environment. Thirteen polymer sheets were evaluated comparing permeation terms for biogeochemically relevant solutes, biodegradability, and mechanical strength. Cellulose-based dialysis membranes are most satisfactory when employed in low temperature anoxic environments. For this membrane, correction terms were calculated to account for diffusion losses during retrieval and sampling. Optimal incubation times can now be predicted from experimentally determined permeation coefficients for several porewater solutes. Dialysis porewater samplers (DPS) have been successfully applied during more than 100 independent experiments for the collection of interstitial water from surface sediments. DPS were used for water depths as deep as 290 meters.     

Efficient DEM‐based overland flow routing using integrated recursive algorithms

Two primary concerns in performing watershed overland flow routing are the numerical instability and computational efficiency. The stability of executing an explicit scheme has to be maintained by observing the Courant–Friedrich–Lewy criterion, which is adopted to confirm that the numerical marching speed is larger than the wave celerity. Moreover, there is another criterion of time step devised in previous studies to avoid back‐and‐forth refluxing between adjacent grids. The situation of refluxing usually occurs on flat regions. In light of this, the selection of a small time increment to honor both restrictions simultaneously is believed to decrease the computational efficiency in performing overland flow routing. This study aims at creating a robust algorithm to relax both restrictions. The proposed algorithm was first implemented on a one‐dimensional overland plane to evaluate the accuracy of the numerical result by comparing it with an analytical solution. Then, the algorithm was further applied to a watershed for 2D runoff simulations. The results show that the proposed integrated algorithm can provide an accurate runoff simulation and achieve satisfactory performance in terms of computational speed.     

基于镭同位素分布的黄海和东海垂直混合速率计算

在黄海和东海采样测定了水体中的镭同位素分布,用平流扩散模型描述镭同位素分布,最小二乘方法计算了垂直涡动扩散系数和上升流或下降流流速.结果给出北黄海中部、南黄海中部、浙江沿岸和台湾北部海域存在上升流,流速分别为0.46×10-3cm·s-1、0.17×10-3~1.39×10-3cm·s-1、2.02×10-3~3.04×10-3 cm·s-1和1.06×10-3~2.51×10-3 cm·s-1.北黄海中部和东海东北部存在下降流.流速分别为-2.30×10-3 cm·s-1和-0.61×10-3~-2.10×10-3 cm·s-1.计算同时给出的垂直涡动扩散系数为5.84~48.2 cm2·s-1,平均值为22.3 cm2·s-1.北黄海和浙江沿岸上升流流速与文献的结果一致;北黄海中部存在下降流与文献的结论一致.本研究结果与文献结果一致是对所建立的方法的肯定,也是对文献研究结果的支持.     

水动力条件对沉积物-水界面氧通量的影响

氧环境决定了水体沉积物中各种生命所需元素的最终归趋,沉积物-水界面是水相与沉积物相氧传递的重要场所,而水动力条件是影响沉积物-水界面氧传递的重要因素.选择三峡库区一级支流御临河为研究对象,根据长年监测数据建立实验室模型,采用声学多普勒流速测试仪及微电极测试系统构建了非侵入式涡度相关测试系统,探究了不同水动力条件对沉积物-水界面氧通量的影响.结果表明：水体静止状态下沉积物-水界面溶解氧浓度随时间的增加而减少,非静止状态下随时间的增加而增加;沉积物-水界面氧通量随水体流速的增加而增加.根据氧通量求解对应流速下垂直涡动扩散系数并进行线性拟合,当水体流速为0.01~0.14 m/s时,垂直涡动扩散系数与水体流速的相关性最好,此时沉积物-水界面氧通量的传递以涡动扩散为主导.     

3D P-WAVE VELOCITY STRUCTURE AT THE NORTHEASTERN MARGIN OF ORDOS BLOCK

Using the 7 100 absolute first arrivals of P waves and 91 513 relative P arrival times of 726 events at the northeastern margin of the Ordos block since 2009, the 3D fine structure of P wave velocity within the depth of 15km in the crust was inverted by the double difference seismic tomography method. The results show that there exist obvious high-speed continuous bodies in the northwest of the study area, and their lateral areas increase gradually with depth, while the velocity of east and south is relatively low. The velocity inhomogeneity exists and differs at different depths. The lateral differences of velocity are related to seismicity and faults. The 5~15km depth profile shows that earthquakes tend to occur in the area with relatively high velocity or high speed transition zones, which to some extent reflects the fragility of regional crustal media and the strong differential movement of faults in vertical and horizontal directions where the crust body is easy to absorb and store strain energy and generate major earthquakes. A "Y"-shape low-velocity channel is present in the lower crust around Liangcheng, corresponding to the NW-trending Heilaoyao-Shahukou fault set, which may reveal the migration path of the Late Tertiary-Quaternary basalt eruption. The Helingeer M6.2 earthquake in 1976 was related to the formation of the locking section of the thermal welding in this area. The three-dimensional fine structure of P wave velocity presented in this paper provides intuitive seismological evidence for physical and chemical properties of crustal media and the deep tectonic environment of earthquake preparation.     

Postglacial topographic evolution of glaciated valleys: a stochastic landscape evolution model

The retreat of valley glaciers has a dramatic effect on the stability of glaciated valleys and exerts a prolonged influence on the subsequent fluvial sediment transport regime. We have studied the evolution of an idealized glaciated valley during the period following retreat of ice using a numerical model. The model incorporates a stochastic process to represent deep‐seated landsliding, non‐linear diffusion to represent shallow landsliding and an approximation of the Bagnold relation to represent fluvial sediment transport. It was calibrated using field data from several recent surveys within British Columbia, Canada. We present ensemble model results and compare them with results from a deterministic linear‐diffusion model to show that explicit representation of large landslides is necessary to reproduce the morphology and channel network structure of a typical postglacial valley. Our model predicts a rapid rate of fluvial sediment transport following deglaciation with a subsequent gradual decline, similar to that inferred for Holocene time. We also describe how changes in the model parameters affect the estimated magnitude and duration of the paraglacial sediment pulse. Copyright © 2005 John Wiley & Sons, Ltd.