2-D crustal Poisson′s ratio from seismic travel time inversion in Changbaishan Tianchi volcanic region

Based on the inversion method of 2D velocity structure and interface, the crustal velocity structures of P-wave and S-wave along the profile L1 are determined simultaneously with deep seismic sounding data in Changbaishan Tianchi volcanic region, and then its Poisson's ratio is obtained. Calculated results show that this technique overcomes some defects of traditional forward calculation method, and it is also very effective to determine Poisson's ratio distribution of deep seismic sounding profile, especially useful for study on volcanic magma and crustal fault zone. Study result indicates that there is an abnormally high Poisson's ratio body that is about 30 km wide and 12 km high in the low velocity region under Tianchi crater. Its value of Poisson's ratio is 8% higher than that of surrounding medium and it should be the magma chamber formed from melted rock with high temperature. There is a high Poisson's ratio zone ranging from magma chamber to the top of crust, which may be the uprise passage of hot substance. The lower part with high Poisson's ratio, which stretches downward to Moho, is possibly the extrusion way of hot substance from the uppermost mantle. The conclusions above are consistent with the study results of both tomographic determination of 3D crustal structure and magnetotelluric survey in this region.

Geoacoustic Inversion Using Time Reversal of Ocean Noise

We present a passive geoacoustic inversion method using two hydrophones, which combines noise interferometry and time reversal mirror （TRM） techniques. Numerical simulations are firstly performed, in which strong fo- cusing occurs in the vicinity of one hydrophone when Green＇s function （GF） is back-propagated from the other hydrophone, with the position and strength of the focus being sensitive to sound speed and density in the bottom. We next extract the GF from the noise cross-correlation function measured by two hydrophones with 8025-m distance in the Shallow Water ＇06 experiment. After realizing the TRM process, sound speed and density in the bottom are inverted by optimizing focusing of the back-propagated GF. The passive inversion method is inherently environmentally friendly and low-cost.

Major methods of seismic anisotropy

Seismic anisotropy reveals that seismic wave velocity, amplitude, and other physical properties show variations in different directions, which can be divided into lattice-preferred orientation(LPO) and shape-preferred orientation(SPO) according to its physical mechanisms. The main methods for studying seismic anisotropy include shearwave splitting analysis, P-wave travel time inversion and surface-wave tomography, etc. There are some differences and correlations among these methods. Seismic anisotropy is an important way to reveal the dynamic processes of crust-mantle evolution, and it is significant for monitoring crustal stress changes and improve seismic exploration studies. With the help of long-term observation, the application of machine learning techniques and combining inversion based on multiple phases would become potential developments in seismic anisotropy studies. This may improve the understanding of complex seismic anisotropic models, such as multiple layers anisotropy with an oblique axis of symmetry.

Suppression to the cross-channel interference based on the short time Fourier transform

A new cross-channel interference suppression method is proposed to decrease the cross-channel interference in beat signals based on the short time Fourier transform （STY3＂） and the inverse short time Fourier transform （ISTFT） when the dual-orthogonal polarimetric frequency-modulated continu- ous wave （FMCW） radar adopts the opposite-slope linear frequency modulation signal pair in the simultaneous measurement mode. The STFT is applied only on the signals in the cross-interference intervals in the four polarimetric channels to decrease the computation complexity. A mask matrix for suppressing the interference is constructed using the constant false alarm ratio （CFAR） detection on the spectrograms by the STFY. The simulative results show that the cross-channel interference is effi- ciently suppressed by the proposed method. The comparison between the proposed method and the rejection method verifies the improved performance of the proposed method.

Tunnel seismic tomography method for geological prediction and its application

Typical existing methods of tunnel geological prediction include negative apparent velocity, horizontal seismic profile, and the Tunnel Seismic Prediction （TSP） method as this technology is under development at home and abroad. Considering simpler observational methods and data processing, it is hard to accurately determine the seismic velocity of the wall rock in the front of the tunnel face. Therefore, applying these defective methods may result in inaccurate geological inferences which will not provide sufficient evidence for classifying the wall rock characteristics. This paper proposes the Tunnel Seismic Tomography （TST） method using a spatial observation arrangement and migration and travel time inversion image processing to solve the problem of analyzing the velocity structure of wall rock in the front of the tunnel face and realize accurate imaging of the geological framework of the tunnel wall rock. This method is very appropriate for geological prediction under complex geological conditions.

The seismic CT method in measuring rock bodies

The seismic computed tomography （CT） method is derived from the basic principles of X-ray section scanning first applied in medical science. The method records P-wave arrivals between shots and receivers in separate boreholes. Using the velocity information from 2D seismic P-wave arrival inversion, we can determine the distribution of velocity in rock and soil bodies. This paper introduces a practical case of using the seismic CT method for detecting the structure of the rocks for tunnel engineering and to utilize SIRT algorithms for doing first arrival time iterative inversion. Compared with other exploration methods, it is more efficient and accurate.

Characteristic of crustal structure beneath the rifts in southern Tibetan plateau

The method of 2-D travel time inversion, which can be applied to determining 2-D velocity structure and interfaces simultaneously, is used in this paper to reprocess the data of Paiku Co-Pumoyingcuo seismic profile across the Nyima-Tingri rift and Shenzha-Dinggye rift. P-wave velocity structure and interfaces beneath the profile are obtained. The interfaces in the crust near Tingri and Dinggye which are located on rifts have a tendency to uplift, and velocities of middle and lower crusts are high. Low velocity layer in upper crust has an offset. Compared with the distribution of the earthquakes in this region, it is speculated that normal faults near Tingri and Dinggye extend to the upper mantle. Apparently it is affected by deep material： the uplift of mantle causes partial melting in the crust, thus the thickness of crust in this area becomes thin, and tension failures occur in this region easily. On the basis of the characteristics of the earthquakes＇ distribution and the structures of the crustal velocity and interfaces, materials from the mantle still uplifts and the failures are still active.

Study on fine structure of crust-mantle transi-tion zone in Yanqing-Hailai basin based on CDP and DSS data

The fine structure of crust mantle transition zone in Yanqing Hailai(Yan Huai basin) basin has been analyzed and discussed by using the data of Yanqing Hailai deep seismic reflection profile and Beijing Huailai Fengzhen deep seismic sounding profile obtained respectively in 1992 and in 1993. The primary model is established based on CDP stacking profile. The phases of seismic refraction waves and wide angle reflection waves are analyzed, travel time inversion is carried out and 2 D ray tracing is computed. Synthetic seismograms are completed by using re flectivity method for waveform fitting of phases P M and P M′, they are reflection waves from both the upper and the lower boundaries of the crust mantle transition zone in this basin. The results show that the P M′ reflection waves are stronger at some points and come from the lower boundary of Moho discontinuity. It is confirmed and inter preted that the Moho discontinuity in Yanqing Hailai basin consists of a group of thin layers with alternatively higher and lower velocities.

Delight and Frustration with Number “Seven” in Plane Geometry and the Regular Heptagon

As starting point for patterns with seven-fold symmetry, we investigate the basic possibility to construct the regular heptagon by bicompasses and ruler. To cover the whole plane with elements of sevenfold symmetry is only possible by overlaps and (or) gaps between the building stones. Resecting small parts of overlaps and filling gaps between the heptagons, one may come to simple parqueting with only a few kinds of basic tiles related to sevenfold symmetry. This is appropriate for parqueting with a center of seven-fold symmetry that is illustrated by figures. Choosing from the basic patterns with sevenfold symmetry small parts as elementary stripes or elementary cells, one may form by their discrete translation in one or two different directions periodic bordures or tessellation of the whole plane but the sevenfold point-group symmetry of the whole plane is then lost and there remains only such symmetry in small neighborhoods around one or more centers. From periodic tiling, we make the transition to aperiodic tiling of the plane. This is analogous to Penrose tiling which is mostly demonstrated with basic elements of fivefold symmetry and we show that this is also possible with elements of sevenfold symmetry. The two possible regular star-heptagons and a semi-regular star-heptagon play here a basic role.

To Construct High Level Secure Communication System： CTMI Is Not Enough

Public key cryptographic （PKC） algorithms, such as the RSA, elliptic curve digital signature algorithm （ECDSA） etc., are widely used in the secure communication sys- tems, such as OpenSSL, and a variety of in- formation security systems. If designer do not securely implement them, the secret key will be easily extracted by side-channel attacks （SCAs） or combinational SCA thus mitigat- ing the security of the entire communication system. Previous countermeasures of PKC im- plementations focused on the core part of the algorithms and ignored the modular inversion which is widely used in various PKC schemes. Many researchers believe that instead of straightforward implementation, constant time modular inversion （CTMI） is enough to resist the attack of simple power analysis combined with lattice analysis. However, we find that the CTMI security can be reduced to a hidden t-bit multiplier problem. Based on this feature, we firstly obtain Hamming weight of interme- diate data through side-channel leakage. Then, we propose a heuristic algorithm to solve the problem by revealing the secret （partial and full） base of CTMI. Comparing previous nec-essary input message for masking filtering, our procedure need not any information about the secret base of the inversion. To our knowl- edge, this is the first time for evaluating the practical security of CTM! and experimental results show the fact that CTMI is not enough for high-level secure communication systems.

MATLAB fitting for overcurrent relay and a new digital algorithm study

The overcurrent (OC) protection limit is set usually accorging to a OC protection setting table on digital integrated protection equipment in mine explode isolation high voltage (HV) vacuum switch. For digital integrated protection equipment, OC protection setting table must be converted to be a microcomputer algorithm. This paper first intro-duced a method of the fitting OC protection setting table to be OC relay inverse time characteristics equations using MATLAB least square fitting. On the basis of analyzing these fitting equations, a notion, “integral limit rate” was put forward initially and a OC in-verse time digital algorithm was developed. MATLAB simulation results and a digital signal processor (DSP) based digital integrated protection equipment running test indicate that this algorithm has less calculation amount, less taking up memory, high control accuracy, implements the no-grade setting of OC delay values, suits for all kinds of low-middle mi-crocomputer system implementation.

Time evolution of atomic inversion in a standing wave light field

The interaction between an atomic beam of two-level atoms and a standing wave light field has been studied by the exact solution of a time-dependent quantum system developed recently. When the initial atomic state is choosen to be ground, we find that with the limit of zero detuning the atoms will oscillate between the upper and the lower levels with a decaying amplitude. The most interesting result obtained in this paper is when the initial atomic state is a particular superposition of the two levels, now the system does not oscillate at any time.

The Analysis of Spatial-Temporal Evolution of City Accessibility Based on Highway Network in Henan Province in China

Accessibility is an important tool to evaluate the maturity of a regional traffic network structure which describes the traffic convenience in the traffic network. The paper defines a new accessibility index by using the resident population weighted average value of the sum of inverse of the traveling time distance and time threshold coming from ordinary traffic network, and then uses this accessibility index to analyze the spatial-temporal characteristics of Henan highway network, as well as its evolution patterns from 2005 to 2020. The results show that with the expansion and improvement of Henan highway network, city accessibility level has been significantly improved, spatial convergence is obvious, the cities in the north central are always High-High aggregation area, the cities in the south are always Low-Low aggregation area, gradually forming the characteristics of Northwest high and Southeast low, relative balance between East and West. There is some non-conforming phenomenon in highway mileage growth and improvement of the city accessibility levels, but this situation is being weakened, the highway network layout is gradually rationalized, the spatial distribution of city accessibility and that of population are beginning to converge.

Some Equations for the External Space and Their Implications on the Inverse Cone Concept and the Geometry of the Physical Universe

In previous works, they were proposed a photonic model of the Big Bang [1] and several parameters derived from the Hubble-Lemaitre equation [2]. Since these parameters result higher than the classical ones and, otherwise, the General Theory of Relativity does not apply far away the Physical Universe, in this paper, it will be revised the adequacy of such parameters in the external Space and their influence on the relativistic concept of the cone of time. As well, it will be intended to define the Physical Universe geometry accordingly to a thermo-dynamical analysis of the Big Bang.

Real-time trajectory planning for UCAV air-to-surface attack using inverse dynamics optimization method and receding horizon control

This paper presents a computationally efficient real-time trajectory planning framework for typical unmanned combat aerial vehicle （UCAV） performing autonomous air-to-surface （A/S） attack. It combines the benefits of inverse dynamics optimization method and receding horizon optimal control technique. Firstly, the ground attack trajectory planning problem is mathematically formulated as a receding horizon optimal control problem （RHC-OCP）. In particular, an approximate elliptic launch acceptable region （LAR） model is proposed to model the critical weapon delivery constraints. Secondly, a planning algorithm based on inverse dynamics optimization, which has high computational efficiency and good convergence properties, is developed to solve the RHCOCP in real-time. Thirdly, in order to improve robustness and adaptivity in a dynamic and uncer- tain environment, a two-degree-of-freedom （2-DOF） receding horizon control architecture is introduced and a regular real-time update strategy is proposed as well, and the real-time feedback can be achieved and the not-converged situations can be handled. Finally, numerical simulations demon- strate the efficiency of this framework, and the results also show that the presented technique is well suited for real-time implementation in dynamic and uncertain environment.