Research on High Speed Communication Method of ELF-EM While Drilling Based on Adaptive Combined Filtering Algorithm

In terms of the electromagnetic wave measurement while drilling,the extremely low frequency electromagnetic(ELF-EM)signal below 20Hz is usually chosen as the downhole data carrier.To improve the transmission rate of ELF-EM signal and the signal to noise ratio(SNR)at the receiving end,the DQPSK modulation was proposed as the modulation method for the communication of electromagnetic wave system.Different from the traditional IQ orthogonal modulation and coherent demodulation methods,the proposed phase selection modulation and correlation algorithm demodulation are easier to implement and more practical.With regard to the communication synchronization,a fast algorithm,which based on the normalized cross-relation number,was used for waveform matching,and the maximum point of the correlation coefficient was used as the starting point of communication synchronization.The communication simulation results show that the proposed DQPSK modulation signal based on the adaptive combined filtering algorithm has better terminal error rate and transmission rate than the traditional modulation method.Under the same carrier frequency and code width,the transmission rate of DQPSK modulation is 4 to 5 times and 2 times that of PPM modulation and 2DPSK modulation respectively.The communication modulation and demodulation modes as well as the decoding algorithm with combined adaptive filter proposed in this paper can effectively solve practical engineering problems.

Comparison of the reaction of polar mesosphere winter echoes and polar mesosphere summer echoes to high-frequency heating in terms of modulated characteristics

In this work,for the first time,we have analyzed and compared the responses of polar mesosphere winter echoes(PMWE)and their summer counterpart,polar mesosphere summer echoes(PMSE),to high-frequency(HF)heating in terms of modulated characteristics(i.e.,backscatter intensity reduction,recovery,and overshoot).Both PMWE and PMSE observations were from the same site(Tromsφ,Norway;69.6°N,19.2°E)and radar(EISCAT[European Incoherent Scatter Scientific Association]very high frequency,224 MHz).The heating patterns of both PMWE and PMSE were found to be similar;however,PMSE was more greatly affected by HF heating.Polar mesosphere summer echoes showed recovery and overshoot more frequently than did PMWE.In addition,the mean recovery and overshoot of PMSE were greater than those of PMWE.The associated electron temperature enhancement was estimated for both PMWE and PMSE and showed that,compared with PMWE,the electron temperature enhancement was more significant in PMSE.The strong heating effects on PMSE may be due to the considerable increase in electron temperature.

A new model for Doppler shift of C-band echoes backscattered from sea surface

Within the framework of the two-scale scattering model,the Doppler shift of C-band radar return signals from the nonlinear sea surface are numerically evaluated.As an analytical approximation method,the Bragg resonance scattering method cannot accurately describe the backscattering field from sea surface.Therefore,in the twoscale scattering model,more accurate scattering coefficient(the normalized radar cross section,NRCS)evaluated by the C-band dual-polarized(HH/VV)empirical geophysical model function(CSAR model)is employed to replace the traditional Bragg NRCS to weight the Doppler shift.The numerical results indicate that there are obvious differences between the Doppler shift weighted by the CSAR NRCS and that weighted by the traditional Bragg NRCS.The hydrodynamic modulation of the large-scale waves is one of the important factors that affect the difference between the Doppler shift predicted in upwind and downwind directions.If the relaxation rate in the hydrodynamic modulation is set to be the angular frequency of the dominant water waves,the Doppler shift predicted by the numerical method can fit the results of the empirical model(C-band empirical geophysical model function,CDOP)well at moderate wind speed.Under low wind condition,the comparison shows that the empirical CDOP model appears to overestimate the Doppler shift.In order to facilitate the application,at the end of this paper a semi-empirical CSAR-DOP model,which is a polynomial fitting formula,is developed for evaluating the Doppler shift of C-band signals from time varying sea surface.

Penetration characteristics of VLF wave from atmosphere into lower ionosphere

The factors affecting the reflection and transmission coefficient of the ionosphere have been analyzed.These factors include wave frequency,incident angle,geomagnetic inclination,electron density and collision frequency in the ionosphere.The ionosphere refractive index is also analyzed.The ionosphere above 70 km is considered to be homogeneous and anisotropic,and the reflection and transmission coefficient matrix is calculated using matrix method.Simultaneously the Booker quartic equation is solved to get the refractive index in the ionosphere.The results show that when the wave frequency is higher,it is easier to penetrate into the ionosphere from its bottom boundary and the propagation attenuation in the ionosphere is smaller.TE（traverse electric） wave and TM（traverse magnetic） wave can both penetrate into the ionosphere with a small incident angle,while TE wave can hardly transmit into the ionosphere when the incident angle is large.The transmission coefficient decreases as the geomagnetic inclination increases.TE and TM wave cannot penetrate into the ionosphere at magnetic equator.When the electron collision frequency is higher,it is easier for VLF wave to penetrate into the ionosphere and the attenuation of ordinary wave is weaker,which may be caused by the energy transportation between the waves and the particles.The ordinary（O） wave experiences severer attenuation than extraordinary（X） wave,and X wave is a penetration mode whereas O wave is a non-penetration mode in the ionosphere.All the results indicate that VLF wave with higher frequency is easier to penetrate into the ionosphere and to be recorded by the satellites at high latitude.It is hard for ULF and the lower frequency VLF wave to transmit into the ionosphere directly for the severe reflection and attenuation.It may transmit into the ionosphere with a small incident angle due to the nonlinear effect,for example,the interaction between the waves and the particles or cross modulation,and then propagate along the whistle duct with small attenuation.This work may be a preliminary theoretical exploration for the future calculation on the response of ground based VLF artificial transmitter in the ionosphere and further study on the seismic ionosphere coupling model.

Review of ionospheric irregularities and ionospheric electrodynamic coupling in the middle latitude region

This paper briefly reviews ionospheric irregularities that occur in the E and F regions at mid-latitudes. Sporadic E(ES) is a common ionospheric irregularity phenomenon that is first noticed in the E layer. ES mainly appears during daytime in summer hemispheres, and is formed primarily from neutral wind shear in the mesosphere and lower thermosphere(MLT) region. Field-aligned irregularity(FAI) in the E region is also observed by Very High Frequency(VHF) radar in mid-latitude regions. FAI frequently occurs after sunset in summer hemispheres, and spectrum features of E region FAI echoes suggest that type-2 irregularity is dominant in the nighttime ionosphere. A close relationship between ES and E region FAI implies that ES may be a possible source of E region FAI in the nighttime ionosphere. Strong neutral wind shear, steep ES plasma density gradient, and a polarized electric field are the significant factors affecting the formation of E region FAI. At mid-latitudes, joint observational experiments including ionosonde, VHF radar, Global Positioning System(GPS) stations, and all-sky optical images have revealed strong connections across different scales of ionospheric irregularities in the nighttime F region, such as spread F(SF), medium-scale traveling ionospheric disturbances(MSTID), and F region FAI.Observations suggest that different scales of ionospheric irregularities are generally attributed to the Perkins instability and subsequently excited gradient drift instability. Nighttime MSTID can further evolve into small-scale structures through a nonlinear cascade process when a steep plasma density gradient exists at the bottom of the F region. In addition, the effect of ionospheric electrodynamic coupling processes, including ionospheric E-F coupling and inter-hemispheric coupling on the generation of ionospheric irregularities, becomes more prominent due to the significant dip angle and equipotentiality of magnetic field lines in the mid-latitude ionosphere. Polarized electric fields can map to different ionospheric regions and excite plasma instabilities which form ionospheric irregularities. Nevertheless,the mapping efficiency of a polarized electric field depends on the ionospheric background and spatial scale of the field.

A Comparison Study of Sensitivity on PJ and NPS Models in China Seas

Evaporation duct is an ubiquitous natural phenomenon over the ocean and can be diagnosed by evaporation duct model.The model proposed by Paulus and Jeske and another model established by the American naval postgraduate school are the most widely used.They are called PJ model and NPS model,respectively.Two methods are used to investigate the global sensitivity of PJ model and NPS model in China Seas.The first method is based on meteorological and oceanographic observation data in China Seas.Considering the system random error caused by sensor measurement inaccuracies,the mean relative error and mean absolute error are used as criterion for sensitivity analysis.The second method,called Extended Fourier Amplitude Sensitivity Test(EFAST),takes into account the interaction between input parameters and is used for sensitivity analysis.The results show that NPS model is more sensitive to the random errors of sensors than PJ model.The mean relative errors of PJ model and NPS model are 11.43%and 14.81%,respectively.The results of global sensitivity parameter analysis indicate that wind speed is the key factor of PJ model,while all input parameter of NPS model have relatively large total sensitivity index.In addition,sensitivity analysis results confirm that wind speed is one of main driving factors for the formation of evaporation duct.These results are valuable for the selection of diagnosis models for evaporation duct,the evaluation of radio wave propagation in the marine atmospheric surface layer,and the prediction technique of evaporation duct based on numerical weather prediction(NWP)in China seas.

Estimation of sea clutter inherent Doppler spectrum from shipborne S-band radar sea echo

Measurement of shipborne radar sea echo is instrumental in collecting the sea clutter data in open sea areas.However,the ship movement would introduce an extra Doppler component into the spectrum of the sea clutter,so the sea clutter inherent spectrum must be estimated prior to investigating the sea clutter Doppler characteristics from the shipborne radar sea echo.In this paper we show some results about a shipborne sea clutter measurement experiment that was conducted in the South China Sea in a period between 2017 and 2018;abundant clutter data have been collected by using a shipborne S-band clutter measurement radar.To obtain the sea clutter inherent Doppler spectrum from these data,an estimation method,based on the mapping relationship between the shipborne clutter spectrum and the inherent clutter spectrum,is proposed.This method is validated by shipborne clutter data sets under the same measuring conditions except for the ship speed.Using this method,the characteristics of the Doppler spectrum lineshapes in the South China Sea are calculated and analyzed according to different sea states,wave directions,and radar resolutions,which can be instrumental in designing the radar target detection algorithms.

Modified GIT model for predicting wind-speed behavior of low-grazing-angle radar sea clutter

A modified GIT model for describing the variational trend of mean clutter reflectivity as a function of wind speed is proposed. It uses two slope adjustment factors and two critical wind-speed factors to define and adjust the increasing slope of reflectivity with respect to wind speed. In addition, it uses a constant factor to compensate the overall amplitude of clutter reflectivity. The performance of the modified GIT model has been verified on the basis of the L-band low-grazing-angle radar sea clutter data. The results are in good agreement with the experimental data, indicating that the model is more effective in predicting the wind-speed behavior of clutter reflectivity than the conventional GIT model, especially for lower and higher wind speeds. We believe that the proposed model can provide deeper insights into the relationship between radar sea clutter reflectivity and sea state conditions.

Preliminary observation results of the Coherent Beacon Systemonboard the China Seismo-Electromagnetic Satellite-1

This paper reports, for the first time, observation results of the Coherent Beacon System(CBS) onboard the China SeismoElectromagnetic Satellite-1(CSES-1). We describe the CBS, and the Computerized Ionospheric Tomography(CIT) algorithm program is validated by numerical experiment. Two examples are shown, for daytime and nighttime respectively. The Equatorial Ionization Anomaly(EIA) can be seen, and the northern crest core is located at ～20°N in the reconstruction image at 07:28 UTC on 20 July 2018(daytime).Disturbances are shown in the reconstruction image at 18:40 UTC on 13 July 2018(nighttime). We find that beacon measurements are more consistent with ionosonde measurements than model results, by comparing Nm F2 at three sites at Lanzhou, Chongqing, and Kunming; consistency with ionosonde measurements validates beacon measurements. Finally, we have studied Vertical Total Electron Content(VTEC) variations from ground to ～500 km(the height of CSES-1 orbit) and ratios of VTEC between beacon measurements and CODE(Center for Orbit Determination in Europe) data. VTEC variation from ground to ～500 km has a range of 7.2–16.5 TECU for the daytime case and a range of 1.1–1.7 TECU for the nighttime case. The Beacon/CODE ratio of VTEC varies with latitude and time. The mean Beacon/CODE ratio is 0.69 for the daytime case and 0.26 for the nighttime case. The fact that the nighttime case yields lower ratios indicates the higher altitude of the ionosphere during nighttime when the ionosphere is assumed to be a thin layer.

Numerical simulation of high-resolution azimuthal resistivity laterolog response in fractured reservoirs

The high-resolution azimuthal resistivity laterolog response in a fractured formation was numerically simulated using a three-dimensional finite element method. Simulation results show that the azimuthal resistivity is determined by fracture dipping as well as dipping direction, while the amplitude differences between deep and shallow laterolog resistivities are mainly controlled by the former. A linear relationship exists between the corrected apparent conductivities and fracture aperture. With the same fracture aperture, the deep and shallow laterolog resistivities present small values with negative separations for low-angle fractures, while azimuthal resistivities have large variations with positive separations for high-angle fractures that intersect the borehole. For dipping fractures, the variation of the azimuthal resistivity becomes larger when the fracture aperture increases. In addition, for high-angle fractures far from the borehole, a negative separation between the deep and shallow resistivities exists when fracture aperture is large as well as high resistivity contrast exists between bedrock and fracture fluid. The decreasing amplitude of dual laterolog resistivity can indicate the aperture of low-angle fractures, and the variation of the deep azimuthal resistivity can give information of the aperture of high-angle fractures and their position relative to the borehole.

Blow-up for a Class of Degenerate Reaction-diffusion Equation with Nonlocal Source

This paper deals with the properties of the solution to a class of nonlocal degenerate reaction-diffusion equation with nonlocal source,subject to the null Dirichlet boundary condition.We first give sufficient conditions for that the solution exists globally or blows up in the finite time.Then the blow-up time is also given.At last,we obtain a property differing from the local source which the blow-up set is the entire interval.

THE EVIDENCE OF LOSSY WEDGE DIFFRACTION COEFFICIENT IN THEORY

Formulas of diffraction field of lossy wedges with less than 180?wedge angle are derived on the basis of the Fresnel-Kirchhoff wave theory and their numerical results are compared with those of the heuristic lossy wedge diffraction coefficient given by Luebbers (1984), showing good agreement between the two types of numerical results which have different bases in theory. The agreement shows that the lossy wedge diffraction coefficient as an extension of UTD is quite reasonable.

Numerical simulation of chorus-driving acceleration of relativistic electrons at extremely low L-shell during geomagnetic storms

During 2018 major geomagnetic storm,relativistic electron enhancements in extremely low L-shell regions(reaching L∼3)have been reported based on observations of ZH-1 and Van Allen probes satellites,and the storm is highly likely to be accelerated by strong whistler-mode waves occurring near very low L-shell regions where the plasmapause was suppressed.It is very interesting to observe the intense chorus-accelerated electrons locating in such low L-shells and filling into the slot region.In this paper,we further perform numerical simulation by solving the two-dimensional Fokker-Planck equation based on the bounce-averaged diffusion rates.Numerical results demonstrate the evolution processes of the chorus-driven electron flux and confirm the flux enhancement in low pitch angle ranges(20◦-50◦)after the wave-particle interaction for tens of hours.The simulation result is consistent with the observation of potential butterfly pitch angle distributions of relativistic electrons from both ZH-1 and Van Allen probes.

Deep-detection of formation boundary using transient multicomponent electromagnetic logging measurements

The maximum distance at which an electromagnetic(EM)logging while drilling(LWD)tool detects an approaching boundary is defined as the depth of detection(DOD).Ultra-deep detection capability of the transient multicomponent EM logging measurement is investigated here.First,we adopt sine and cosine transform to calculate the transient multicomponent responses.Compared to the cosine transform,sine transform is more accurate in solving late-time responses.Then,a time-domain geosignal is introduced to sense the boundary.Results show that DOD of this transient EM measurement can be up to tens of meters,including directionally sensitivity.Additionally,by studying the decay characteristics of different components with time,cross component is confirmed to decay much faster than the coaxial/coplanar components in the formation coordinate system.A pseudo-inversion is thereby proposed to determine the dip angle of anisotropic stratified formation.Theoretical simulation results indicate that this algebraic method can determine the true dip at some particular moments.It is still stable and valid even when considering random measurement errors.Moreover,we establish the linear relationship between the time at which the half-point of geosignal curve appears and the distance to boundary(DTB),which would assist in the quick determination of DTB.

Analysis of inversion error characteristics of stellar occultation simulation data

Atmospheric stellar occultation observation technology is an advanced space-based detection technology that can measure the vertical distribution of trace gas composition,temperature,and aerosol content in a planet’s atmosphere.In this study,an inversion algorithm of the onion-peeling method was constructed to invert the transmittance obtained from the forward mask.The method used a three-dimensional ray-tracing simulation to obtain the transmission path of the light in the Earth’s atmosphere.The relevant parameters were then combined in the high-resolution transmission molecular absorption(HITRAN)database,and line-by-line integration was performed to calculate the atmospheric transmittance.The transmittance value was then used as an input to calculate the vertical distribution of oxygen molecules when using the single-wavelength inversion of the onion-peeling method.Finally,the oxygen molecule content was compared with the value attained by the Mass Spectrometer and Incoherent Scatter Radar Extended(MSISE00)atmospheric model to determine the relative error of our model.The maximum error was found to be 0.3%,which is low enough to verify the reliability of our algorithm.Using Global-scale Observations of the Limb and Disk(GOLD)measured data to invert the oxygen number density,we calculated its relative deviation from the published result to further verify the algorithm.The inversion result was affected by factors such as prior data,the absorption spectral line type,the ellipticity of the Earth,and the accuracy of the orbit.Analysis of these error-influencing factors showed that the seasons and the Earth’s ellipticity affected the accuracy of the model only 0.001%and could therefore be ignored.However,latitude and solar activity had a greater impact on accuracy,on the order of 0.1%.The absorption line type affected the accuracy of the model by as much as 1%.All three of these factors therefore need to be considered during the inversion process.

Reflection Characteristics of Layered Media in Polar Mesopause Related to Polar Mesosphere Summer Echoes

Polar mesosphere summer echoes （PMSE） are very strong radar echoes from alti- tudes close to the polar summer mesopause. The data from sounding rocket campaigns indicate that the radar signal to noise ratio （SNR）, electron density and dust charge density of polar meso- sphere in summer show obvious layered structure. In this paper the theory of wave propagation in layered media is used to study the reflectance and SNR at each layer in polar mesosphere. The calculated SNR using theory of dusty plasma is found in good agreement with the experimental result, which may imply that the intensity of the radar echoes reflected by the layered structure in polar mesosphere where polar mesosphere summer echoes used to occur can interpret partially the phenomenon of PMSE. In other words, reflection may play an important role in the occurrence of PMSE.

Study on the Relation Between the Reflectivity and Frequency in Dusty Plasma of the Polar Summer Mesopause

Relation between the volume reflectivity and frequency in view of the polar mesosphere summer echoes （PMSE） is studied using the small-scale structure of electron density caused by charged dust particles. A theoretical expression for the radar volume reflectivity is derived, which agrees with the statistical result. Both the theoretical and statistical results are confirmed with the data obtained from simultaneous observations at three frequencies. Hence the small-scale structure caused by charged dust particles may be a useful tool for the study on the generation mechanism of PMSE.

The microwave scattering characteristics of sea ice in the Bohai Sea

Microwave remote sensing has become the primary means for sea-ice research, and has been supported by a great deal of field experiments and theoretical studies regarding sea-ice microwave scattering. However, these studies have been barely carried in the Bohai Sea. The sea-ice microwave scattering mechanism was first developed for the thin sea ice with slight roughness in the Bohai Sea in the winter of 2012, and included the backscattering coefficients which were measured on the different conditions of three bands（L, C and X）, two polarizations（HH and VV）, and incident angles of 20° to 60°, using a ground-based scatterometer and the synchronous physical parameters of the sea-ice temperature, density, thickness, salinity, and so on. The theoretical model of the sea-ice electromagnetic scattering is obtained based on these physical parameters. The research regarding the sea-ice microwave scattering mechanism is carried out through two means, which includes the comparison between the field microwave scattering data and the simulation results of the theoretical model, as well as the feature analysis of the four components of the sea-ice electromagnetic scattering. It is revealed that the sea-ice microwave scattering data and the theoretical simulation results vary in the same trend with the incident angles. Also, there is a visible variant in the sensitivity of every component to the different bands.For example, the C and X bands are sensitive to the top surface, the X band is sensitive to the scatterers, and the L and C bands are sensitive to the bottom surface, and so on. It is suggested that the features of the sea-ice surfaces and scatterers can be retrieved by the further research in the future. This experiment can provide an experimental and theoretical foundation for research regarding the sea-ice microwave scattering characteristics in the Bohai Sea.

A comparison of MLT wind between meteor radar chain data and SDWACCM results

A meteor radar chain located along the 120°E meridian in the Northern Hemisphere from low to middle latitudes provides longterm horizontal wind observations of the mesosphere and lower thermosphere(MLT)region.In this study,we report a seasonal variation and its latitudinal feature in the horizontal mean wind in the MLT region observed by six meteor radar instruments located at Mohe(53.5°N,122.3°E),Beijing(40.3°N,116.2°E),Mengcheng(33.4°N,116.5°E),Wuhan(30.6°N,114.4°E),Kunming(25.6°N,108.3°E),and Fuke(19.5°N,109.1°E)stations.In addition,we compare the wind in the MLT region measured by the meteor radar stations with those simulated by the Whole Atmosphere Community Climate Model(WACCM).In general,the WACCM appears to capture well the seasonal and latitudinal variations in the zonal wind component.In particular,the temporal evolution of the eastward zonal wind maximum shifts from July to May as the latitude decreases.However,the simulated WACCM meridional wind exhibits differences from the meteor radar observations.

CHARACTERISTICS OF CAVITY-BACKED NARROW SLOT ANTENNA IN A LOSSY MEDIUM

The characteristics of a cavity-backed using the moment method and Dyadic Green's funused to convert the double series of the elements inintegration middle value theorem and elliptic integr}the elements. Numerical results show the effects of tthe input impedance, which could be used in designing the Cavity-Backed Slot Antenna