Diffusive dehydration and bubble resorption during open-system degassing of rhyolitic melts

Degassing through open paths such as bubble and/or fracture networks is considered to be the principal mode of degassing in silicic magmas. However, its detailed mechanisms remain unclear. To investigate the behavior of bubbles in a hypothetical open-system condition, we performed a series of vesiculation experiments on natural rhyolitic obsidian using a newly designed semipermeable cell, which artificially maintains a pressure difference between its inside and outside. The thick-wall cell maintains a constant volume within the sample chamber, while allowing water vapor to escape the cell during the experimental runs. The cells containing obsidian cores with ca. 0.66 wt.% initial water content were externally heated to 1000 °C for a period of 1–288 h. The run charges generally showed a zonal structure composed of two contrasting regions: a central region within which the bubbles were concentrated (bubble-rich core, BC) and a bubble-free melt region surrounding the BC (bubble-free margin, BFM). With increased heating duration, the thickness of the BFM increased via dissolution of the outermost bubbles in the BC. The water content was nearly uniform throughout the BC, whereas an outward-decreasing gradient was observed in the BFM. We found that diffusive dehydration occurred from the sample surface, and the bubbles were resorbed into the melt. Thus, the BFM–BC boundary moved inwards. These processes were modelled numerically, and the calculation results were in good agreement with the experimental data. If a “lifetime” of open paths is approximated as the relaxation time of a melt in a shallow volcanic environment, then the paths have to be pinched off quickly (1.2 h at maximum pressure difference between the open path and the melt) and thus the thickness of the bubble-free layer reaches at most ~ 0.1 mm. For the formation of bubble-free obsidian layers with a width of a few millimeters, which are often observed in natural obsidian flows, open paths should be maintained for at least a few hundred hours.     

A finite element solution for the fractional advection–dispersion equation

The fractional advection–dispersion equation (FADE) known as its non-local dispersion, has been proven to be a promising tool to simulate anomalous solute transport in groundwater. We present an unconditionally stable finite element (FEM) approach to solve the one-dimensional FADE based on the Caputo definition of the fractional derivative with considering its singularity at the boundaries. The stability and accuracy of the FEM solution is verified against the analytical solution, and the sensitivity of the FEM solution to the fractional order α and the skewness parameter β is analyzed. We find that the proposed numerical approach converge to the numerical solution of the advection–dispersion equation (ADE) as the fractional order α equals 2. The problem caused by using the first- or third-kind boundary with an integral-order derivative at the inlet is remedied by using the third-kind boundary with a fractional-order derivative there. The problems for concentration estimation at boundaries caused by the singularity of the fractional derivative can be solved by using the concept of transition probability conservation. The FEM solution of this study has smaller numerical dispersion than that of the FD solution by Meerschaert and Tadjeran (J Comput Appl Math 2004). For a given α, the spatial distribution of concentration exhibits a symmetric non-Fickian behavior when β = 0. The spatial distribution of concentration shows a Fickian behavior on the left-hand side of the spatial domain and a notable non-Fickian behavior on the right-hand side of the spatial domain when β = 1, whereas when β = −1 the spatial distribution of concentration is the opposite of that of β = 1. Finally, the numerical approach is applied to simulate the atrazine transport in a saturated soil column and the results indicat that the FEM solution of the FADE could better simulate the atrazine transport process than that of the ADE, especially at the tail of the breakthrough curves.     

辽宁地区中等地震活动增强与后续最大震级的关系探讨

本文定义了描述中等地震活动增强的三个参量E、T和R,通过对辽宁地区13次ML≥5.0地震的统计分析,分别给出了辽西-辽蒙交界、营-海-岫和辽南及两侧海域地区主震震级对E、T、R的回归方程并与实际震级进行对比分析,结果表明,中等地震活动增强参量E、T、R与未来主震震级间存在一定的定量关系,在震级预测上有一定指导意义。     

Study on forward modeling of through-casing resistivity logging

Through-casing resistivity logging (TCRL) is a new kind of logging method for formation evaluation,reservoir monitoring and saturation changes tracking. In this paper,the basic principle of TCRL is de-tailed,and the modified transmission line equation method (MTLEM) is first proposed in China. The MTLEM has been employed to simulate the responses of TCRL,on the basis of which,the effects of formation resistivity,cement and casing on the responses of TCRL are analyzed,and the signals in-duced in different formations are also calculated,which can help the design of TCRL tool and data processing and interpretation in China.     

剪胀性砂土本构模型的研究

对于中密砂和密砂,剪胀现象是一个比较显著的特征,相变线是描述剪胀现象的一条特征线,它比较容易确定。因此,以相变线作为状态参考线,提出并定义了一种状态参量,在已有弹塑性模型基础上,将状态参量引入剪胀方程和塑性硬化模量表达式中,建立了一个新的考虑砂土剪胀性的本构模型。该模型尤其适合于中密砂和密砂,对于松砂则退化为原弹塑性模型。模型形式较简单,包含9个材料参数,且比较容易确定。数值模拟结果与试验数据吻合较好。     

绳索取心钻杆的瞬间液相扩散连接

采用瞬间液相扩散连接方法焊接绳索取心钻杆。试验表明, 瞬间液相扩散连接接头拉伸时断在45MnMoB杆体侧, 面弯和背弯180°未开裂, 且焊缝内外表面无余高, 成形美观, 满足绳索取心钻杆焊接接头的技术要求。扭转试验中, 管体首先开裂, 焊缝未开裂, 扭转强度远远高于螺纹连接扭转强度。表明瞬间液相扩散连接技术是满足绳索取心钻杆焊接要求的一种先进方法。      

两淮地区煤炭资源地质勘探类型定量综合评价方法

安全、高效的现代化矿井建设,要求煤炭资源勘探速度快、精度高。通过总结安徽两淮地区近几十年来煤炭资源综合勘探经验,围绕“高精度、快速勘探”的主题,采用恰当的评价指标,对煤层稳定性和构造复杂程度进行定量评价,并提出了本地区煤炭资源地质勘探类型定量综合评价原则和方法,开发了定量评价软件系统。实际应用效果证明,该方法可行,结果可靠。      

The integrated risk assessment and zoning of meteorological disaster in Heilongjiang Province of China

This paper describes studies on the integrated risk assessment and zoning of meteorological disaster in Heilongjiang Province,in northeastern China,by using information-diffusion theory and cluster method with 35 years of summer temperature and precipitation data from 74 meteorological stations from 1971 through 2005.The information-diffusion theory has been used extensively in risk assessment,yet almost no one has done research about risk assessment by information-diffusion theory based on meteorological disaster standards.Some research results are as follows:the risk probability of low temperature and cold damage in the northern region is higher than that in the southern region;the risk probability of general low temperature and cold damage in the southwestern region is the highest;the risk probability of serious low temperature and cold damage in the northern region is the highest,followed by the central and southeast region;the high-risk region of arid disaster in Heilongjiang Province was primarily located in the southwestern,central,and southern parts of the province;the high-intensity arid disaster was located in the south-eastern region;the high-risk region of flood in Heilongjiang Province was primarily located from the southwest and then across the central part to the western part of Heilongjiang Province;the high-intensity flood disasters were located in almost every part of Heilongjiang Province.We can conclude from the integrated meteorological disaster risk zoning that the high-risk region of mete-orological disaster is primarily located in the southern and northern part of the province,the moderate-risk region is distributed in the central southern region and western region,the low-risk region is located in the eastern part,and the light-risk region is located in the central western part of Heilongjiang Province.     

M5 model tree based modelling of reference evapotranspiration

This paper investigates the potential of M5 model tree based regression approach to model daily reference evapotranspiration using climatic data of Davis station maintained by California irrigation Management Information System (CIMIS). Four inputs including solar radiation, average air temperature, average relative humidity, and average wind speed whereas reference evapotranspiration calculated using a relation provided by the CIMIS was used as output. To compare the performance of M5 model tree in predicting the reference evapotranspiration, FAO–56 Penman–Monteith equation and calibrated Hargreaves–Samani relation was used. A comparison of results suggests that M5 model tree approach works well in comparison to both FAO–56 and calibrated Hargreaves–Samani relations. To judge the generalization capability of M5 model tree approach, model created by using the Davis data set was tested with the datasets of four different sites. Results from this part of the study suggest that M5 model tree could successfully be employed in modeling the reference evapotranspiration. Further, sensitivity analysis with M5 model tree approach suggests the suitability of solar radiation, average air temperature, average relative humidity, and average wind speed as input parameters to model the reference evapotranspiration Copyright © 2009 John Wiley & Sons, Ltd.     

曲率模态小波法用于网壳结构损伤的识别和定位

工程结构损伤的识别与定位研究以往主要针对梁、框架等结构形式,根据大跨度空间结构杆件和节点繁多等特点,提出用曲率模态和小波混合方法对空间结构的损伤进行识别和定位.以跨度100 m的Schwedler网壳结构损伤前、后的曲率模态作为标识量,分别通过离散和连续小波变换,判断网壳结构有无损伤和损伤位置,统计了小波系数差与结构损伤的图形关系,计算了各种损伤工况下该方法判断损伤的准确程度.结果发现基于曲率模态和小波方法的大跨度网壳结构损伤定位精度很高,充分证明该方法对此类结构损伤定位具有有效性和实用性.