外差激光雷达测量水体布里渊散射可行性研究

该文分析外差激光雷达测量水体布里渊散射的可行性。借助测量水体布里渊散射的频移 ,可得到水体声速或温度。文中利用两条 12 7I2 吸收线锁定激光输出频率 ,产生稳定的频差约为7.5 GHz发射激光和激光本振。通过对 12 7I2 吸收线稳定度、激光频率稳定度、系统信噪比和反演精度的分析和计算机模拟 ,得出水体深度 2 5 m以内水体布里渊散射频移的测量误差可以控制在2 .5 MHz以内 ,达到水中声速 0 .6 m/ s或温度 0 .1℃的测量精度。     

Use of unit response functions for management of regional multilayered aquifers: application to the North Aquitaine Tertiary system (France)

Although there are many numerical applications used in aquifer management, few methodologies seem to be sufficiently adapted to provide solutions for regional-scale problems. Their applications for multilayered aquifers are also too limited. A new groundwater management tool adapted to such large physical systems has been developed, using the unit response function (URF) approach. This tool is based on the link between a numerical groundwater flow model and a mathematical optimisation tool. Pre- and post processing are based on a geographical information system (GIS). The developed tool has been applied to a real management case—the regional Aquitaine multilayered aquifer (France). A representative hydrodynamic model was built using MODFLOW 2000 code. First tests show that the linearity condition needed for the application of URF methodology is respected despite the high complexity of the system. Two management scenarios were tested. Optimal solutions thus obtained allow for viable exploitation alternatives to be proposed, which integrate complex environmental constraints.

---

Résumé Malgré un nombre important d’applications existant dans le domaine de la gestion des aquifères, peu de méthodologies semblent adaptées pour proposer des solutions à des problèmes réels de dimension régionale. Leurs application pour des systèmes multicouches sont souvent trop limitées. Un outil de gestion adapté à ces systèmes de large taille a été développé, basé sur la théorie des Fonctions de Réponses Unitaires (FRU). Cet outil est basé sur le couplage d’un modèle numérique d’écoulement et d’un code d’optimisation mathématique. Les entrées et sorties du modèle sont réalisées à l’aide d’un Système d’Information Géographique (SIG). L’outil développé est appliqué à un cas de gestion réel—le système aquifère régional multicouche aquitain (France). Un modèle hydrodynamique représentatif a été construit, basé sur le code MODFLOW 2000. Les premiers tests montrent que les conditions de linéarité nécessaires à l’application de la méthode FRU sont respectées, malgré l’importante complexité du système. Deux scénarii de gestion sont testés. Les solutions optimales ainsi obtenues permettent de proposer des solutions d’exploitation alternatives viables, intégrant des contraintes environnementales complexes.

---

Resumen Aunque hay muchas aplicaciones numéricas usadas en la gestión de acuíferos, unas pocas metodologías parecen estar adaptadas suficientemente para dar las soluciones a los problemas de escala regional. También son limitadas sus aplicaciones para los acuíferos multicapa. Se ha desarrollado una nueva herramienta para el manejo del agua subterránea, adaptada a tales sistemas físicos grandes, usando el avance de la Función de Respuesta Unitaria (URF). Esta herramienta se basa en la combinación entre un modelo numérico de flujo de agua subterránea y una herramienta de optimización matemática. El pre-proceso y el post-proceso se basan en un Sistema de Información Geográfico (GIS). La herramienta desarrollada se ha aplicado a un caso real de manejo en el acuífero regional multicapa de Aquitania (Francia). Se construyó un modelo hidrodinámico representativo, que usó el código MODFLOW 2000. Las primeras pruebas muestran que la condición de linearidad, necesaria para la aplicación de la metodología URF, se respeta a pesar de la complejidad alta del sistema. Se probaron dos escenarios de manejo. Las soluciones óptimas así obtenidas permiten la proyección de alternativas viables de explotación, las cuales integran restricciones medioambientales complejas.

     

Weak-contrast approximation of the elastic scattering matrix in anisotropic media

The plane-wave reflection and transmission coefficients at a plane interface between two anisotropic media constitute the elements of the elastic scattering matrix. For a 1-D anisotropic medium the eigenvector decomposition of the system matrix of the transformed elasto-dynamic equations is used to derive a general expression for the scattering matrix. Depending on the normalization of the eigenvectors, the expressions give scattering coefficients for amplitudes or for vertical energy flux.Computing the vertical slownesses and the corresponding polarizations, the eigenvector matrix and its inverse can be found. We give a simple formula for the inverse, regardless of the normalization of the eigenvectors. When the eigenvectors are normalized with respect to amplitudes of displacement (or velocity), the calculation of the scattering matrix for amplitudes is simplified.When the relative changes in all parameters are small, a weak-contrast approximation of the scattering matrix, based on the exactly determined polarization vectors in an average medium, is obtained. The same approximation is also derived directly from the transformed elasto-dynamic equations for a smooth vertically inhomogeneous medium, proving the consistency of the approximation.For monoclinic media, with the mirror symmetry plane parallel to the interface, the approximative scattering matrix is given in terms of analytic expressions for the non-normalized eigenvectors and vertical slownesses. For transversely isotropic media with a vertical axis of symmetry (VTI) and isotropic media, explicit solutions for the weak-contrast approximations of the scattering matrices have been obtained. The scattering matrix for amplitudes for isotropic media is well known. The scattering matrix for vertical energy flux may have applications in AVO analysis and inversion due to the reciprocity of the reflection coefficients for converted waves.Numerical examples for monoclinic and VTI media provide good agreement between the approximative and the exact reflection matrices. It is, however, expected that the approximations cannot be used when the symmetry properties of the two media are very different. This is because the approximation relies on a small relative contrast between the eigenvectors in the two media.Presented at the Workshop Meeting on Seismic Waves in Laterally Inhomogeneous Media, Castle of Trest, Czech Republic, May 22–27, 1995.     

The elastic properties of single-crystal fayalite as determined by dynamical measurement techniques

We present new elasticity measurements on single-crystal fayalite and combine our results with other data from resonance, pulse superposition interferometry, and Brillouin scattering to provide a set of recommended values for the adiabatic elastic moduliC _ij and their temperature variations. We use a resonance method (RPR) with specimens that were previously investigated by pulse superposition experiments. The nineC _ij of fayalite are determined from three new sets of measurements. One set of our newC _ij data is over the range 300–500 K. We believe that the relatively large discrep ncies found in someC _ij are due in large part to specimen inhomogeneities (chemical and microstructural) coupled with differences in the way various techniques sample, rather than only systematic errors associated with experimental procedures or in the preparations of the specimens.Our recommendeaC _ij's (GPa) and (C _ij/T) _p (GPa/K) are: The resulting values for the isotropic bulk and shear moduli,K _s and , and their temperature derivatives are:K _s=134(4) GPa; =50.7(0.3) GPa; (K _s/T) _p =–0.024(0.005) GPa/K; and (/T) _p =–0.013(0.001) GPa/K. An important conclusion is thatK _s increases as the Fe/(Fe+Mg) ratio in olivine is increased.     

Scattering of elastic waves by a fracture zone containing randomly distributed cracks

We theoretically study the scattering ofP, SV andSH waves by a zonal distribution of cracks, which simulates a fault fracture zone. An investigation is conducted how the geometrical properties of the crack distribution and the frictional characteristics of the crack surface are reflected in the attenuation and dispersion of incident waves, as well as in the amplitudes of the transmitted and reflected waves from the zone. If the crack distribution within the fault zone changes temporally during the preparation process of the expected earthquake, it will be important for earthquake prediction to monitor it, utilizing the scattering-induced wave phenomena.We consider the two-dimensional problem. Aligned cracks with the same length are assumed to be randomly distributed in a zone with a finite width, on which elastic waves are assumed to be incident. The distribution of cracks is assumed to be homogeneous and sparse. The crack surface is assumed to be stress-free, or to undergo viscous friction; the latter case simulates fluid-filled cracks. The opening displacement of the crack is assumed to be negligibly small. The idea of the mean wave formalism is employed in the analysis, and Foldy's approximation is assumed.When the crack surface is stress-free, it is commonly observed for every wave mode (P, SV andSH) that the attenuation coefficientQ ^–1 peaks aroundka1, the phase velocity is almost independent ofk in the rangeka<1 and it increases monotonically withk in the rangeka>1, wherek is the intrinsicS wavenumber anda is the half length of the crack. The effect of the friction is to shift the peak ofQ ^–1 and the corner of the phase velocity curve to the low wavenumber range. The high wavenumber asymptote ofQ ^–1 is proportional tok ^–1 independently of model parameters and the wave modes. If the seismological observation thatQ ^–1 ofS waves has a peak at around 0.5 Hz in the earth's crust is combined with our results, the upper limit of crack size within the crust is estimated about 4 km. The information regarding the transmitted and reflected waves, such as the high wavenumber limit of the amplitude of the transmitted wave etc., allows estimation of the strength of the friction.     

Water temperature behaviour in the River Loire since 1976 and 1881

Analysis of monthly mean river temperatures, recorded on an hourly basis in the middle reaches of the Loire since 1976, allows reconstruction by multiple linear regression of the annual, spring and summer water temperatures from equivalent information on air temperatures and river discharge. Since 1881, the average annual and summer temperatures of the Loire have risen by approximately 0.8?°C, this increase accelerating since the late 1980s due to the rise in air temperature and also to lower discharge rates. In addition, the thermal regime in the Orleans to Blois reach is considerably affected by the inflow of groundwater from the Calcaires de Beauce aquifer, as shown by the summer energy balance. To cite this article: F. Moatar, J. Gailhard, C. R. Geoscience 338 (2006).