Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.     

An examination of directional discontinuities and magnetic polarity changes around interplanetary sector boundaries usingE > 2 keV electrons

Past studies of interplanetary magnetic sector boundaries have been based on the assumption that one can determine the field polarities by comparing the field directions with those of the nominal Parker spiral angles. Previous investigators have found evidence for decreases of B, the magnitude of the magnetic fieldB, and increases of , the angle betweenB and the ecliptic plane, at sector boundaries. Others have argued that the characteristic thickness of sector boundaries exceeds that of tangential discontinuities, making sector boundaries a separate class of structures.We use a simple technique for inferring the polarities of interplanetary magnetic fields based on the assumption thatE > 2 keV electrons are always flowing along the magnetic field away from the Sun. Electron data from the UC Berkeley experiment on the ISEE-3 spacecraft are used to examine periods around several apparent sector boundaries in 1978 and 1979. We compare properties of (a) boundaries with field polarity changes and (b) large-angle ( > 60°) directional discontinuities with no field polarity changes. We find no significant differences between the sector boundaries and the directional discontinuities in terms of associated decreases in B or of values of . These results suggest no significant difference between sector boundaries and directional discontinuities other than the change in field polarities. Within limited statistics we find that about half the polarity changes would not have been identified using a requirement that > 90° and that half of the > 120° discontinuities would have been misidentified as polarity changes.     

Doppler velocity measurements made with a scanning photoelectric magnetograph

We discuss the problems connected with the measurements and evaluation of line-of-sight velocities, obtained with a scanning photoelectric magnetograph using a line-shifter with enhanced sensitivity. We bring arguments for the validity of the results of our photoelectric Doppler velocity recordings. We have found a network of cellularly shaped patterns in the distribution of photo-electrically measured line-of-sight motions, upflowing in the magnetically quiet (blue-shifted) and downflowing in magnetically active (red-shifted) areas of the photosphere, if the mean velocity level is estimated for a sufficiently large measured area. The features of both directions are mutually complementary. We demonstrate the effect of the shift of the reference zero velocity level on the topology of the line-of-sight velocity maps, and the dependence of this level on the size of the area from which it is estimated.     

Solar flare activity and the structure of the coronal neutral line

In this study we analyse the positions of major flares from 1978 and 1979, with respect to the magnetic structure of the solar corona, as described by a potential field model. We find that major flares exhibit no strong association with the neutral line at the chromospheric level. However, when we calculate the neutral line's position at higher and higher altitudes in the corona, we find that major flares show an increasing tendency to be found close to these high-altitude coronal neutral lines. The correlation between flares and higher-altitude coronal neutral lines reaches a maximum at an altitude of 0.35R , and thereafter decreases as the neutral line is moved out to the source surface at an altitude of 1.50R . This indicates that major flares are strongly associated with coronal structure at the 0.35R level ( 250 000 km) - an altitude surprisingly high in the corona. This reinforces the idea that flares are associated with large-scale coronal magnetic fields and also indicates that the region of coronal magnetic topology important to solar flare processes may be larger than previously thought.     

Dynamical evolution of NEAs: Close encounters, secular perturbations and resonances

We discuss the main mechanisms affecting the dynamical evolution of Near-Earth Asteroids (NEAs) by analyzing the results of three numerical integrations over 1 Myr of the NEA (4179) Toutatis. In the first integration the only perturbing planet is the Earth. So the evolution is dominated by close encounters and looks like a random walk in semimajor axis and a correlated random walk in eccentricity, keeping almost constant the perihelion distance and the Tisserand invariant. In the second integration Jupiter and Saturn are present instead of the Earth, and the 3/1 (mean motion) and v ₆ (secular) resonances substantially change the eccentricity but not the semimajor axis. The third, most realistic, integration including all the three planets together shows a complex interplay of effects, with close encounters switching the orbit between different resonant states and no approximate conservation of the Tisserand invariant. This shows that simplified 3-body or 4-body models cannot be used to predict the typical evolution patterns and time scales of NEAs, and in particular that resonances provide some fast-track dynamical routes from low-eccentricity to very eccentric, planet-crossing orbits.On leave from the Department of Mathematics, University of Pisa, Via Buonarroti 2, 56127 Pisa, Italy, thanks to the G. Colombo fellowships of the European Space Agency.     

Shear-induced instability and arch filament eruption: A magnetohydrodynamic (MHD) numerical simulation

We investigate, via a two-dimensional (nonplanar) MHD simulation, a situation wherein a bipolar magnetic field embedded in a stratified solar atmosphere (i.e., arch-filament-like structure) undergoes symmetrical shear motion at the footpoints. It was found that the vertical plasma flow velocities grow exponentially leading to a new type of global MHD-instability that could be characterized as a Dynamic Shearing Instability, with a growth rate of about 8{ovV} _A a, where {ovV} _A is the average Alfvén speed and a ^–1 is the characteristic length scale. The growth rate grows almost linearly until it reaches the same order of magnitude as the Alfvén speed. Then a nonlinear MHD instability occurs beyond this point. This simulation indicates the following physical consequences: the central loops are pinched by opposing Lorentz forces, and the outer closed loops stretch upward with the vertically-rising mass flow. This instability may apply to arch filament eruptions (AFE) and coronal mass ejections (CMEs).To illustrate the nonlinear dynamical shearing instability, a numerical example is given for three different values of the plasma beta that span several orders of magnitude. The numerical results were analyzed using a linearized asymptotic approach in which an analytical approximate solution for velocity growth is presented. Finally, this theoretical model is applied to describe the arch filament eruption as well as CMEs.     

Comparison of structural and least-squares lines for estimating geologic relations

Two different goals in fitting straight lines to data are to estimate a true linear relation (physical law) and to predict values of the dependent variable with the smallest possible error. Regarding the first goal, a Monte Carlo study indicated that the structural-analysis (SA) method of fitting straight lines to data is superior to the ordinary least-squares (OLS) method for estimating true straight-line relations. Number of data points, slope and intercept of the true relation, and variances of the errors associated with the independent (X) and dependent (Y) variables influence the degree of agreement. For example, differences between the two line-fitting methods decrease as error in X becomes small relative to error in Y. Regarding the second goal—predicting the dependent variable—OLS is better than SA. Again, the difference diminishes as X takes on less error relative to Y. With respect to estimation of slope and intercept and prediction of Y, agreement between Monte Carlo results and large-sample theory was very good for sample sizes of 100, and fair to good for sample sizes of 20. The procedures and error measures are illustrated with two geologic examples.     

复杂片岩地区桩基工程实践

复杂片岩地区的桩基工程常出现桩基承载力不足的问题。通过对一典型工程实例的分析，探讨了引起桩基承载力不足的主要原因，及在类似地层桩基施工中应注意的问题。