Modeling of ice phenomena in the mouth of the Vistula River

The mouth of the Vistula River, which is a river outlet located in tideless area, is analyzed. The Vistula River mouth is a man-made, artificial channel which was built in the 19th century in order to prevent the formation of ice jams in the natural river delta. Since the artificial river outlet was constructed, no severe ice-related flood risk situations have ever occurred. However, periodic ice-related phenomena still have an impact on the river operation. In the paper, ice processes in the natural river delta are presented first to refer to the historical jams observed in the Vistula delta. Next, the calibrated mathematical model was applied to perform a series of simulations in the Vistula River mouth for winter storm condition to determine the effects of ice on the water level in the Vistula River and ice jam potential of the river outlet.     

The application of two‐dimensional hydrodynamic models for underwater archaeological finds from the Vistula River in Warsaw,Poland

Low water levels occurring on the Warsaw section of the Vistula River during droughts in 2012 and 2015 stirred considerable interest by exposing parts of the riverbed. The river's low flow discharge created favourable conditions for archaeologists who, from the river's bottom, managed to salvage sculptures, architectural artefacts, and military accessories dating back to the middle of the 17th century. Literature shows that archaeological finds at the bottom of a major lowland river are unique and are attributed to very specific hydrological and geological conditions. Embankments and works meant to regulate the channel of the Vistula River in Warsaw narrowed the channel and accelerated the erosion of the river's bottom. A comparison of rating curves created between 1919 and 2015 shows that the river's bottom had lowered by 225 cm. The diameter of sediment particles that could be moved at maximum shear stress has been calculated using the Meyer–Peter and Muller formula. According to the calculated shear stress corresponding to the highest historical flood of 1844 particles smaller than 0.05 m were transported down the river, and it is true that artefacts of that size had not been collected at the site. Taking into consideration the very specific geological and hydrological conditions, we believe that finding so many 17th century artefacts with a well‐documented history at a site at the bottom of major lowland river with intensive channel processes is indeed unique. Additionally, it is the first time that results of hydrodynamic modelling have been successfully used for archaeological survey to define the critical shear stress and the smallest diameter of artefacts, which could remain stable on a river's bottom.     

煤巷小构造Rayleigh型槽波超前探测数值模拟

对煤巷小构造地震波场进行了数值模拟研究,分析了层状煤层中地震波的传播特征.研究表明:(1)在煤巷迎头前方煤层内以纵波震源激发的Rayleigh型槽波相对于体波能量较强,波列较长,波速较低.(2)沿煤层传播的Rayleigh型槽波在小构造面上产生Rayleigh型槽波反射波,反射Rayleigh型槽波垂直分量相对于水平分量能量较强.沿煤层反向传播的反射Rayleigh型槽波在煤巷迎头面上转换为沿煤巷底板传播的Rayleigh面波.沿煤巷底板可以接收到能量较强的反射Rayleigh型槽波产生的Rayleigh面波,其可以作为超前探测小构造面的特征波.在地震记录上反射Rayleigh型槽波产生的Rayleigh面波波至最迟,在时间域与其他波列时间间隔较大,其垂直分量能量相对于水平分量较强,在地震记录上容易识别.(3)在相同的地质条件下应用反射地震超前探测方法,标志煤巷迎头前方存在小构造面的反射地震波能量较弱,受煤巷顶、底板界面和采煤迎头面的强反射波干扰,在地震记录中难以识别.     

CROSS-BOREHOLE OBSERVATION OF MODE CONVERSION FROM BOREHOLE STONELEY WAVES TO CHANNEL WAVES AT A COAL LAYER1

Conversion of borehole Stoneley waves to channel waves was observed in data from a seismic cross-borehole experiment conducted between wellbores penetrating a thin coal layer at 2022 m depth, near Rifle, Colorado. Traveltime moveout observations show that borehole Stoneley waves underwent partial conversion to channel waves at the coal layer. The channel waves were detected directly in an adjacent borehole 35 m away at receiver positions within the coal. Stoneley waves, subsequently produced by partial conversion of channel waves, were also detected at receiver positions located up to 50 m above and below the coal layer in the adjacent borehole. We infer the channel wave to be the first-higher Rayleigh mode by comparing the observed group velocity with theoretically derived dispersion curves. Identifying the conversion between borehole and stratigraphically guided waves is significant because coal penetrated by multiple wells may be detected without placing a transmitter or receiver in the coal itself.     

0.9 m薄煤层SH型槽波频散特征及波形模式

在0.9 m薄煤层中使用放炮方法做微震震源的条件,通过对所采集到的薄煤层槽波信号进行分析,发现薄煤层槽波在频域中存在高频和低频两个独立并且不连续的波段.其中高频区大约以2000 Hz为中心,低频域的中心频率约为490 Hz,并且高频域和低频域的能量差异不显著.通过时频分析,可以清晰地看到高、低两个频域几乎在同一时刻触发,并且其小波相关系数在这两个域中的分布规律表现出一定的相似性.通过对0.9 m薄煤层槽波频散曲线的理论分析可知,现场观测到的Airy震相的频率及速度和其理论值较为接近.震源置于煤层中心,且炸药能量对顶底板的扰动,对第二阶对称波形模式下槽波高频部分的形成起着关键作用.在这一对称高阶波形模式下槽波的波速基本上和煤层顶底板中S波的波速一致.由于这两个触发的波形模式在时间域中具有相似的特征,建议在高频域和低频域同时发育较好的薄煤层槽波勘探中,可以采用在同一时间域中高低频相结合的方法提高利用槽波勘探分析的效果.     

Numerical simulations of full-wave fields and analysis of channel wave characteristics in 3-D coal mine roadway models

Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling twodimensional wave fields and therefore cannot accurately simulate three-dimensional (3-D) full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity–stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity V_x, V_y, and V_z for the same node in 3-D staggered-grid finite difference models by calculating the average value of V_y, and V_z of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways. Instead, the Rayleigh channel wave can be used for the advance detection of geological structures in coal mine roadways.     

HTI煤层介质槽波波场与频散特征初步研究

煤层内裂隙较为发育,具有明显的各向异性.目前槽波理论研究以各向同性介质为主,对HTI介质中槽波及其频散性质研究很少.本文以弱各向异性、含垂直裂隙HTI煤层介质为研究对象,研究了HTI煤层介质中的三维槽波波场,采用交错网格高阶有限差分法模拟槽波,推导了三层水平层状HTI煤层介质的Love型槽波理论频散公式和振幅深度分布,分析了HTI各弹性参数对频散曲线的影响.HTI介质和各向同性介质基阶Love槽波频散曲线差异较小,高阶较大;煤厚主要影响Airy相频率,而Airy相速度不变;煤层vs对Airy相速度影响很大;煤层γ对基阶Love槽波影响很小,高阶稍大.各波偏振方向不再与波的传播方向平行或垂直,而是呈一定夹角.利用基阶Love槽波频散曲线推测裂隙发育较为困难,可利用高阶频散曲线.     

Dispersion features of transmitted channel waves and inversion of coal seam thickness

In-seam seismic survey currently is a hot geophysical exploration technology used for the prediction of coal seam thickness in China. Many studies have investigated the relationship between the group velocity of channel wave at certain frequency and the actual thickness of exposed coal beds. But these results are based on statistics and not universally applicable to predict the thickness of coal seams. In this study, we first theoretically analyzed the relationship between the depth and energy distribution of multi-order Love-type channel waves and found that when the channel wave wavelength is smaller than the thickness of the coal seam, the energy is more concentrated, while when the wavelength is greater than the thickness, the energy reduces linearly. We then utilized the numerical simulation technology to obtain the signal of the simulated Love-type channel wave, analyzed its frequency dispersion, and calculated the theoretical dispersion curves. The results showed that the dispersion characteristics of the channel wave are closely related to the thickness of coal seam, and the shear wave velocity of the coal seam and its surrounding rocks. In addition, we for the first time realized the joint inversion of multi-order Love-type channel waves based on the genetic algorithm and inversely calculated the velocities of shear wave in both coal seam and its surrounding rocks and the thickness of the coal seam. In addition, we found the group velocity dispersion curve of the single-channel transmitted channel wave using the time–frequency analysis and obtained the phase velocity dispersion curve based on the mathematical relationship between the group and phase velocities. Moreover, we employed the phase velocity dispersion curve to complete the inversion of the above method and obtain the predicted coal seam thickness. By comparing the geological sketch of the coal mining face, we found that the predicted coal seam thickness is in good agreement with the actual thickness. Overall, adopting the channel wave inversion method that creatively uses the complete dispersion curve can obtain the shear wave velocities of the coal and its surrounding rocks, and analyzing the depth of the abruptly changed shear wave velocity can accurately obtain the thickness of the coal seam. Therefore, our study proved that this inversion method is feasible to be used in both simulation experiments and actual detection.     

地震槽波的数学-物理模拟初探

针对地震槽波在低速层的传播特性,开展了煤层内地震槽波勘探的数值模拟和物理模拟研究的初探工作.在数值模拟研究方面,采用交错网格有限差分法对煤层中的地震槽波进行三分量全波场模拟.基于波场快照和人工合成地震记录研究了不同模型中的波场特征和各种波型的传播规律.在物理模拟方面,通过选用不同配比的环氧树脂和硅橡胶类材料构建地震槽波物理模型,利用透射法和反射法观测系统获得了清晰的地震槽波记录以研究槽波的地震学特征.研究表明,在煤层内槽波的地震波场中,Love型槽波的能量小于Rayleigh型槽波的SV分量,大于Rayleigh型槽波的SH分量.相对于Love型槽波和Rayleigh型槽波的SH分量,Rayleigh型槽波的SV分量在围岩中的泄露能量较强.在煤层界面附近的围岩中,地震波仍以槽波形式传播,随着距离的增加能量逐渐衰减.随着煤层变薄,煤层槽波主频向高频方向移动,频散现象增强,传播速度增大.     

槽波层析成像方法在煤田勘探中的应用——以河南义马矿区为例

在地下煤田的开发中,工作面内的小构造、异常体、煤层厚度变化等是需要解决的关键问题,而槽波探测则为这些问题的解决提供了重要的物探方法.本文对河南义马矿区11061工作面进行槽波透射法测量,在巷道显示的煤层厚度变化为1.5~8 m,从理论频散曲线分析速度与厚度关系,确定125 Hz频率槽波主要用于观测厚度约为2~5 m的煤层厚度变化;有效提取了684个频散曲线,并分别拾取了125 Hz时槽波群速度与走时,采用走时层析成像方法获得工作面内煤层速度、厚度以及高应力区分布特征,回采验证了结果的正确性.     

Adsorbable Organic Halogens (AOX) in Polish Rivers – Levels and Changes

This paper presents results of a study on the concentration of adsorbable organic halogens (AOX) in water samples of the Oder River and its tributaries collected from 1998 to 2002. The Vistula River water for AOX was analysed in the period 2001–2002. The investigated waters of the Vistula near the mouth (Kiezmark) were more contaminated with AOX than the Oder near the mouth (Widuchowa, Krajnik Dolny). In general, AOX levels in the Polish rivers are comparable to the values found in the European rivers, e.g. the Rhine and the Ton. The research has shown that AOX were correlated with pH, total suspended solids, and the riverine flow. No relation between dissolved oxygen and redox potential and AOX was found out. During the flood period AOX, inorganic halogens (TIX), and mineral substances (sulfates) were analysed.     

Trace elements in alluvia of the upper Vistula as indicators of palaeohydrology

In the Vistula valley near Góra, at the northern foreland of the Carpathians Mountains, there are young alluvia deposited in the historical period. These alluvia are difficult to date precisely because of the lack of organic interbedded material. In the Vistula drainage basin upstream of Góra since the second half of 18th century small centres of metal, textile, motor, and petrochemical industries have been developed. Concentrations of trace elements such as: Cr, Mn, Co, Ni, Cu, Zn, Cd, Pb, and of some macroelements were used as indicators of the age of different alluvial fills. The main part of the Vistula valley floor was built up by sediment accumulation before the development of industry but after the medieval phase of rural colonization of the Carpathians. Deposits filling the abandoned channels dissecting the valley floor have an increased content of some trace elements and deposits date from the second half of the 18th century, when the first foundries and blacksmith's shops were developed. A drastic increase of concentration of the most abundant trace elements in the youngest alluvial fill indicates that the fill dates from the last decade of rapid industrial development.     

Numerical study of seismic scattering and waveguide excitation in faulted coal seams

Finite‐difference P‐SV simulations of seismic scattering characteristics of faulted coal‐seam models have been undertaken for near‐surface P‐ and S‐wave sources in an attempt to understand the efficiency of body‐wave to channel‐wave mode conversion and how it depends on the elastic parameters of the structure. The synthetic seismograms clearly show the groups of channel waves generated at the fault: one by the downgoing P‐wave and the other by the downgoing S‐wave. These modes travel horizontally in the seam at velocities less than the S‐wavespeed of the rock. A strong Airy phase is generated for the fundamental mode. The velocity contrast between the coal and the host rock is a more important parameter than the density contrast in controlling the amplitude of the channel waves. The optimal coupling from body‐wave energy to channel‐wave energy occurs at a velocity contrast of 1.5. Strong guided waves are produced by the incident S‐sources for source angles of 75° to 90° (close to the near‐side face of the fault). As the fault throw increases, the amplitude of the channel wave also increases. The presence of a lower‐velocity clay layer within the coal‐seam sequence affects the waveguiding characteristics. The displacement amplitude distribution is shifted more towards the lower‐wavespeed layer. The presence of a ‘washout’ zone or a brecciated zone surrounding the fault also results in greater forward scattering and channel‐wave capture by the coal seam.     

含小断层煤层Rayleigh型槽波波场和频散分析

煤层隐伏小断层在煤田地质勘探期间能否查明,是影响煤矿安全生产的主要地质因素.煤层中传播的地震槽波,非常适用于探测煤田异常构造.本文利用谱元法模拟计算了含断层模型Rayleigh型地震槽波的产生、传播过程,对于直达槽波、反射槽波、透射槽波进行了波场特征分析,频散特征分析,频谱特征分析.根据不同小落差断层模型（垂直断距分别为1/4、1/2、3/4、1个煤厚;断层倾角分别为90°、60°、45°）,分析了Rayleigh型槽波传播过程中断层对反射槽波、透射槽波能量、频谱的影响效应.基于以上分析,对于小断层模型中的Rayleigh型地震槽波取得了一些规律性认识,并对Rayleigh型地震槽波的勘探应用做了一些探讨.     

Effects of heated water discharge on the evaporation from a river surface

The investigations concerned the effects of thermal waters on the river evaporation process and the development of formulae, based on standard hydrological and meteorological data, for determination of evaporation from the Vistula river.     

Modified image algorithm to simulate seismic channel waves in 3D tunnel model with rugged free surfaces

The existence of rugged free‐surface three‐dimensional tunnel conditions in the coal seams, caused either by geological or mining processes, will inevitably influence wave propagation characteristics when the seismic waves go through the coal mines. Thus, a modified image algorithm has been developed to account for seismic channel waves propagating through this complicated topography with irregular free surfaces. Moreover, the seismic channel waves commonly exhibit damped and dispersive signatures, which is not only because of their own unique sandwich geometry of rock–coal–rock but also because of the viscoelastic behavior of coal. Considering the complexity of programming in three‐dimensional tunnel models with rugged free surfaces, an optimized vacuum grid search algorithm, enabling to model highly irregular topography and to compute efficiently, is also proposed when using high‐order staggered finite‐difference scheme to simulate seismic channel wave propagations in viscoelastic media. The numerical simulations are implemented to investigate the accuracy and stability of the method and the impact of coal's viscoelastic behavior on seismic channel wave propagation characteristics. The results indicate that the automatic vacuum grid search algorithm can be easily merged into high‐order staggered finite‐difference scheme, which can efficiently be applied to calculate three‐dimensional tunnel models with rugged free surfaces in the viscoelastic media. The simulation also suggests that the occurrence of a three‐dimensional tunnel with free surfaces has a remarkable influence on the seismic channel wave propagation characteristics and elastic energy distribution.     

THE PROPAGATION OF ATTENUATED SH CHANNEL WAVES*

The problem of detecting discontinuities which interrupt coal seams is of great importance to the coal mining industry. One possible method of detecting such discontinuities is that of in-seam seismology where both source and detectors are placed under- ground on the coal face. In this paper the propagation of SH waves from a line source in the seam is investigated. There exists a particular set of waves–channel waves–, confined to the coal seam. These waves are dispersive and have an associated Airy phase whose geometrical attenuation is least. However, if absorption of energy within the coal is included then the Airy phase is the dominant wavegroup only for distances less than a certain maximum. If a detection criterion proposed by Dresen and Freystätter applies, then there is a lower limit to the frequency that should be used for exploration. This requirement may be in conflict with attenuation considerations.     

Storage of sediment‐associated heavy metals along the channelized Odra River,Poland

The importance of long‐term storage of heavy metals in groyne fields, functioning over 150 years, is investigated for the River Odra (Oder), western Poland. Construction of groynes along the Odra preceded rapid development of heavy industrialization in the largest coal mine districts in Poland and the Czech Republic that resulted in persistent riverine pollution. The 187 km long Middle Odra reach was repeatedly channelized from the first half of the eighteenth century to the turn of the twentieth century, during which time partially filled groyne fields were dissected by new bank lines and groyne systems, with older groyne fields partially keyed into the floodplain. Consequently, concentrations of zinc, lead, cadmium, and copper within historically deposited groyne field sediments exceed local geochemical background levels by more than 60, 40, 15 and 10 times, respectively. Sediments contaminated with heavy metals occur within three distinctive geomorphic zones: zone I is up to 250 m wide and furthest from the present channel, comprising decimeter‐thick polluted sediments, overlying eighteenth century sand and gravel bars; zone II represents the former nineteenth century groyne fields, with widths between 10 and 100 m, filled with as much as 3 m of polluted sediments; zone III represents the twentieth century groyne fields, which are several to a dozen metres wide and filled with polluted sediments averaging depths of more than 2 m. This investigation indicates that large and extensive sediment quantities of moderately polluted sediments are stored immediately along the banks of the River Odra. These sediments could be a significant secondary pollution source and therefore careful maintenance of contemporary bank protection structures is required. Copyright © 2009 John Wiley & Sons, Ltd.     

EFFECT OF COAL MINING ON RIVER SEDIMENTATION-A Case Study From The Shengfu-Dongsheng Coal FieldIn The Loess Plateau, North China

IINTRODUCTIONTheWulanmulunheRiver(theupperreachoftheKuyeheRiver)isatributaryoftheYellowRiver(Fig.I).Ithasawatershedareaof3839km2upstreamfromtheoutletgaugingstationatWangdaohengta(ithasbeensuccessfollyoperatedsince1958).ThisisafieldinthefragileecologicalareaoftheLoessPlateau.AnnualrainfallintheWulanmulenheRiverbasinrangesfrom124mmto640rumwithanaverageof343mm.Themaincharacteristicsofrainfallarehighintensityandspatialandtemporalheterogenety.Nearly70%ofthetotalannualrainfallfallsinthefl…     

Temporal and spatial variability of extreme river flow quantiles in the Upper Vistula River basin,Poland

Temporal and spatial variability in extreme quantile anomalies of seasonal and annual maximum river flows was studied for 41 gauging stations at rivers in the Upper Vistula River basin, Poland. Using the quantile perturbation method, the temporal variability in anomalies was analysed. Interdecadal oscillating components were extracted from the series of anomalies using the Hilbert‐Huang transform method. Period length, part of variance of each component, and part of unexplained variance were assessed. Results show an oscillating pattern in the temporal occurrence of extreme flow quantiles with clusters of high values in the 1960–1970s and since the late 1990s and of low values in the 1980s and at the beginning of the 1990s. The anomalies show a high variability on the right bank of the Upper Vistula River basin during the summer season with the highest values in catchments located in the western and south‐western parts of the basin. River flow extreme quantiles were found to be associated with large‐scale climatic variables from the regions of the North Atlantic Ocean, Scandinavia, Eastern Europe, Asia, and, to a lesser extent, the Pacific Ocean. Similarities between temporal variability of river flows and climatic factors were revealed. Results of the study are important for flood frequency analysis because a long observation period is necessary to capture clusters of high and low river flows.