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.     

Radiation parameters of an especially great radio microwave burst

By using the data of an especially great solar radio microwave burst occurring at 1700 UT on 30 October, 1992 in a very active active region NOAA/USAF Region 7321, on the basis of the theoretical emission mechanism of nonthermal gyro-synchrotron radiation, we have derived the radiation parameters such as the energy distribution of energetic electrons in the source region of the burst, the emission coefficient, absorption coefficient and source function of the gyro-synchrotron radiation and made a simple analysis.     

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.     

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.     

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.     

Modeling the seismic response of basin-edge structure by a sequential waveform method

A sequential waveform method is developed to simulate the seismic response of basin-edge structure excited by a plane incident P-wave. The full procedure involves: (a) a previous parameterization of the investigated model using the seismic wave velocities and depths of the sedimentary stratifications; (b) an input motion determined from the records at stations installed on hard rock; (c) forward computation of the P-SV elastic wave field by means of a two-dimensional finite difference (FD) method; (d) the optimization of the model vector using simulated annealing technique and comparing the simulated seismic response of the tested structure with the observed wave field; (e) the correction of the initial model by trial-and-error by testing the differences between synthetics and observed data, and (f) the final solution obtained by iteration using the conjugate gradient algorithm. The search of an optimal basin-edge model has been parallel processed by varying the shapes and velocities of strata on the basis of the fitting of relative timing, amplitude and phase between the output and the observed data. The input motion and sensitivity have been checked and the validity of the method has been demonstrated by numeric analysis. Using the teleseismic records generated by 7 earthquakes recorded at 26 broadband seismic stations, we have studied the seismic velocity structure of the southern edge of the Jiyang depression located in the Bohai Bay basin, northern China. Two cross sections show an agreement between the velocity results and the geological sections available in the region. In addition, we obtain evidence of three hidden faults under the sections and features that suggest major extensions at the Paleogene.