用小波变换方法提取地下流体观测异常信息

首先对小波分析方法作了简单介绍,然后运用该方法,对1992～2002年上海地区9口水氡观测井的资料进行了6阶小波变换,提取的地震异常信息与实际地震有较好的对应关系;对1983年山东荷泽5.9级地震和1998河南清丰4.0级地震前聊古1井的氡观测资料进行了小波分析,发现1～6阶小波系数均出现异常;对1999年山东苍山5.2级地震前枣庄十里泉井的水汞观测资料进行了小波分析,发现1～6阶小波系数也出现异常.因此,运用小波分析方法可以有效地提取地震前兆异常信息.     

EWT-Elman组合模型短期电离层TEC预报

针对电离层总电子含量(TEC)非线性、高噪声的特点,建立基于经验小波变换(EWT)和Elman神经网络的短期电离层组合预报模型.运用该模型对不同地磁环境的电离层TEC时间序列进行建模预报,结果表明,EWT-Elman组合模型可反映电离层TEC的变化特征,地磁平静期预测平均相对精度为93％,均方根误差为1.04 TECu...     

Wavelet analyses of neural networks based river discharge decomposition

The problem of discharge forecasting using precipitation as input is still very active in Hydrology, and has a plethora of approaches to its solution. But, when the objective is to simulate discharge values without considering the phenomenology behind the processes involved, Artificial Neural Networks, ANN give good results. However, the question of how the black box internally solve this problem remains open. In this research, the classical rainfall-runoff problem is approached considering that the total discharge is a sum of components of the hydrological system, which from the ANN perspective is translated to the sum of three signals related to the fast, middle and slow flow. Thus, the present study has two aims (a) to study the time-frequency representation of discharge by an ANN hydrologic model and (b) to study the capabilities of ANN to additively decompose total river discharge. This study adds knowledge to the open problem of the physical interpretability of black-box models, which remains very limited. The results show that total discharge is adequately simulated in the time frequency domain, although less power spectrum is evident during the rainy seasons in the ANN model, due to fast flow underestimation. The wavelet spectrum of discharge represents well the slow, middle and fast flow components of the system with transit times of 256, 12–64 and 2–12 days, respectively. Interestingly, these transit times are remarkably similar to those of the soil water reservoirs of the studied system, a small headwater catchment in the tropical Andes. This result needs further research because it opens the possibility of determining MMT on a fraction of the cost of isotopic based methods. The cross-power spectrum indicates that the error in the simulated discharge is more related to the misrepresentation of the fast and the middle flow components, despite limitations in the recharge period of the slow flow component. With respect to the representation of individual signals of the slow, middle and fast flows components, the three neurons were uncapable to individually represent such flows. However, the combination of pairs of these signals resemble the dynamics and the spectral content of the aforementioned flows signals. These results show some evidence that signal processing techniques may be used to infer information about the hydrological functioning of a basin.     

A coupled 3‐dimensional bonded discrete element and lattice Boltzmann method for fluid‐solid coupling in cohesive geomaterials

This paper presents a 3D bonded discrete element and lattice Boltzmann method for resolving the fluid‐solid interaction involving complicated fluid‐particle coupling in geomaterials. In the coupled technique, the solid material is treated as an assembly of bonded and/or granular particles. A bond model accounting for strain softening in normal contact is incorporated into the discrete element method to simulate the mechanical behaviour of geomaterials, whilst the fluid flow is solved by the lattice Boltzmann method based on kinetic theory and statistical mechanics. To provide a bridge between theory and application, a 3D algorithm of immersed moving boundary scheme was proposed for resolving fluid‐particle interaction. To demonstrate the applicability and accuracy of this coupled method, a benchmark called quicksand, in which particles become fluidised under the driving of upward fluid flow, is first carried out. The critical hydraulic gradient obtained from the numerical results matches the theoretical value. Then, numerical investigation of the performance of granular filters generated according to the well‐acknowledged design criteria is given. It is found that the proposed 3D technique is promising, and the instantaneous migration of the protected soils can be readily observed. Numerical results prove that the filters which comply with the design criteria can effectively alleviate or eliminate the appearance of particle erosion in dams.     

Canopy cover effects on local soil water dynamics in a tropical agroforestry system: Evaporation drives soil water isotopic enrichment

Despite the widely held assumption that trees negatively affect the local water budget in densely planted tree plantations, we still lack a clear understanding of the underlying processes by which canopy cover influences local soil water dynamics in more open, humid tropical ecosystems. In this study, we propose a new conceptual model that uses a combination of stable isotope and soil moisture measurements throughout the soil profile to assess potential mechanisms by which evaporation (of surface soil water and of canopy‐intercepted rainfall) affects the relationship between surface soil water isotopic enrichment (lc‐excess) and soil water content. Our conceptual model was derived from soil water data collected under deciduous and evergreen plants in a shade grown coffee agroforestry system in Costa Rica. Reduced soil moisture under shade trees during the “drier” season, coinciding when these trees were defoliated, was largely the result of increase soil water evaporation as indicated by the positive relationship between soil water content and lc‐excess of surface soil water. In contrast, the evergreen coffee shrubs had a higher leaf area index during the “drier” season, leading to enhanced rainfall interception and a negative relationship between lc‐excess and soil water content. During the wet season, there was no clear relationship between soil water content and between lc‐excess of surface soil water. Greater surface soil water under coffee during the dry season may, in part, explain greater preferential flow under coffee compared with under trees in conditions of low rainfall intensities. However, with increasing rainfall intensities during the wet season, there was no obvious difference in preferential flow between the two canopy covers. Results from this study indicate that our new conceptual model can be used to help disentangling the relative influence of canopy cover on local soil water isotopic composition and dynamics, yet also stresses the need for additional measurements to better resolve the underlying processes by which canopy structure influences local water dynamics.     

Spectral-type simulation of spatially correlated fracture set properties

Fracture set properties such as orientation, spacing, trace length, and waviness tend to be spatially correlated. These properties can be efficiently simulated by spectral analysis procedures that take advantage of the computational speed of the fast Fourier transform. The covariance function of each property to be simulated is obtained from the variogram function estimated from mapped fracture set data and is typically referenced to the mean vector of the set. Simulation procedures for normally and exponentially distributed data involve generating uncorrelated Fourier coefficients that are assigned proper variance according to the spectral density, which is the Fourier transform of the covariance function. These coefficients are then reverse Fourier transformed to produce simulated set properties that have the desired variance and variogram function.     

On a class of variational equations transformable to the Gauss hypergeometric equation

A new class of linear ordinary differential equations with periodic coefficients is found which can be transformed to the Gauss hypergeometric equation, and therefore the monodromy matrices are computable explicitly. These equations appear as the variational equations around a straight-line solution in Hamiltonian systems of the form H = T(p) + V(q), where T(p) and V(q) are homogeneous functions of p and q, respectively.     

单元刚度矩阵的几何量分析

本文利用向量和矩阵的运算推导了单元刚度矩阵的几何量表示,分析了受三角形单元和矩形单元几何形状影响的单元刚度矩阵的性质,并将有限元网格与总体刚度矩阵相联系,直接说明了合理剖分区域Ω的基本原则。     

Wavelet analysis and applications to some dynamical systems

For about five years the wavelet transform technique has given interesting results in a great number of different fields such as mathematics, quantum mechanics, signal analysis and image processing, cluster analysis... The wavelet transform appears as a new time-frequency method which is particularly well-suited to detect and to localize discontinuities and scaling behaviours in signals. The main properties of the wavelet transform and its improvements over classical analyzing methods are summarized. Some results among the first applications to the dynamical systems are presented: solution of partial differential equations, fractal and turbulence characterization, and asteroid family determination from cluster analysis.     

Estimating macropore and matrix flow using the hydrograph separation procedure in an experimental forest plot

This study was conducted to estimate macropore space, macropore flow and matrix flow in an experimental forest plot in the Ouachita Mountains of Arkansas. Lateral soil water fluxes and soil capillary potentials were observed in the isolated plot during applied rainfall experiments. Rainfalls were applied 17 times during the period 17 July to 10 October 1991. The subsurface hydrograph separation technique was used to estimate macropore space, macropore flux and matrix flux. The boundary between macropore and matrix flow was statistically determined by covariance analysis. The maximum estimated lateral macropore space was approximately 0.006 (cm³ cm^?3). The maximum estimated lateral macropore and matrix flow were 0.042 and 0.00066 cm s^?1, respectively. This report also emphasizes the need for further research on the hydrograph separation procedure for estimating macropores and macropore flow.