In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid colu...In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns (radius -27 mm). Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole (the difference is 〉30 times), the stability and efficiency of the perfectly matched layer (PML) scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain (FDTD) in the ALWD simulation, including field-splitting PML (SPML), multiaxial PML(M- PML), non-splitting PML (NPML), and complex frequency-shifted PML (CFS-PML). The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor (a0) and the scaling factor (β0) depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.展开更多
Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumb...Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumber-domain wavefield separation are not defined at the same grid point with the staggered-grid finite-difference method for elastic wavefield simulation, we propose the wavenumber-domain interpolation method to estimate the required values at the common grid points prior to the wavenumber-domain true-amplitude wavefield separation. Moreover, numerical experiments show that the wavenumber-domain interpolation method has high interpolation accuracy and the trueamplitude wavefield separation method shows good amplitude preservation. The application of the proposed methodology to elastic reverse-time migration can obtain good amplitudepreserved images even in the case of some velocity error.展开更多
A better understanding of the regional disparity and imbalance characteristics of China's urbanization development is the important premise for constituting correct policy and strategy and promoting the healthy an...A better understanding of the regional disparity and imbalance characteristics of China's urbanization development is the important premise for constituting correct policy and strategy and promoting the healthy and sustainable development of urbanization in the 21st century. The regional differences of China's urbanization level have close relations with natural conditions of landform and climate etc.,the urbanization level reduces with the eleva-tion of topography and decrease of precipitation. According to the statistical data set of ur-banization in 1950-2006,the temporal change course of inter-provincial disparity of Chinese urbanization level since the founding of New China in 1949 was studied,and then the inter-regional and intra-regional disparities of urbanization development were analyzed by the Theil index and its nested decomposition method,to grasp the dynamic change of spatial disparities of China's urbanization level on the whole. Using the imbalance index model,the imbalance status of urban population distribution relative to total population,grain output,total agricultural output value,gross output value of industry,tertiary industrial output value as well as gross regional product was discussed,to hold the balance characteristics of urbanization development relative to the regional development conditions from the macroscopic scales.展开更多
The full-space transient electromagnetic response of water-filled goaves in coal mines were numerically modeled. Traditional numerical modeling methods cannot be used to simulate the underground full-space transient e...The full-space transient electromagnetic response of water-filled goaves in coal mines were numerically modeled. Traditional numerical modeling methods cannot be used to simulate the underground full-space transient electromagnetic field. We used multiple transmitting loops instead of the traditional single transmitting loop to load the transmitting loop into Cartesian grids. We improved the method for calculating the z-component of the magnetic field based on the characteristics of full space. Then, we established the full- space 3D geoelectrical model using geological data for coalmines. In addition, the transient electromagnetic responses of water-filled goaves of variable shape at different locations were simulated by using the finite-difference time-domain (FDTD) method. Moreover, we evaluated the apparent resistivity results. The numerical modeling results suggested that the resistivity differences between the coal seam and its roof and floor greatly affect the distribution of apparent resistivity, resulting in nearly circular contours with the roadway head at the center. The actual distribution of apparent resistivity for different geoelectrical models of water in goaves was consistent with the models. However, when the goal water was located in one side, a false low-resistivity anomaly would appear on the other side owing to the full-space effect but the response was much weaker. Finally, the modeling results were subsequently confirmed by drilling, suggesting that the proposed method was effective.展开更多
The mode characteristics for twin coupled microdisks are investigated by finite-difference time-domain technique. In the coupled microdisks,the same order whispering-gallery(WG) modes can form coupled modes with split...The mode characteristics for twin coupled microdisks are investigated by finite-difference time-domain technique. In the coupled microdisks,the same order whispering-gallery(WG) modes can form coupled modes with split mode wavelengths. We find that the coupled fundamental and first order WG modes can have anticrossing mode coupling as their wavelengths approach the same value in some case,which prevents the cross of the coupled mode wavelengths. The anticrossing mode coupling greatly reduces the coupled mode Q-factor,because the coupled mode field distribution transfers between the fundamental and the first order WG modes.展开更多
This paper utilizes the Theil and decoupling indices to analyze variation in carbon productivity as well as the factors that influence regional carbon productivity in China and proposes carbon emission reduction count...This paper utilizes the Theil and decoupling indices to analyze variation in carbon productivity as well as the factors that influence regional carbon productivity in China and proposes carbon emission reduction countermeasures. The authors conclude that most provinces exhibit year-on-year rising carbon productivity, a trend which decreases moving from east to western China. When applied to carbon productivity, the Theil index presents distinct regional differences. Moreover, the regional variance in carbon productivity is consistently reduced in eastern China and becomes smaller in central China. The difference, however, grows in western China. Carbon productivity grows with the highest speed in central China and the lowest speed in western China. Overall variation in carbon productivity mainly arises from intra-regional difference, whereas inter-regional difference mainly contributed by eastern China. In recent years; both the decoupling index, a dynamic value equal to the rate of change rate in carbon emissions divided by the rate of change in GDP during a given period of time, and carbon productivity vary in different economic development stages. Even if under the same decoupling state, carbon productivity remains different in three regions, i.e., that of the eastern region is higher than the other two regions .展开更多
Modeling technique for electromagnetic fields excited by antennas is an important topic in computational electromagnetics, which is concerned with the numerical solution of Maxwell's equations. In this paper, a no...Modeling technique for electromagnetic fields excited by antennas is an important topic in computational electromagnetics, which is concerned with the numerical solution of Maxwell's equations. In this paper, a novel hybrid technique that combines method of moments(MoM) with finite-difference time-domain(FDTD) method is presented to handle the problem. This approach employed Huygen's principle to realize the hybridization of the two classical numerical algorithms. For wideband electromagnetic data, the interpolation scheme is used in the MoM based on the dyadic Green's function. On the other hand, with the help of equivalence principle, the scattered electric and magnetic fields on the Huygen's surface calculated by MoM are taken as the sources for FDTD. Therefore, the electromagnetic fields in the environment can be obtained by employing finite-difference time-domain method. Finally, numerical results show the validity of the proposed technique by analyzing two canonical samples.展开更多
基金supported by NSFC(No.41174118)one of the major state S&T special projects(No.2008ZX05020-004)+1 种基金a Postdoctoral Fellowship of China(No.2013M530106)China Scholarship Council(No.2010644006)
文摘In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns (radius -27 mm). Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole (the difference is 〉30 times), the stability and efficiency of the perfectly matched layer (PML) scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain (FDTD) in the ALWD simulation, including field-splitting PML (SPML), multiaxial PML(M- PML), non-splitting PML (NPML), and complex frequency-shifted PML (CFS-PML). The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor (a0) and the scaling factor (β0) depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.
基金supported by the National Science Foundation of China(No.41174100)the Large-scale Oil and Gas Field and Coalbed Methane Development Major Projects(No.2011ZX05019-008-08)the China National Petroleum Corporation(No.2014A-3609)
文摘Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumber-domain wavefield separation are not defined at the same grid point with the staggered-grid finite-difference method for elastic wavefield simulation, we propose the wavenumber-domain interpolation method to estimate the required values at the common grid points prior to the wavenumber-domain true-amplitude wavefield separation. Moreover, numerical experiments show that the wavenumber-domain interpolation method has high interpolation accuracy and the trueamplitude wavefield separation method shows good amplitude preservation. The application of the proposed methodology to elastic reverse-time migration can obtain good amplitudepreserved images even in the case of some velocity error.
基金National Science and Technology Supporting Program of "the Eleventh Five-Year Plan",No.2006BAJ05A06 2006BAJ14B03-01Innovation Project of CAS, No.KZCX2-YW-307-02
文摘A better understanding of the regional disparity and imbalance characteristics of China's urbanization development is the important premise for constituting correct policy and strategy and promoting the healthy and sustainable development of urbanization in the 21st century. The regional differences of China's urbanization level have close relations with natural conditions of landform and climate etc.,the urbanization level reduces with the eleva-tion of topography and decrease of precipitation. According to the statistical data set of ur-banization in 1950-2006,the temporal change course of inter-provincial disparity of Chinese urbanization level since the founding of New China in 1949 was studied,and then the inter-regional and intra-regional disparities of urbanization development were analyzed by the Theil index and its nested decomposition method,to grasp the dynamic change of spatial disparities of China's urbanization level on the whole. Using the imbalance index model,the imbalance status of urban population distribution relative to total population,grain output,total agricultural output value,gross output value of industry,tertiary industrial output value as well as gross regional product was discussed,to hold the balance characteristics of urbanization development relative to the regional development conditions from the macroscopic scales.
基金supported by the National Key Scientific Instrument and Equipment Development Project(No.2011YQ03013307)the Project Funded by the Priority Academic Program Development of Jiangsu Higher Education InstitutionsKey Laboratory of Coal Resources Exploration and Comprehensive Utilization,Ministry of Land and Resources
文摘The full-space transient electromagnetic response of water-filled goaves in coal mines were numerically modeled. Traditional numerical modeling methods cannot be used to simulate the underground full-space transient electromagnetic field. We used multiple transmitting loops instead of the traditional single transmitting loop to load the transmitting loop into Cartesian grids. We improved the method for calculating the z-component of the magnetic field based on the characteristics of full space. Then, we established the full- space 3D geoelectrical model using geological data for coalmines. In addition, the transient electromagnetic responses of water-filled goaves of variable shape at different locations were simulated by using the finite-difference time-domain (FDTD) method. Moreover, we evaluated the apparent resistivity results. The numerical modeling results suggested that the resistivity differences between the coal seam and its roof and floor greatly affect the distribution of apparent resistivity, resulting in nearly circular contours with the roadway head at the center. The actual distribution of apparent resistivity for different geoelectrical models of water in goaves was consistent with the models. However, when the goal water was located in one side, a false low-resistivity anomaly would appear on the other side owing to the full-space effect but the response was much weaker. Finally, the modeling results were subsequently confirmed by drilling, suggesting that the proposed method was effective.
文摘The mode characteristics for twin coupled microdisks are investigated by finite-difference time-domain technique. In the coupled microdisks,the same order whispering-gallery(WG) modes can form coupled modes with split mode wavelengths. We find that the coupled fundamental and first order WG modes can have anticrossing mode coupling as their wavelengths approach the same value in some case,which prevents the cross of the coupled mode wavelengths. The anticrossing mode coupling greatly reduces the coupled mode Q-factor,because the coupled mode field distribution transfers between the fundamental and the first order WG modes.
文摘This paper utilizes the Theil and decoupling indices to analyze variation in carbon productivity as well as the factors that influence regional carbon productivity in China and proposes carbon emission reduction countermeasures. The authors conclude that most provinces exhibit year-on-year rising carbon productivity, a trend which decreases moving from east to western China. When applied to carbon productivity, the Theil index presents distinct regional differences. Moreover, the regional variance in carbon productivity is consistently reduced in eastern China and becomes smaller in central China. The difference, however, grows in western China. Carbon productivity grows with the highest speed in central China and the lowest speed in western China. Overall variation in carbon productivity mainly arises from intra-regional difference, whereas inter-regional difference mainly contributed by eastern China. In recent years; both the decoupling index, a dynamic value equal to the rate of change rate in carbon emissions divided by the rate of change in GDP during a given period of time, and carbon productivity vary in different economic development stages. Even if under the same decoupling state, carbon productivity remains different in three regions, i.e., that of the eastern region is higher than the other two regions .
基金Supported in part by China Postdoctoral Science Foundation under Grant No.201M550839in part by the Key Research Program of the Chinese Academy of Sciences under Grant No.KGZD-EW-603
文摘Modeling technique for electromagnetic fields excited by antennas is an important topic in computational electromagnetics, which is concerned with the numerical solution of Maxwell's equations. In this paper, a novel hybrid technique that combines method of moments(MoM) with finite-difference time-domain(FDTD) method is presented to handle the problem. This approach employed Huygen's principle to realize the hybridization of the two classical numerical algorithms. For wideband electromagnetic data, the interpolation scheme is used in the MoM based on the dyadic Green's function. On the other hand, with the help of equivalence principle, the scattered electric and magnetic fields on the Huygen's surface calculated by MoM are taken as the sources for FDTD. Therefore, the electromagnetic fields in the environment can be obtained by employing finite-difference time-domain method. Finally, numerical results show the validity of the proposed technique by analyzing two canonical samples.