电离层参量M(3000)F2和hmF2的经验正交函数方法建模研究

The ionospheric parameter M(3000)F2 (the so-called transmission factor or the propagation factor) is important not only in practical applications such as frequency planning for radio-communication but also in ionospheric modeling. This parameter is strongly anti-correlated with the ionospheric F2-layer peak height hmF2，a parameter often used as a key anchor point in some widely used empirical models of the ionospheric electron density profile (e.g., in IRI and NeQuick models). Since hmF2 is not easy to obtain from measurements and M(3000)F2 can be routinely scaled from ionograms recorded by ionosonde/digisonde stations distributed globally and its data has been accumulated for a long history, usually the value of hmF2 is calculated from M(3000)F2 using the empirical formula connecting them. In practice, CCIR M(3000)F2 model is widely used to obtain M(3000)F2 value. However, recently some authors found that the CCIR M(3000)F2 model has remarkable discrepancies with the measured M(3000)F2, especially in low-latitude and equatorial regions. For this reason, the International Reference Ionosphere (IRI) research community proposes to improve or update the currently used CCIR M(3000)F2 model. Any efforts toward the improvement and updating of the current M(3000)F2 model or newly development of a global hmF2 model are encouraged. In this dissertation, an effort is made to construct the empirical models of M(3000)F2 and hmF2 based on the empirical orthogonal function (EOF) analysis combined with regression analysis method. The main results are as follows: 1. A single station model is constructed using monthly median hourly values of M(3000)F2 data observed at Wuhan Ionospheric Observatory during the years of 1957–1991 and compared with the IRI model. The result shows that EOF method is possible to use only a few orders of EOF components to represent most of the variance of the original data set. It is a powerful method for ionospheric modeling. 2. Using the values of M(3000)F2 observed by ionosondes distributed globally, data at grids uniformly distributed globally were obtained by using the Kriging interpolation method. Then the gridded data were decomposed into EOF components using two different coordinates: (1) geographical longitude and latitude; (2) modified dip (Modip) and local time. Based on the EOF decompositions of the gridded data under these two coordinates systems, two types of the global M(3000)F2 model are constructed. Statistical analysis showed that the two types of the constructed M(3000)F2 model have better agreement with the observational M(3000)F2 than the M(3000)F2 model currently used by IRI. The constructed models can represent the global variations of M(3000)F2 better. 3. The hmF2 data used to construct the hmF2 model were converted from the observed M(3000)F2 based on the empirical formula connecting them. We also constructed two types of the global hmF2 model using the similar method of modeling M(3000)F2. Statistical analysis showed that the prediction of our models is more accurate than the model of IRI. This demonstrated that using EOF analysis method to construct global model of hmF2 directly is feasible. The results in this thesis indicate that the modeling technique based on EOF expansion combined with regression analysis is very promising when used to construct the global models of M(3000)F2 and hmF2. It is worthwhile to investigate further and has the potential to be used to the global modeling of other ionospheric parameters.

西南极南设得兰群岛中新生代火山岩基底沉积学及其大地构造制约

Livingston Island, the second island of South Shetland Island, constains Mesozoic-Cenozoic basement, Mesozoic-Cenozoic volcanic sequences, plutonic intrusions and post-subduction volcanic rocks, which document the history and evolution of an important part of the South Shetland Islands magmatic arc. The sedimentary sequence is named the Miers Bluff Formation (MBF) and is interpreted as turbidite since the first geological study on South Shetland Islands, and is interpreted as turbidite. It base and top are not exposed, but a thickness of more than 3000m has been suggested and seems plausible. The turbidite is overlain by Mid - Cretaceous volcanic rocks and intruded by Eocene tonalites. The age of the Miers Bluff Formation is poorly constrained Late Carboniferous -Early Triassic. Sedimentary Environment, tectonic setting and forming age of sedimentary rocks of the Miers Bluff Formation were discussed by means of the methods of sedimentology, petrography and geochemistry, combinig with the study of trace fossils and microfossil plants. The following conclusions are obstained. A sedimentary geological section of Johnsons Dock is made by outside measuring and watching, and then according the section, the geological map near the Spanish Antarctic station was mapped. Four pebbly mudstone layers are first distinguished, which thickness is about 10m. The pebbly mudstone is the typical rock of debris flow, and the depostional environment of pebbly mudstone may be the channel of mid fan of submarine fan. The sedimentsry structural characteristics and size analysis of sandstones show the typical sedimentary feature of turbidity flow and the Miers Bluff Formation is a deep-water turbidite (include some gravity-flow sediments). The materials of palaeocurrents suggest the continental slope dip to southeast, and indicate the provenance of turbidity sediment in the northwest area. By facies analysis, six main facies which include seven subfacies were recognized, which are formed in mid-fan and lower-fan of submarine, meanwhile, the sedimentary features of each facies and subfacies are summarized. The study of clastic composition, major elements, trace elements and rare earth elements indicates the forming setting of the Miers Bluff Formaton is active continental margin and continental island arc and the provenance is dissected magmatic arc which main composition is felsic gneiss. Many trace fossils of the whole succession were found in the turbidites of the Miers Bluff Formation. All these trace fossils are deep sea ichnofossils. There are fifteen ichnogenus, sixteen ichnospecies. Moreover, a new trace fossil was found and a new ichnogenus and new ichnospecies was proposed - Paleaichnus antarctics ichnogen, et ichnosp, nov.. Except the new ichnogenus and ichnospecies, others had been found in deep-sea flysch turbidites. Some are in mudstone and are preserved in the cast convex of overlying sandstone sole, they formed before turbidity flows occurred and belong to the high-different Graphoglyptida of fiysch mudstone. Others as Fucusopsis and Neonereites are preserved in sandstones and stand for trace assemblages after turbidity sedimentation. These trace fossils are typical members of abyssal "Nereites" ichnofacies, and provide for the depositional environment of the Miers Bluff Formation. Fairly diverse microfossil plants have been recovered from the Miers Bluff Formation, Livingston Island, including spores, pollen, acritarchs, wood fragments and cuticles. Containing a total of about 45 species (forms) of miospores, the palynofiora is quantitatively characterized by the dominance of non-striate bisaccate pollen, but spores of pteridophytes and pollen of gymnosperms are proportionate in diversity. It is somewhat comparable to the subzone C+D of the Alisporites zone of Antarctica, and the upper Craterisporites rotundus zone and the lower Polycingulatisporites crenulatus zone of Australia, suggesting a Late Triassic (possibly Norian-Rhaetian) age, as also evidenced by the sporadic occurrence of Aratrisporites and probable Classopollis as well as the complete absence of bisaccate Striatiti. The parent vegetation and paleoclimate are preliminarily deduced. At last, the paper prooses the provenance of sedimentary rocks of the Miers Bluff Formation locates in the east part to the southern Chile(or Southern South American). In the Triassic period, contrasting with New Zealand, Australia and South American of the Pacific margin of Gondwanaland, the Miers Bluff Formation is deposited in the fore-arc basin or back-arc basin of magmatic arc.

地质体科学计算可视化研究及其应用

Describing visually space-time properties of geological phenomena consists of one of the most important parts in geology research. Such visual images are of usually helpful for analyzing geological phenomena and for discovering the regulations behind geological phenomena. This report studies mainly three application problems of scientific visualization in geology: (Dvisualizing geological body A new geometric modeling technique with trimmed surface patches has been eveloped to visualize geological body. Constructional surfaces are represented as trimmed surfaces and a constructional solid is represented by the upper and lower surface composed of trimmed surface patches from constructional surfaces. The technique can completely and definitely represent the structure of geological body. It has been applied in visualization for the coal deposit in Huolinhe, the aquifer thermal energy storage in Tianjin and the structure of meteorite impact in Cangshan et al. (2)visualizing geological space field Efficient visualization methods have been discussed. Marching-Cube algorithm used has been improved and is used to extract iso~surface from 3D data set, iso-line from 2D data set and iso-point from ID data set. The improved method has been used to visualize distribution and evolution of the abnormal pressures in Zhungaer Basin. (3)visualizing porous space a novel way was proposed to define distance from any point to a convex set. Thus a convex set skeleton-based implicit surface modeling technique is developed and used to construct a simplified porous space model. A Buoyancy Percolation numerical simulation platform has been developed to simulate the process of migration of oil in the porous media saturated with water.

开发早期油藏动态地质模型与油藏工程优化研究——以埕岛油田馆上段油藏为例

In order to developing reservoir of Upper of Ng at high-speed and high-efficient in Chengdao oilfield which is located in the bally shallow sea, the paper builds up a series of theory and means predicting and descripting reservoir in earlier period of oilfield development. There are some conclusions as follows. 1. It is the first time to form a series of technique of fine geological modeling of the channel-sandy reservoir by means of mainly seismic methods. These technique include the logging restriction seismic inversion, the whole three dimension seismic interpretation, seismic properties analysis and so on which are used to the 3-dimension distributing prediction of sandy body, structure and properties of the channel reservoir by a lot of the seismic information and a small quantity of the drilling and the logging information in the earlier stage of the oil-field development. It is the first time that these methods applied to production and the high-speed development of the shallow sea oilfield. The prediction sandy body was modified by the data of new drilling, the new reservoir prediction thinking of traced inversion is built. The applied effect of the technique was very well, according to approximately 200 wells belonging to 30 well groups in Chengdao oilfield, the drilling succeeded rate of the predicting sandy body reached 100%, the error total thickness only was 8%. 2. The author advanced the thinking and methods of the forecasting residual-oil prediction at the earlier stage of production. Based on well data and seismic data, correlation of sediment units was correlated by cycle-correlation and classification control methods, and the normalization and finely interpretation of the well logging and sedimentation micro-facies were acquired. On the region of poor well, using the logging restriction inversion technique and regarding finished drilling production well as the new restriction condition, the sand body distributing and its property were predicted again and derived 3-dimension pool geologic model including structure, reservoir, fluid, reservoir engineering parameter and producing dynamic etc. According to the reservoir geologic model, the reservoir engineering design was optimized, the tracking simulation of the reservoir numerical simulation was done by means of the dynamic data (pressure, yield and water content) of development well, the production rule and oil-water distributing rule was traced, the distributing of the remaining oil was predicted and controlled. The dynamic reservoir modeling method in metaphase of development was taken out. Based on the new drilling data, the static reservoir geologic model was momentarily modified, the research of the flow units was brought up including identifying flow units, evaluating flow units capability and establishing the fine flow units model; according to the dynamic data of production and well testing data, the dynamic tracing reservoir description was realized through the constant modification of the reservoir geologic model restricted these dynamic data by the theory of well testing and the reservoir numerical simulation. It was built the dynamic tracing reservoir model, which was used to track survey of the remaining oil on earlier period. The reservoir engineering tracking analysis technique on shallow sea oilfield was founded. After renewing the structure history since tertiary in Chengdao area by the balance section technique and estimating the activity character of the Chengbei fault by the sealing fault analysis technique, the meandering stream sediment pattern of the Upper of Ng was founded in which the meandering border was the uppermost reservoir unit. Based on the specialty of the lower rock component maturity and the structure maturity, the author founded 3 kinds of pore structure pattern in the Guanshang member of Chengdao oil-field in which the storing space mainly was primary (genetic) inter-granular pore, little was secondary solution pore and the inter-crystal pore tiny pore, and the type of throat mainly distributed as the slice shape and the contract neck shape. The positive rhythmic was briefly type included the simple positive rhythm, the complex positive rhythm and the compound rhythm. Interbed mainly is mudstone widely, the physical properties and the calcite interbed distribute localized. 5. The author synthetically analyzed the influence action of the micro-heterogeneity, the macro-heterogeneity and the structure heterogeneity to the oilfield water flood development. The efficiency of water flood is well in tiny structure of convex type or even type at top and bottom in which the water breakthrough of oil well is soon at the high part of structure when inject at the low part of structure, and the efficiency of water flood is poor in tiny structure of concave type at top and bottom. The remaining oil was controlled by sedimentary facies; the water flooding efficiency is well in the border or channel bar and is bad in the floodplain or the levee. The separation and inter layer have a little influence to the non-obvious positive rhythm reservoir, in which the remaining oil commonly locate within the 1-3 meter of the lower part of the separation and inter layer with lower water flooding efficiency.

Graphitic nanofilaments as novel support of Ru-Ba catalysts for ammonia synthesis

A novel graphitic-nanofilament-(GNF-) supported Ru-Ba catalyst is prepared and used in ammonia synthesis reaction. The Ru-Ba/GNFs catalyst shows remarkably high activity and stability for ammonia synthesis, which can be attributed to high purity and graphitization of GNFs with unique structure. TEM micrographs of the Ru-Ba/GNFs catalysts show that Ru metal particles uniformly disperse on the outer wall of GNFs, and the particles become bigger than that before ammonia synthesis reaction after 50 h of operation at 500degreesC and 7.0 MPa, probably due to the Ru crystals covered by promoter and support materials and/or sintering of Ru crystals. (C) 2002 Elsevier Science (USA).

Investigation of platinum utilization and morphology in catalyst layer of polymer electrolyte fuel cells

Platinum utilization in the gas-diffusion catalyst layer and thin-film catalyst layer is investigated. The morphology of PTFE and Nafion in a simulated catalyst layer is examined by scanning electronmicroscopy (SEM) and transmission electron microscopy (TEM). The results show that the platinum utilization of the thin-film catalyst layer containing only Pt/C and Nafion is 45.4%. The low utilization is attributed to the fact that the electron conduction of many catalyst particles is impaired by some thick Nafion layers or clumps. For the gas-diffusion (E-TEK) electrode, the platinum utilization is mainly affected by the proton conduction provided by Nafion. The blocking effect of PTFE on the active sites is not serious. When the electrode is sufficiently impregnated with Nafion by an immersion method, the platinum utilization can reach 77.8%. Transmission electron micrographs reveal that although some thick Nafion layers and clumps are observed in the Pt/C + Nafion layer, the distribution of Nafion in the catalyst layer is basically uniform. The melted PTFE disperses in the catalyst layer very uniformly. No large PTFE clumps or wide net-like structure is observed. The reactant gas may have to diffuse evenly in the catalyst layer. (C) 1999 Elsevier Science S.A. All rights reserved.

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.

Curved Inertia Frames: Visual Attention and Perceptual Organization Using Convexity and Symmetry

In this paper we present an approach to perceptual organization and attention based on Curved Inertia Frames (C.I.F.), a novel definition of "curved axis of inertia'' tolerant to noisy and spurious data. The definition is useful because it can find frames that correspond to large, smooth, convex, symmetric and central parts. It is novel because it is global and can detect curved axes. We discuss briefly the relation to human perception, the recognition of non-rigid objects, shape description, and extensions to finding "features", inside/outside relations, and long- smooth ridges in arbitrary surfaces.

Perceptual Organization, Figure-Ground, Attention and Saliency

Notions of figure-ground, inside-outside are difficult to define in a computational sense, yet seem intuitively meaningful. We propose that "figure" is an attention-directed region of visual information processing, and has a non-discrete boundary. Associated with "figure" is a coordinate frame and a "frame curve" which helps initiate the shape recognition process by selecting and grouping convex image chunks for later matching- to-model. We show that human perception is biased to see chunks outside the frame as more salient than those inside. Specific tasks, however, can reverse this bias. Near/far, top/bottom and expansion/contraction also behave similarly.

Extracting Salient Curves from Images: An Analysis of the Saliency Network

The Saliency Network proposed by Shashua and Ullman is a well-known approach to the problem of extracting salient curves from images while performing gap completion. This paper analyzes the Saliency Network. The Saliency Network is attractive for several reasons. First, the network generally prefers long and smooth curves over short or wiggly ones. While computing saliencies, the network also fills in gaps with smooth completions and tolerates noise. Finally, the network is locally connected, and its size is proportional to the size of the image. Nevertheless, our analysis reveals certain weaknesses with the method. In particular, we show cases in which the most salient element does not lie on the perceptually most salient curve. Furthermore, in some cases the saliency measure changes its preferences when curves are scaled uniformly. Also, we show that for certain fragmented curves the measure prefers large gaps over a few small gaps of the same total size. In addition, we analyze the time complexity required by the method. We show that the number of steps required for convergence in serial implementations is quadratic in the size of the network, and in parallel implementations is linear in the size of the network. We discuss problems due to coarse sampling of the range of possible orientations. We show that with proper sampling the complexity of the network becomes cubic in the size of the network. Finally, we consider the possibility of using the Saliency Network for grouping. We show that the Saliency Network recovers the most salient curve efficiently, but it has problems with identifying any salient curve other than the most salient one.

Near-Optimal Distributed Failure Circumscription

Small failures should only disrupt a small part of a network. One way to do this is by marking the surrounding area as untrustworthy --- circumscribing the failure. This can be done with a distributed algorithm using hierarchical clustering and neighbor relations, and the resulting circumscription is near-optimal for convex failures.

Mid-Level Vision and Recognition of Non-Rigid Objects

We address mid-level vision for the recognition of non-rigid objects. We align model and image using frame curves - which are object or "figure/ground" skeletons. Frame curves are computed, without discontinuities, using Curved Inertia Frames, a provably global scheme implemented on the Connection Machine, based on: non-cartisean networks; a definition of curved axis of inertia; and a ridge detector. I present evidence against frame alignment in human perception. This suggests: frame curves have a role in figure/ground segregation and in fuzzy boundaries; their outside/near/top/ incoming regions are more salient; and that perception begins by setting a reference frame (prior to early vision), and proceeds by processing convex structures.

Recognizing 3-D Objects Using 2-D Images

We discuss a strategy for visual recognition by forming groups of salient image features, and then using these groups to index into a data base to find all of the matching groups of model features. We discuss the most space efficient possible method of representing 3-D models for indexing from 2-D data, and show how to account for sensing error when indexing. We also present a convex grouping method that is robust and efficient, both theoretically and in practice. Finally, we combine these modules into a complete recognition system, and test its performance on many real images.

A Modern Differential Geometric Approach to Shape from Shading

How the visual system extracts shape information from a single grey-level image can be approached by examining how the information about shape is contained in the image. This technical report considers the characteristic equations derived by Horn as a dynamical system. Certain image critical points generate dynamical system critical points. The stable and unstable manifolds of these critical points correspond to convex and concave solution surfaces, giving more general existence and uniqueness results. A new kind of highly parallel, robust shape from shading algorithm is suggested on neighborhoods of these critical points. The information at bounding contours in the image is also analyzed.

Image Chunking: Defining Spatial Building Blocks for Scene Analysis

Rapid judgments about the properties and spatial relations of objects are the crux of visually guided interaction with the world. Vision begins, however, with essentially pointwise representations of the scene, such as arrays of pixels or small edge fragments. For adequate time-performance in recognition, manipulation, navigation, and reasoning, the processes that extract meaningful entities from the pointwise representations must exploit parallelism. This report develops a framework for the fast extraction of scene entities, based on a simple, local model of parallel computation.sAn image chunk is a subset of an image that can act as a unit in the course of spatial analysis. A parallel preprocessing stage constructs a variety of simple chunks uniformly over the visual array. On the basis of these chunks, subsequent serial processes locate relevant scene components and assemble detailed descriptions of them rapidly. This thesis defines image chunks that facilitate the most potentially time-consuming operations of spatial analysis---boundary tracing, area coloring, and the selection of locations at which to apply detailed analysis. Fast parallel processes for computing these chunks from images, and chunk-based formulations of indexing, tracing, and coloring, are presented. These processes have been simulated and evaluated on the lisp machine and the connection machine.