单（多）井抽灌对浅部地温场的影响研究

The technique of energy extraction using groundwater source heat pumps, as a sustainable way of low-grade thermal energy utilization, has widely been used since mid-1990's. Based on the basic theories of groundwater flow and heat transfer and by employing two analytic models, the relationship of the thermal breakthrough time for a production well with the effect factors involved is analyzed and the impact of heat transfer by means of conduction and convection, under different groundwater velocity conditions, on geo-temperature field is discussed.A mathematical model, coupling the equations for groundwater flow with those for heat transfer, was developed. The impact of energy mining using a single well system of supplying and returning water on geo-temperature field under different hydrogeological conditions, well structures, withdraw-and-reinjection rates, and natural groundwater flow velocities was quantitatively simulated using the finite difference simulator HST3D. Theoretical analyses of the simulated results were also made. The simulated results of the single well system indicate that neither the permeability nor the porosity of a homogeneous aquifer has significant effect on the temperature of the production segment provided that the production and injection capability of each well in the aquifers involved can meet the designed value. If there exists a lower permeable interlayer, compared with the main aquifer, between the production and injection segments, the temperature changes of the production segment will decrease. The thicker the interlayer and the lower the interlayer permeability, the longer the thermal breakthrough time of the production segment and the smaller the temperature changes of the production segment. According to the above modeling, it can also be found that with the increase of the aquifer thickness, the distance between the production and injection screens, and/or the regional groundwater flow velocity, and/or the decrease of the production-and-reinjection rate, the temperature changes of the production segment decline. For an aquifer of a constant thickness, continuously increase the screen lengths of production and injection segments may lead to the decrease of the distance between the production and injection screens, and the temperature changes of the production segment will increase, consequently.According to the simulation results of the single well system, the parameters, that can cause significant influence on heat transfer as well as geo-temperature field, were chosen for doublet system simulation. It is indicated that the temperature changes of the pumping well will decrease as the aquifer thickness, the distance between the well pair and/or the screen lengths of the doublet increase. In the case of a low permeable interlayer embedding in the main aquifer, if the screens of the pumping and the injection wells are installed respectively below and above the interlayer, the temperature changes of the pumping well will be smaller than that without the interlay. The lower the permeability of the interlayer, the smaller the temperature changes. The simulation results also indicate that the lower the pumping-and-reinjection rate, the greater the temperature changes of the pumping well. It can also be found that if the producer and the injector are chosen reasonably, the temperature changes of the pumping well will decline as the regional groundwater flow velocity increases. Compared with the case that the groundwater flow direction is perpendicular to the well pair, if the regional flow is directed from the pumping well to the injection well, the temperature changes of the pumping well is relatively smaller.Based on the above simulation study, a case history was conducted using the data from an operating system in Beijing. By means of the conceptual model and the mathematical model, a 3-D simulation model was developed and the hydrogeological parameters and the thermal properties were calibrated. The calibrated model was used to predict the evolution of the geo-temperature field for the next five years. The simulation results indicate that the calibrated model can represent the hydrogeological conditions and the nature of the aquifers. It can also be found that the temperature fronts in high permeable aquifers move very fast and the radiuses of temperature influence are large. Comparatively, the temperature changes in clay layers are smaller and there is an obvious lag of the temperature changes. According to the current energy mining load, the temperature of the pumping wells will increase by 0.7°C at the end of the next five years. The above case study may provide reliable base for the scientific management of the operating system studied.

高温高压下部分熔融岩石和岩石玻璃弹性波速研究

人类对地球深部结构的认识主要依赖于天然地震的观测资料，高温高压下矿物、岩石和岩浆玻璃的弹性波速测量，是对野外地震波探测资料进行物质反演的重要依据，也为建立地球内部结构模型和地球动力学研究提供重要的实验数据。大量研究证实，高温高压下岩石的部分熔融将形成地震波低速层。然而，前人的研究中，多以岩石的纵波波速（Vp）测量结果来讨论区域地壳结构和低速层的成因，而且很少对实验过程中的中间产物进行观察分析。另一方面，由于玻璃在高温高压下具有特殊的弹性性质，一些科学家推测地球内部岩石的非晶质化也将导致地震波低速层形成。但目前这一推测尚缺乏充分的实验数据支持。为此，作者依托YJ-3000吨大腔体高压实验技术平台，利用脉冲反射法和透射－反射法，完成了： （Ⅰ）三江地区花岗岩和角闪斜长片麻岩在最高压力2.0GPa、最高温度1200℃下的纵波波速（Vp）和最高温度600℃下的横波波速（Vs）研究，并通过岩石物态变化过程中的取样实验，综合探讨岩石中矿物脱水、固－固相变、部分熔融对其弹性波速的影响。获得以下主要结论： ① 花岗岩和角闪斜长片麻岩的Vp和Vs随压力及温度的变化趋势基本一致。室温下岩石的Vp和Vs随压力升高而升高，岩石波速具明显的各向异性，而且其各向异性随压力增大到约0.5GPa后逐渐趋于一恒定值； ② 恒定压力下，岩石的Vp和Vs首先随温度升高近线性缓慢降低，当750℃ 950℃后，石英相变完成，岩石的波速由于熔体含量增加又快速降低； ③ 高温高压下岩石的Vp和Vs研究显示了对三江地区地壳结构一致的约束结果，即该区花岗岩主要分布在上地壳，而角闪斜长片麻岩从上地壳底部到中地壳底部均有分布，这一结果与前人利用Vp研究建立的该区地壳模型基本一致； ④ 在三江地区中上地壳高石英含量的岩石中，石英的α－β相变是地壳地震波低速层形成的主要因素，而随岩石中石英含量的变化，高温高压下岩石的部分熔融及岩石的波速各向异性也可能形成低速层。 （Ⅱ）化学成分从基性到酸性的7种岩石的熔体玻璃在1.0GPa和2.0GPa，最高温度1000℃下的Vp研究和最高温度730℃下的Vs研究。获得以下主要结论： ① 与岩石波速随压力增大而增大不同，室温、0.4－2.0GPa压力下，除两种基性岩石（正长辉石岩和粗面玄武岩）的熔体玻璃外，其它5种中酸性岩石的熔体玻璃的Vp均随压力增大而减小，而这7种玻璃的Vs全部随压力增大而减小。而且，玻璃波速随压力增加而异常降低的幅度随样品中SiO2含量的增高逐渐增大； ② 恒定压力下，随实验温度升高，各种玻璃的弹性波速缓慢降低。当温度高于玻璃转变温度Tg后，玻璃弹性波速的温度系数（∂V/∂T）增大约3倍以上； ③ 研究证实了一种新的地震波低速层成因模式，即上地幔岩石中如果含有超过10vol％的玄武岩玻璃，将会形成地震波低速层；假如下地壳基性岩石中中酸性非晶质体含量超过20vol％，也可能导致地壳低速层的形成。

高温高压下地球深部物质弹性波速测量技术研究和应用

岩石、矿物的弹性波速对地球的地震波速度进行物质成分和晶体结构的解释是十分关键的.该次论文工作的目的是在YJ-3000六面顶压机原有弹性波速测量技术的基础上建立一套P波和S波的弹性波速测量的新方法,可以在模拟地球上地幔的高温高压条件下获取物质的弹性参数.通过对超声学的基本概念的引入和超声测量基本技术的分析,特别是超声探头的采用,结合超声测量中透射法和反射法各自的优点,创造性提出了透射-反射联用的新方法.该方法的使用可以达到较高的走时测量精度,避免了温度和压力梯度对高温高压下弹性波波速测量的影响,同时还结束了多年来该压机上不能进行S波波速测量的历史.采用该方法在高压或同时高温高压下对不同的样品进行了弹性波速的测定,测量结果与前人相同样品的数据吻合得较好,表明了这种新方法的可行性和可靠性.

Direct Computation of 3D Shape Invariants and the Focus of Expansion

Structure from motion often refers to the computation of 3D structure from a matched sequence of images. However, a depth map of a surface is difficult to compute and may not be a good representation for storage and recognition. Given matched images, I will first show that the sign of the normal curvature in a given direction at a given point in the image can be computed from a simple difference of slopes of line-segments in one image. Using this result, local surface patches can be classified as convex, concave, parabolic (cylindrical), hyperbolic (saddle point) or planar. At the same time the translational component of the optical flow is obtained, from which the focus of expansion can be computed.

Colliding and Moving Bose-Einstein Condensates: Studies of superfluidity and optical tweezers for condensate transport

In this thesis, two different sets of experiments are described. The first is an exploration of the microscopic superfluidity of dilute gaseous Bose- Einstein condensates. The second set of experiments were performed using transported condensates in a new BEC apparatus. Superfluidity was probed by moving impurities through a trapped condensate. The impurities were created using an optical Raman transition, which transferred a small fraction of the atoms into an untrapped hyperfine state. A dramatic reduction in the collisions between the moving impurities and the condensate was observed when the velocity of the impurities was close to the speed of sound of the condensate. This reduction was attributed to the superfluid properties of a BEC. In addition, we observed an increase in the collisional density as the number of impurity atoms increased. This enhancement is an indication of bosonic stimulation by the occupied final states. This stimulation was observed both at small and large velocities relative to the speed of sound. A theoretical calculation of the effect of finite temperature indicated that collision rate should be enhanced at small velocities due to thermal excitations. However, in the current experiments we were insensitive to this effect. Finally, the factor of two between the collisional rate between indistinguishable and distinguishable atoms was confirmed. A new BEC apparatus that can transport condensates using optical tweezers was constructed. Condensates containing 10-15 million sodium atoms were produced in 20 s using conventional BEC production techniques. These condensates were then transferred into an optical trap that was translated from the ‘production chamber’ into a separate vacuum chamber: the ‘science chamber’. Typically, we transferred 2-3 million condensed atoms in less than 2 s. This transport technique avoids optical and mechanical constrainsts of conventional condensate experiments and allows for the possibility of novel experiments. In the first experiments using transported BEC, we loaded condensed atoms from the optical tweezers into both macroscopic and miniaturized magnetic traps. Using microfabricated wires on a silicon chip, we observed excitation-less propagation of a BEC in a magnetic waveguide. The condensates fragmented when brought very close to the wire surface indicating that imperfections in the fabrication process might limit future experiments. Finally, we generated a continuous BEC source by periodically replenishing a condensate held in an optical reservoir trap using fresh condensates delivered using optical tweezers. More than a million condensed atoms were always present in the continuous source, raising the possibility of realizing a truly continuous atom lase.

Interaction between coronal mass ejections and the solar wind

Jones, R. A.; Breen, A. R.; Fallows, R. A.; Canals, A.; Bisi, M. M.; Lawrence, G. (2007). Interaction between coronal mass ejections and the solar wind, Journal of Geophysical Research, 112, Issue A8 RAE2008

Newton's aerodynamic problem in media of chaotically moving particles

Plakhov, A.Y.; Torres, D., (2005) 'Newton's aerodynamic problem in media of chaotically moving particles', Sbornik: Mathematics 196(6) pp.885-933 RAE2008

The Impact of Sunspots on the Interpretation of Coronal Observations of the O VI Doublet

Morgan, Huw; Habbal, S. R., 'The impact of sunspots on the interpretation of coronal observations of the OVI doublet', The Astrophysical Journal (2005) 630(2) pp.L189-L192 RAE2008

A Self-Organizing Neural Model for Motor Equivalent Phoneme Production

This paper describes a model of speech production called DIVA that highlights issues of self-organization and motor equivalent production of phonological units. The model uses a circular reaction strategy to learn two mappings between three levels of representation. Data on the plasticity of phonemic perceptual boundaries motivates a learned mapping between phoneme representations and vocal tract variables. A second mapping between vocal tract variables and articulator movements is also learned. To achieve the flexible control made possible by the redundancy of this mapping, desired directions in vocal tract configuration space are mapped into articulator velocity commands. Because each vocal tract direction cell learns to activate several articulator velocities during babbling, the model provides a natural account of the formation of coordinative structures. Model simulations show automatic compensation for unexpected constraints despite no previous experience or learning under these constraints.

Unsupervised Neural Network for the Control of a Mobile Robot

This article introduces an unsupervised neural architecture for the control of a mobile robot. The system allows incremental learning of the plant during robot operation, with robust performance despite unexpected changes of robot parameters such as wheel radius and inter-wheel distance. The model combines Vector associative Map (VAM) learning and associate learning, enabling the robot to reach targets at arbitrary distances without knowledge of the robot kinematics and without trajectory recording, but relating wheel velocities with robot movements.

The Perception of Globally Coherent Motion

How do human observers perceive a coherent pattern of motion from a disparate set of local motion measures? Our research has examined how ambiguous motion signals along straight contours are spatially integrated to obtain a globally coherent perception of motion. Observers viewed displays containing a large number of apertures, with each aperture containing one or more contours whose orientations and velocities could be independently specified. The total pattern of the contour trajectories across the individual apertures was manipulated to produce globally coherent motions, such as rotations, expansions, or translations. For displays containing only straight contours extending to the circumferences of the apertures, observers' reports of global motion direction were biased whenever the sampling of contour orientations was asymmetric relative to the direction of motion. Performance was improved by the presence of identifiable features, such as line ends or crossings, whose trajectories could be tracked over time. The reports of our observers were consistent with a pooling process involving a vector average of measures of the component of velocity normal to contour orientation, rather than with the predictions of the intersection-of-constraints analysis in velocity space.

Optimisation of granola breakfast cereal manufacturing process by wet granulation and pneumatic conveying

This study has considered the optimisation of granola breakfast cereal manufacturing processes by wet granulation and pneumatic conveying. Granola is an aggregated food product used as a breakfast cereal and in cereal bars. Processing of granola involves mixing the dry ingredients (typically oats, nuts, etc.) followed by the addition of a binder which can contain honey, water and/or oil. In this work, the design and operation of two parallel wet granulation processes to produce aggregate granola products were incorporated: a) a high shear mixing granulation process followed by drying/toasting in an oven. b) a continuous fluidised bed followed by drying/toasting in an oven. In high shear granulation the influence of process parameters on key granule aggregate quality attributes such as granule size distribution and textural properties of granola were investigated. The experimental results show that the impeller rotational speed is the single most important process parameter which influences granola physical and textural properties. After that binder addition rate and wet massing time also show significant impacts on granule properties. Increasing the impeller speed and wet massing time increases the median granule size while also presenting a positive correlation with density. The combination of high impeller speed and low binder addition rate resulted in granules with the highest levels of hardness and crispness. In the fluidised bed granulation process the effect of nozzle air pressure and binder spray rate on key aggregate quality attributes were studied. The experimental results show that a decrease in nozzle air pressure leads to larger in mean granule size. The combination of lowest nozzle air pressure and lowest binder spray rate results in granules with the highest levels of hardness and crispness. Overall, the high shear granulation process led to larger, denser, less porous and stronger (less likely to break) aggregates than the fluidised bed process. The study also examined the particle breakage of granola during pneumatic conveying produced by both the high shear granulation and the fluidised bed granulation process. Products were pneumatically conveyed in a purpose built conveying rig designed to mimic product conveying and packaging. Three different conveying rig configurations were employed; a straight pipe, a rig consisting two 45° bends and one with 90° bend. Particle breakage increases with applied pressure drop, and a 90° bend pipe results in more attrition for all conveying velocities relative to other pipe geometry. Additionally for the granules produced in the high shear granulator; those produced at the highest impeller speed, while being the largest also have the lowest levels of proportional breakage while smaller granules produced at the lowest impeller speed have the highest levels of breakage. This effect clearly shows the importance of shear history (during granule production) on breakage during subsequent processing. In terms of the fluidised bed granulation, there was no single operating parameter that was deemed to have a significant effect on breakage during subsequent conveying. Finally, a simple power law breakage model based on process input parameters was developed for both manufacturing processes. It was found suitable for predicting the breakage of granola breakfast cereal at various applied air velocities using a number of pipe configurations, taking into account shear histories.

A miniaturised arrow ballistic measurement system

A novel miniaturised system for measurement of the in-flight characteristics of an arrow is introduced in this paper. The system allows the user to measure in-flight parameters such as the arrow’s speed, kinetic energy and momentum, arrow drag and vibrations of the arrow shaft. The system consists of electronics, namely a three axis accelerometer, shock switch, microcontroller and EEPROM memory embedded in the arrow tip. The system also includes a docking station for download and processing of in-flight ballistic data from the tip to provide the measured values. With this system, a user can evaluate and optimize their archery equipment setup based on measured ballistic values. Recent test results taken at NIST show the accuracy of the launch velocities to be within +/- 0.59%, when compared with NIST’s most accurate ballistic chronograph.

A mathematical theory of stochastic microlensing. II. Random images, shear, and the Kac-Rice formula

Continuing our development of a mathematical theory of stochastic microlensing, we study the random shear and expected number of random lensed images of different types. In particular, we characterize the first three leading terms in the asymptotic expression of the joint probability density function (pdf) of the random shear tensor due to point masses in the limit of an infinite number of stars. Up to this order, the pdf depends on the magnitude of the shear tensor, the optical depth, and the mean number of stars through a combination of radial position and the star's mass. As a consequence, the pdf's of the shear components are seen to converge, in the limit of an infinite number of stars, to shifted Cauchy distributions, which shows that the shear components have heavy tails in that limit. The asymptotic pdf of the shear magnitude in the limit of an infinite number of stars is also presented. All the results on the random microlensing shear are given for a general point in the lens plane. Extending to the general random distributions (not necessarily uniform) of the lenses, we employ the Kac-Rice formula and Morse theory to deduce general formulas for the expected total number of images and the expected number of saddle images. We further generalize these results by considering random sources defined on a countable compact covering of the light source plane. This is done to introduce the notion of global expected number of positive parity images due to a general lensing map. Applying the result to microlensing, we calculate the asymptotic global expected number of minimum images in the limit of an infinite number of stars, where the stars are uniformly distributed. This global expectation is bounded, while the global expected number of images and the global expected number of saddle images diverge as the order of the number of stars. © 2009 American Institute of Physics.

Bunching-induced optical nonlinearity and instability in cold atoms [Invited].

We report a new nonlinear optical process that occurs in a cloud of cold atoms at low-light-levels when the incident optical fields simultaneously polarize, cool, and spatially-organize the atoms. We observe an extremely large effective fifth-order nonlinear susceptibility of χ(⁵) = 7.6 × 10⁻¹⁵ (m/V)⁴, which results in efficient Bragg scattering via six-wave mixing, slow group velocities (∼ c/10⁵), and enhanced atomic coherence times (> 100 μs). In addition, this process is particularly sensitive to the atomic temperatures, and provides a new tool for in-situ monitoring of the atomic momentum distribution in an optical lattice. For sufficiently large light-matter couplings, we observe an optical instability for intensities as low as ∼ 1 mW/cm² in which new, intense beams of light are generated and result in the formation of controllable transverse optical patterns.