Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid...Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.展开更多
The subsidence of the mining area was monitored by analyzing the phase of permanent scatters （PS） which maintained high coherence in magnitude of SAR images.A hew method of spatial unwrapping was presented which use...The subsidence of the mining area was monitored by analyzing the phase of permanent scatters （PS） which maintained high coherence in magnitude of SAR images.A hew method of spatial unwrapping was presented which used the subsidence rates calculated on comer reflector （CR） points as constraints for PS network to perform the spatial unwrapping using the parametric adjustment method.The algorithm achieved the integration of CR data and PSInSAR algorithm.The colliery dense distributed area around Baisha reservoir was chosen as the study area in the experiment.The time series of subsidence from February in 2007 to February in 2010 is successfully inversed by using the periodic function to simulate the linear and nonlinear components of the deformation.The simulation results show that the accuracy can be ± 2.1 mm with the leveling data being used as the external validation data.展开更多
Atmospheric effects on interferometric synthetic aperture radar(InSAR) measurements are quantitatively studied based on a tandem pair of SAR data and a month-long continuous GPS tracking data obtained at six stations....Atmospheric effects on interferometric synthetic aperture radar(InSAR) measurements are quantitatively studied based on a tandem pair of SAR data and a month-long continuous GPS tracking data obtained at six stations. Differential atmospheric signals extracted from the SAR data for two selected areas show apparent power law characteristics. The RMS values of the signals are 2.04 and 3.66 rad respectively for the two areas. These differential delays can potentially cause in the two areas peak-to-peak deformation errors of 3.64 and 6.52cm, respectively, at the 95% confidence level and Gaussian distribution. The respective potential peak-to-peak DEM errors are 123 and 221 m. The GPS tropospheric total zenith delays estimate indicates that a peak-to-peak error of about 7.8cm can potentially be caused in a SAR interferogram with only 1 d interval at the 95% confidence level. The error increases to about 9.6cm for 10 d interval. The potential peak-to-peak DEM and deformation errors estimated from GPS total zenith delay measurements are however quite similar to those estimated from InSAR data. This provides us with a useful tool to pre-estimate the potential atmospheric effects in a SAR interferogram before we order the SAR images. Nevertheless, the results reveal that even in a small area the atmospheric delays can obscure centimetre level ground displacements and introduce a few hundred meters of errors to derived DEM.展开更多
The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of c...The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
It is an attractive method to combine GPS observations with the information from other surveying system to improve the ambiguity resolution. This research is conducted to investigate how to obtain the prior height inf...It is an attractive method to combine GPS observations with the information from other surveying system to improve the ambiguity resolution. This research is conducted to investigate how to obtain the prior height information in bathymetric surveying by GPS positioning and how to use the prior height information and to obtain a robust result. The authors deal with the collection and the description of the prior height and the method using height validation to improve the ambiguity resolution. The principle of the method, the relationships between the height threshold and the ambiguity search space are presented. A method to determine the threshold for the height validation is suggested. The field tests are carried out to show the feasibility of the proposed methods.展开更多
A new seismic ray-tracing method is put forward based on parabolic travel-time interpolation(PTI) method, which is more accurate than the linear travel-time interpolation (LTI) method. Both PTI method and LTI method a...A new seismic ray-tracing method is put forward based on parabolic travel-time interpolation(PTI) method, which is more accurate than the linear travel-time interpolation (LTI) method. Both PTI method and LTI method are used to compute seismic travel-time and ray-path in a 2-D grid cell model. Firstly, some basic concepts are introduced. The calculations of travel-time and ray-path are carried out only at cell boundaries. So, the ray-path is always straight in the same cells with uniform velocity. Two steps are applied in PTI and LTI method, step 1 computes travel-time and step 2 traces ray-path. Then, the derivation of LTI formulas is described. Because of the presence of refraction wave in shot cell, the formula aiming at shot cell is also derived. Finally, PTI method is presented. The calculation of PTI method is more complex than that of LTI method, but the error is limited. The results of numerical model show that PTI method can trace ray-path more accurately and efficiently than LTI method does.展开更多
Based on the hazard development mechanism, a water solution area is closely related to the supporting effect of pressure-bearing water, the relaxing and collapsing effect of orebody interlayer, the collapsing effect o...Based on the hazard development mechanism, a water solution area is closely related to the supporting effect of pressure-bearing water, the relaxing and collapsing effect of orebody interlayer, the collapsing effect of thawless material in orebody, filling effect caused by cubical expansibility of hydrate crystallization and uplifting effect of hard rock layer over cranny belt. The movement and deformation of ground surface caused by underground water solution mining is believed to be much weaker than that caused by well lane mining, which can be predicted by the stochastic medium theory method. On the basis of analysis on the engineering practice of water solution mining, its corresponding parameters can be obtained from the in-site data of the belt water and sand filling mining in engineering analog approach.展开更多
The grey forecasting model has been successfully applied to many fields. However, the precision of GM(1,1) model is not high. In order to remove the seasonal fluctuations in monitoring series before building GM(1,1) m...The grey forecasting model has been successfully applied to many fields. However, the precision of GM(1,1) model is not high. In order to remove the seasonal fluctuations in monitoring series before building GM(1,1) model, the forecasting series of GM(1,1) was built, and an inverse process was used to resume the seasonal fluctuations. Two deseasonalization methods were presented , i.e., seasonal index-based deseasonalization and standard normal distribution-based deseasonalization. They were combined with the GM(1,1) model to form hybrid grey models. A simple but practical method to further improve the forecasting results was also suggested. For comparison, a conventional periodic function model was investigated. The concept and algorithms were tested with four years monthly monitoring data. The results show that on the whole the seasonal index-GM(1,1) model outperform the conventional periodic function model and the conventional periodic function model outperform the SND-GM(1,1) model. The mean Absolute error and mean square error of seasonal index-GM(1,1) are 30.69% and 54.53% smaller than that of conventional periodic function model, respectively. The high accuracy, straightforward and easy implementation natures of the proposed hybrid seasonal index-grey model make it a powerful analysis technique for seasonal monitoring series.展开更多
A new non-invasive blood glucose measuring apparatus (NBGMA) made up of MSP430F149 SCM (single chip micyoco) was developed,which can measure blood glucose level (BGL) frequently,conveniently and painlessly. The hardwa...A new non-invasive blood glucose measuring apparatus (NBGMA) made up of MSP430F149 SCM (single chip micyoco) was developed,which can measure blood glucose level (BGL) frequently,conveniently and painlessly. The hardware and software of this apparatus were designed,and detecting algorithms based on conservation of energy method (COEM) were presented. According to the law of conservation of energy that the energy derived by human body equals energy consumed by metabolism,and the relationship between convection,evaporation,radiation and the BGL was established. The sensor module was designed. 20 healthy volunteers were involved in the clinical experiment. The BGL measured by an automatic biochemical analyzer (ABA) was set as the reference. Regression analysis was performed to compare the conservation of energy method with the biochemical method,using the 20 data points with blood glucose concentrations ranging from 680 to 1 100 mg/L. Reproducibility was measured for healthy fasting volunteers. The results show that the means of BGL detected by NBGMA and ANA are very close to each other,and the difference of standard deviation (SD) is 24.7 mg/L. The correlative coefficient is 0.807. The coefficient of variation (CV) is 4% at 921.6 mg/L. The resultant regression is evaluated by the Clarke error grid analysis (EGA) and all data points are included in the clinically acceptable regions (region A:100%,region B:0%). Accordingly,it is feasible to measure BGL with COEM.展开更多
Three dimensional geophysical models were abstracted and established according to characteristics of oil and gas reservoir.Then direct current fields for different models were simulated with finite element software(fi...Three dimensional geophysical models were abstracted and established according to characteristics of oil and gas reservoir.Then direct current fields for different models were simulated with finite element software(finite element program generator) by hole-to-surface resistivity method.Numerical solution was compared with analytical one for the homogeneity earth model.And a new parameter of deviation rate was proposed by analyzing different plot curves.The results show that the relative error of solution for homogeneity earth model may attain to 0.043%.And deviation rate decreases from 18% to 1% and its anomaly range becomes wide gradually when the depth of oil and gas reservoir increases from 200 to 1 500 m.If resistivity ratio of oil and gas reservoir to sur-rounding rock decreases from 100 to 10 for the resistive oil and gas reservoir,the amplitude attenuation of deviation rate nearly reaches 8%.When there exists stratum above oil and gas reservoir,and influence of resistive stratum may be eliminated or weakened and anomaly of oil and gas reservoir can be strengthened.展开更多
The characteristic of geographic information system（GfS） spatial data operation is that query is much more frequent than insertion and deletion, and a new hybrid spatial clustering method used to build R-tree for GI...The characteristic of geographic information system（GfS） spatial data operation is that query is much more frequent than insertion and deletion, and a new hybrid spatial clustering method used to build R-tree for GIS spatial data was proposed in this paper. According to the aggregation of clustering method, R-tree was used to construct rules and specialty of spatial data. HCR-tree was the R-tree built with HCR algorithm. To test the efficiency of HCR algorithm, it was applied not only to the data organization of static R-tree but also to the nodes splitting of dynamic R-tree. The results show that R-tree with HCR has some advantages such as higher searching efficiency, less disk accesses and so on.展开更多
Due to the complicated background of objectives and speckle noise, it is almost impossible to extract roads directly from original synthetic aperture radar(SAR) images. A method is proposed for extraction of road netw...Due to the complicated background of objectives and speckle noise, it is almost impossible to extract roads directly from original synthetic aperture radar(SAR) images. A method is proposed for extraction of road network from high-resolution SAR image. Firstly, fuzzy C means is used to classify the filtered SAR image unsupervisedly, and the road pixels are isolated from the image to simplify the extraction of road network. Secondly, according to the features of roads and the membership of pixels to roads, a road model is constructed, which can reduce the extraction of road network to searching globally optimization continuous curves which pass some seed points. Finally, regarding the curves as individuals and coding a chromosome using integer code of variance relative to coordinates, the genetic operations are used to search global optimization roads. The experimental results show that the algorithm can effectively extract road network from high-resolution SAR images.展开更多
Numerical calculation for two integral transforms in 2.5-D transient electromagnetic forward is a difficult and key task, namely, the inverse Fourier transform and the inverse Laplace transform. Some effective algorit...Numerical calculation for two integral transforms in 2.5-D transient electromagnetic forward is a difficult and key task, namely, the inverse Fourier transform and the inverse Laplace transform. Some effective algorithms for them were described. Based on the known algorithms in DC resistivity on wave-number distribution and selection, we proposed a principle on how to choose the least wave-number concerning the central-loop transient electromagnetic method. First, observe the behavior of transformation function curve with regard to wave-number in Fourier domain. In the light of its asymptote, ascertain the coverage scope of wave-number. Compared with analytic solution, the least wave-number in Fourier domain can be derived. Furthermore, the Laplace numerical inversion algorithm which needs only a few Laplace variables in pure real domain was also introduced here. The procedure was applied to forward modeling on transient electromagnetic field of a vertical magnetic dipole over uniform half-space to demonstrate them effectiveness and general applicability.展开更多
In allusion to the difficulty of integrating data with different models in integrating spatial information, the characteristics of raster structure, vector structure and mixed model were analyzed, and a hierarchical v...In allusion to the difficulty of integrating data with different models in integrating spatial information, the characteristics of raster structure, vector structure and mixed model were analyzed, and a hierarchical vector-raster integrative full feature model was put forward by integrating the advantage of vector and raster model and using the object-oriented method. The data structures of the four basic features, i.e. point, line, surface and solid, were described. An application was analyzed and described, and the characteristics of this model were described. In this model, all objects in the real world are divided into and described as features with hierarchy, and all the data are organized in vector. This model can describe data based on feature, field, network and other models, and avoid the disadvantage of inability to integrate data based on different models and perform spatial analysis on them in spatial information integration.展开更多
By analyzing the abrasive theory of concrete diamond sawblade, the proposal that the diamond should be selected by its function in cutting concrete is presented. The part of the big grit diamonds cut rock, and the par...By analyzing the abrasive theory of concrete diamond sawblade, the proposal that the diamond should be selected by its function in cutting concrete is presented. The part of the big grit diamonds cut rock, and the part of the small grit diamonds improve the wearability of the matrix. The contrast tests are done with different shapes of sawbaldes in split segment, slant "U" slot segment, sandwich segment, turbo segment and three-slot segment. The special shapes of sawblades can improve the effect of cooling and the removing ability of the rock powder. The data of tests show that the efficiency of cutting and the life of sawblades are improved by designing the diamond prescription and using the especial geometry of segment.展开更多
基金Projects(2006AA06Z105, 2007AA06Z134) supported by the National High-Tech Research and Development Program of ChinaProjects(2007, 2008) supported by China Scholarship Council (CSC)
文摘Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.
基金Project(2006AA12Z156)supported by the National High Technology Research and Development Program of ChinaProjects(40974006,40774003)supported by the National Natural Science Foundation of China+1 种基金Project supported by Construct Program of the Key Discipline in Hunan Province,ChinaProjects(2011GK3106,2011SK3077)supported by Science and Technology Department of Hunan Province,China
文摘The subsidence of the mining area was monitored by analyzing the phase of permanent scatters （PS） which maintained high coherence in magnitude of SAR images.A hew method of spatial unwrapping was presented which used the subsidence rates calculated on comer reflector （CR） points as constraints for PS network to perform the spatial unwrapping using the parametric adjustment method.The algorithm achieved the integration of CR data and PSInSAR algorithm.The colliery dense distributed area around Baisha reservoir was chosen as the study area in the experiment.The time series of subsidence from February in 2007 to February in 2010 is successfully inversed by using the periodic function to simulate the linear and nonlinear components of the deformation.The simulation results show that the accuracy can be ± 2.1 mm with the leveling data being used as the external validation data.
文摘Atmospheric effects on interferometric synthetic aperture radar(InSAR) measurements are quantitatively studied based on a tandem pair of SAR data and a month-long continuous GPS tracking data obtained at six stations. Differential atmospheric signals extracted from the SAR data for two selected areas show apparent power law characteristics. The RMS values of the signals are 2.04 and 3.66 rad respectively for the two areas. These differential delays can potentially cause in the two areas peak-to-peak deformation errors of 3.64 and 6.52cm, respectively, at the 95% confidence level and Gaussian distribution. The respective potential peak-to-peak DEM errors are 123 and 221 m. The GPS tropospheric total zenith delays estimate indicates that a peak-to-peak error of about 7.8cm can potentially be caused in a SAR interferogram with only 1 d interval at the 95% confidence level. The error increases to about 9.6cm for 10 d interval. The potential peak-to-peak DEM and deformation errors estimated from GPS total zenith delay measurements are however quite similar to those estimated from InSAR data. This provides us with a useful tool to pre-estimate the potential atmospheric effects in a SAR interferogram before we order the SAR images. Nevertheless, the results reveal that even in a small area the atmospheric delays can obscure centimetre level ground displacements and introduce a few hundred meters of errors to derived DEM.
文摘The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
文摘It is an attractive method to combine GPS observations with the information from other surveying system to improve the ambiguity resolution. This research is conducted to investigate how to obtain the prior height information in bathymetric surveying by GPS positioning and how to use the prior height information and to obtain a robust result. The authors deal with the collection and the description of the prior height and the method using height validation to improve the ambiguity resolution. The principle of the method, the relationships between the height threshold and the ambiguity search space are presented. A method to determine the threshold for the height validation is suggested. The field tests are carried out to show the feasibility of the proposed methods.
文摘A new seismic ray-tracing method is put forward based on parabolic travel-time interpolation(PTI) method, which is more accurate than the linear travel-time interpolation (LTI) method. Both PTI method and LTI method are used to compute seismic travel-time and ray-path in a 2-D grid cell model. Firstly, some basic concepts are introduced. The calculations of travel-time and ray-path are carried out only at cell boundaries. So, the ray-path is always straight in the same cells with uniform velocity. Two steps are applied in PTI and LTI method, step 1 computes travel-time and step 2 traces ray-path. Then, the derivation of LTI formulas is described. Because of the presence of refraction wave in shot cell, the formula aiming at shot cell is also derived. Finally, PTI method is presented. The calculation of PTI method is more complex than that of LTI method, but the error is limited. The results of numerical model show that PTI method can trace ray-path more accurately and efficiently than LTI method does.
基金Project(40404001) supported by the National Natural Science Foundation of China
文摘Based on the hazard development mechanism, a water solution area is closely related to the supporting effect of pressure-bearing water, the relaxing and collapsing effect of orebody interlayer, the collapsing effect of thawless material in orebody, filling effect caused by cubical expansibility of hydrate crystallization and uplifting effect of hard rock layer over cranny belt. The movement and deformation of ground surface caused by underground water solution mining is believed to be much weaker than that caused by well lane mining, which can be predicted by the stochastic medium theory method. On the basis of analysis on the engineering practice of water solution mining, its corresponding parameters can be obtained from the in-site data of the belt water and sand filling mining in engineering analog approach.
文摘The grey forecasting model has been successfully applied to many fields. However, the precision of GM(1,1) model is not high. In order to remove the seasonal fluctuations in monitoring series before building GM(1,1) model, the forecasting series of GM(1,1) was built, and an inverse process was used to resume the seasonal fluctuations. Two deseasonalization methods were presented , i.e., seasonal index-based deseasonalization and standard normal distribution-based deseasonalization. They were combined with the GM(1,1) model to form hybrid grey models. A simple but practical method to further improve the forecasting results was also suggested. For comparison, a conventional periodic function model was investigated. The concept and algorithms were tested with four years monthly monitoring data. The results show that on the whole the seasonal index-GM(1,1) model outperform the conventional periodic function model and the conventional periodic function model outperform the SND-GM(1,1) model. The mean Absolute error and mean square error of seasonal index-GM(1,1) are 30.69% and 54.53% smaller than that of conventional periodic function model, respectively. The high accuracy, straightforward and easy implementation natures of the proposed hybrid seasonal index-grey model make it a powerful analysis technique for seasonal monitoring series.
基金Project(07JJ6133) supported by the Natural Science Foundation of Hunan Province, China
文摘A new non-invasive blood glucose measuring apparatus (NBGMA) made up of MSP430F149 SCM (single chip micyoco) was developed,which can measure blood glucose level (BGL) frequently,conveniently and painlessly. The hardware and software of this apparatus were designed,and detecting algorithms based on conservation of energy method (COEM) were presented. According to the law of conservation of energy that the energy derived by human body equals energy consumed by metabolism,and the relationship between convection,evaporation,radiation and the BGL was established. The sensor module was designed. 20 healthy volunteers were involved in the clinical experiment. The BGL measured by an automatic biochemical analyzer (ABA) was set as the reference. Regression analysis was performed to compare the conservation of energy method with the biochemical method,using the 20 data points with blood glucose concentrations ranging from 680 to 1 100 mg/L. Reproducibility was measured for healthy fasting volunteers. The results show that the means of BGL detected by NBGMA and ANA are very close to each other,and the difference of standard deviation (SD) is 24.7 mg/L. The correlative coefficient is 0.807. The coefficient of variation (CV) is 4% at 921.6 mg/L. The resultant regression is evaluated by the Clarke error grid analysis (EGA) and all data points are included in the clinically acceptable regions (region A:100%,region B:0%). Accordingly,it is feasible to measure BGL with COEM.
基金Projects(2006AA06Z105,2007AA06Z134) supported by the National High-Tech Research and Development Program of China
文摘Three dimensional geophysical models were abstracted and established according to characteristics of oil and gas reservoir.Then direct current fields for different models were simulated with finite element software(finite element program generator) by hole-to-surface resistivity method.Numerical solution was compared with analytical one for the homogeneity earth model.And a new parameter of deviation rate was proposed by analyzing different plot curves.The results show that the relative error of solution for homogeneity earth model may attain to 0.043%.And deviation rate decreases from 18% to 1% and its anomaly range becomes wide gradually when the depth of oil and gas reservoir increases from 200 to 1 500 m.If resistivity ratio of oil and gas reservoir to sur-rounding rock decreases from 100 to 10 for the resistive oil and gas reservoir,the amplitude attenuation of deviation rate nearly reaches 8%.When there exists stratum above oil and gas reservoir,and influence of resistive stratum may be eliminated or weakened and anomaly of oil and gas reservoir can be strengthened.
文摘The characteristic of geographic information system（GfS） spatial data operation is that query is much more frequent than insertion and deletion, and a new hybrid spatial clustering method used to build R-tree for GIS spatial data was proposed in this paper. According to the aggregation of clustering method, R-tree was used to construct rules and specialty of spatial data. HCR-tree was the R-tree built with HCR algorithm. To test the efficiency of HCR algorithm, it was applied not only to the data organization of static R-tree but also to the nodes splitting of dynamic R-tree. The results show that R-tree with HCR has some advantages such as higher searching efficiency, less disk accesses and so on.
文摘Due to the complicated background of objectives and speckle noise, it is almost impossible to extract roads directly from original synthetic aperture radar(SAR) images. A method is proposed for extraction of road network from high-resolution SAR image. Firstly, fuzzy C means is used to classify the filtered SAR image unsupervisedly, and the road pixels are isolated from the image to simplify the extraction of road network. Secondly, according to the features of roads and the membership of pixels to roads, a road model is constructed, which can reduce the extraction of road network to searching globally optimization continuous curves which pass some seed points. Finally, regarding the curves as individuals and coding a chromosome using integer code of variance relative to coordinates, the genetic operations are used to search global optimization roads. The experimental results show that the algorithm can effectively extract road network from high-resolution SAR images.
基金Project(40344022) supported by the National Natural Science Foundation of China
文摘Numerical calculation for two integral transforms in 2.5-D transient electromagnetic forward is a difficult and key task, namely, the inverse Fourier transform and the inverse Laplace transform. Some effective algorithms for them were described. Based on the known algorithms in DC resistivity on wave-number distribution and selection, we proposed a principle on how to choose the least wave-number concerning the central-loop transient electromagnetic method. First, observe the behavior of transformation function curve with regard to wave-number in Fourier domain. In the light of its asymptote, ascertain the coverage scope of wave-number. Compared with analytic solution, the least wave-number in Fourier domain can be derived. Furthermore, the Laplace numerical inversion algorithm which needs only a few Laplace variables in pure real domain was also introduced here. The procedure was applied to forward modeling on transient electromagnetic field of a vertical magnetic dipole over uniform half-space to demonstrate them effectiveness and general applicability.
基金Project (40473029) supported bythe National Natural Science Foundation of China project (04JJ3046) supported bytheNatural Science Foundation of Hunan Province , China
文摘In allusion to the difficulty of integrating data with different models in integrating spatial information, the characteristics of raster structure, vector structure and mixed model were analyzed, and a hierarchical vector-raster integrative full feature model was put forward by integrating the advantage of vector and raster model and using the object-oriented method. The data structures of the four basic features, i.e. point, line, surface and solid, were described. An application was analyzed and described, and the characteristics of this model were described. In this model, all objects in the real world are divided into and described as features with hierarchy, and all the data are organized in vector. This model can describe data based on feature, field, network and other models, and avoid the disadvantage of inability to integrate data based on different models and perform spatial analysis on them in spatial information integration.
文摘By analyzing the abrasive theory of concrete diamond sawblade, the proposal that the diamond should be selected by its function in cutting concrete is presented. The part of the big grit diamonds cut rock, and the part of the small grit diamonds improve the wearability of the matrix. The contrast tests are done with different shapes of sawbaldes in split segment, slant "U" slot segment, sandwich segment, turbo segment and three-slot segment. The special shapes of sawblades can improve the effect of cooling and the removing ability of the rock powder. The data of tests show that the efficiency of cutting and the life of sawblades are improved by designing the diamond prescription and using the especial geometry of segment.
