Near-source noise suppression of AMT by compressive sensing and mathematical morphology filtering

In deep mineral exploration, the acquisition of audio magnetotelluric (AMT) data is severely affected by ambient noise near the observation sites; This near-field noise restricts investigation depths. Mathematical morphological filtering (MMF) proved effective in suppressing large-scale strong and variably shaped noise, typically low-frequency noise, but can not deal with pulse noise of AMT data. We combine compressive sensing and MMF. First, we use MMF to suppress the large-scale strong ambient noise; second, we use the improved orthogonal match pursuit (IOMP) algorithm to remove the residual pulse noise. To remove the noise and protect the useful AMT signal, a redundant dictionary that matches with spikes and is insensitive to the useful signal is designed. Synthetic and field data from the Luzong field suggest that the proposed method suppresses the near-source noise and preserves the signal well; thus, better results are obtained that improve the output of either MMF or IOMP.     

2017年8月9日精河6.6级地震前LURR异常时空演化特征

基于加卸载响应比理论,对2017年8月9日新疆精河6.6级地震前加卸载响应比时空演化特征和异常区的孕震积分时程曲线进行分析。结果表明,在精河6.6级地震前,研究区出现了明显的加卸载响应比高值异常区域,异常区呈现环状特征;而且,异常区域随时间逐渐向震源附近迁移,异常演化过程中异常区域内部和周围发生了6次5级以上地震,其中,2次6级地震发生在异常高值减弱过程中,异常峰值到发震的时间分别为11、22个月,大于以往研究的理论发震时间,这可能是受构造应力区域性特征影响的缘故。     

Analytical approach for sheet flow transport in purely acceleration-skewed oscillatory flow

An instantaneous analytical approach is developed to predict sheet flow transport in purely acceleration-skewed oscillatory flow. The approach is derived from exponential approximations of velocity and concentration profiles above a mobile seabed, and it particularly considers factors of phase lead; phase lag (i.e. phase residual and phase shift); acceleration modification; and asymmetries in shear stress, roughness height, and boundary layer development. The approach can predict net boundary layer flow above a mobile seabed, and can revert to the classical bedload model. Instantaneous and net sediment transport rates are studied using the approach. The instantaneous sediment transport rate in an onshore flow stage can be approximated by a power function of velocity in which the exponent is confirmed to range between 1 and 5 with a decrease in the phase residual. The net sediment transport rate predicted using the approach is validated using a considerable amount of measured data, and compared with existing instantaneous and half-period type models that consider the phase lag or acceleration modification. For the net sediment transport rate in purely acceleration-skewed oscillatory flow, the phase residual is less important than the acceleration-skewed boundary layer difference between onshore and offshore acceleration stages.     

Conductivity model for pyrite-bearing laminated and dispersed shaly sands based on a differential equation and the generalized Archie equation

The conductance of pyrite-bearing laminated and dispersed shaly sands is not well understood and resistivity models for pyrite-bearing shaly sands are nonexistent. Thus, we first synthesize clean pyrite-matrix samples, and quartz-matrix samples with variable laminated shale, dispersed shale, and pyrite content and then perform petrophysics experiments to assess the effect of pyrite content on the conductivity of pyrite-bearing shaly sands. Second, based on the differences in conductivity and conduction pathways and geometries because of the variable composition of the pyrite-bearing laminated and dispersed shaly sands, we divide the shaly sands into their components, i.e., laminated shale, quartz grains, pyrite grains, hydrocarbon, dispersed shale, microscopic capillary water, and mobile water. A generalized resistivity model is proposed to describe the conductivity of pyrite-bearing laminated and dispersed shaly sands, based on the combined conductivity differential equation and generalized Archie equation. In the generalized resistivity model, the conductivity differential equation is used to describe the conductivity of dispersed inclusions in a host, whereas the generalized Archie equation is used to describe the conductivity of two conducting phases. Moreover, parallel conductance theory is used to describe the conductivity of dispersed shaly sands and laminated shale. Theoretical analysis suggests that the proposed model satisfies the physical constraints and the model and experimental results agree. The resistivity and resistivity index of shaly sands decrease with increasing conductivity and pyrite. Finally, the accuracy of the resistivity model is assessed based on experimental data from 46 synthetic core samples with different oil saturation. The model can describe the conductivity of clean pyrite-matrix samples, and quartz-matrix samples with different volumes of laminated shale, dispersed shale, and pyrite. An accurate saturation model of pyrite-bearing laminated and dispersed shaly sands is thus obtained and the log data interpretation in complex shaly sands can improve with the proposed model.     

Study and application of PS-wave pre-stack migration in HTI media and an anisotropic correction method

Anisotropy correction is necessary during the processing of converted PSwave seismic data to achieve accurate structural imaging, reservoir prediction, and fracture detection. To effectively eliminate the adverse effects of S-wave splitting and to improve PSwave imaging quality, we tested methods for pre-stack migration imaging and anisotropic correction of PS-wave data. We based this on the propagation rules of seismic waves in a horizontal transverse isotropy medium, which is a fractured medium model that reflects likely subsurface conditions in the field. We used the radial (R) and transverse (T) components of PS-wave data to separate the fast and slow S-wave components, after which their propagation moveout was effectively extracted. Meanwhile, corrections for the energies and propagation moveouts of the R and T components were implemented using mathematical rotation. The PS-wave imaging quality was distinctly improved, and we demonstrated the reliability of our methods through numerical simulations. Applying our methods to three-dimensional and three-component seismic field data from the Xinchang-Hexingchang region of the Western Sichuan Depression in China, we obtained high-quality seismic imaging with continuous reflection wave groups, distinct structural features, and specific stratigraphic contact relationships. This study provides an effective and reliable approach for data processing that will improve the exploration of complex, hidden lithologic gas reservoirs.     

An acoustic-wave theory for casing bond evaluation using viscoelastic slip boundary modeling

Slip boundary condition is commonly utilized to model elastic wave propagation through layered earth media. The same approach is used here to characterize acoustic wave propagation along a cased borehole with various cement bond conditions. By modeling the cement layer between casing and formation as a viscoelastic slip interface with complex coupling rigidity parameters, one can not only reduce the complexity in the classical elastic wave modeling of the problem, but also efficiently model various complicated wave phenomena that are difficult for the existing modeling. More specifically, the new theory can well describe the effect of the cement bond condition change and the location of the change (i.e., whether it is in the first interface between casing and cement, or the second interface between cement and formation) on the acoustic waves, demonstrating the good modeling capability and predicting power. Application of the theory to field data shows that the theory can correctly model the acoustic wave characteristics and interpret the cement bond condition, thus providing a useful fundament theory for casing bond evaluation using acoustic logging.     

Field survey and analysis of hillslopes on tower karst in Guilin,southern China

Limestone dissolution in tropical and subtropical humid southern China created residual hills with steep slopes, a landform that is referred to as tower karst. Two types of tower karst landform feature, fenglin or peak forest (isolated towers) and fengcong or peak cluster (linked‐base towers), were identified in Guilin. Previous studies proposed two hypotheses regarding their origin and evolution. One is the sequential evolution model from peak cluster to peak forest. The other is a parallel development model, which postulates that both peak cluster and peak forest have developed simultaneously. Through detailed field survey and analysis of slope forms on tower karst in Guilin, it was found that the mean slope angle of the towers is very high (62·4°) and ranges from 60° to 75°. There is no significant difference in mean slope angle and slope angle distribution between towers in the peak cluster basin and peak forest floodplain areas. Mean slope angle increases with intensified fluvial dissection. Three levels of caves in the towers of the peak forest in Guilin were identified in previous research. The isolated towers of the peak forest as well as scattered residuals of peak cluster are generally distributed in the centre of the Guilin syncline. Favourable circumstances of allogenic water concentration indicate that development of the peak forest resulted from the combined effects of subcutaneous and subterranean dissolution as well as subsequent collapse and recession by fluvial erosion after uplifting. By contrast, peak clusters generally occur on the limbs of the syncline or at the periphery of the Guilin basin with relatively higher elevations. The thick vadose zone and predominantly vertical flow suggests that peak clusters are mainly formed by the combination of intensive uplifting and the enhancement of original dolines. The evidence of slope survey and slope analysis suggests that both isolated towers and linked‐base towers developed simultaneously but by different mechanisms of formation and different combinations of development processes. Copyright © 2000 John Wiley & Sons, Ltd.     

东线南水北调对黄浦江水质的影响

本文对长江枯季(12—3月)径流的特点、变化规律、河口盐水入侵以及对黄浦江水质的影响等进行了初步分析。长江水量虽然丰富,但枯季径流并不大,枯季径流一般占年径流总量的15%左右。频率50%的平水年各月平均流量都在15000秒立米以下,频率75%的一般少水年枯季1、2月份的月平均流量不足10000秒立米。河口盐度变化随上游来水多少而异,大通站流量在10000秒立米时河口盐度变化敏感,流量再减少河口盐度将急剧增加,便会严重污染河口和黄浦江水质,影响上海市和河口地区的工农业和城市生活用水。因此,在规划长江流域用水时,应使大通站流量在15000秒立米以上,若有困难,至少不能低于10000秒立米。     

帕米尔高原Rushan-Pshart古特提斯洋或中特提斯洋争论的岩石地球化学证据

帕米尔高原上广泛分布的加里东期火山岩,印支期火山岩与燕山期火山岩被认为是块体依次向北俯冲拼贴到欧亚板块上的产物。特别是燕山期火山岩被认为是Rushan-Pshart中特提斯洋闭合的产物。但近年来地层古生物、火山岩证据不断表明Rushan-Pshart洋闭合时间在晚三叠一早侏罗世,由此限定Rushan-Pshart古特提斯洋性质,而帕米尔高原上广泛分布的燕山期火山岩是更南部的Shyok中特提斯洋闭合的产物。目前,中国境内Rushan-Pshart缝合带属性的研究工作展开较晚,研究程度较低。我们对塔什库尔干明铁盖沟一线燕山期火山岩带展开的工作发现零星分布的印支期花岗岩。印支期花岗岩错石U・Pb年龄显示岩体侵位时间在201 Ma左右。全岩主量元素特征表明岩石为I型高钾钙碱性闪长花岗岩;稀土元素在球粒陨石标准化图解中呈轻稀土元素相对富集,重稀土元素相对亏损的右倾海鸥型。明铁盖岩体的微量元素显示大离子亲石元素明显富集,而高场强元素明显亏损。稀土微量元素特征倾向花岗岩属性为岛弧型花岗岩,Sr-Nd-Hf同位素比值特征显示岩体形成于下地壳部分熔融环境。岩体的地化特征表明岩体形成的构造环境为板块汇聚的洋壳俯冲阶段,结合区域地质特征,我们将花岗岩体归为Rushan-Pshart古特提斯洋壳俯冲消减的产物。Rushan-Pshart缝合带传统上认为是中特提斯带,近年的研究进展认为其为古特提斯缝合带,本文的工作支持这种观点。Rushan-Pshart古特提斯缝合带的确立对帕米尔高原与青藏高原主体的块体对比提供了可信的对比方案,并对青藏高原新生代陆内变形方式的争论提供了可靠的证据。     

Weakening laws of slip zone soils during wetting–drying cycles based on fractal theory: a case study in the Three Gorges Reservoir (China)

Acta Geotechnica - Slip zone soils play a crucial role in the evolution of landslides. In this paper, slip zone soils of the Huangtupo Landslide in the Three Gorges Reservoir area, which is an...