An improved numerical simulation method is presented to calculate the downhole temperature distribution for multiple pay zones in producing oil wells. Based on hydrodynamics and heat transfer theory, a 2-D temperature...An improved numerical simulation method is presented to calculate the downhole temperature distribution for multiple pay zones in producing oil wells. Based on hydrodynamics and heat transfer theory, a 2-D temperature field model in cylindrical coordinates is developed. In the model, we considered general heat conduction as well as the heat convection due to fluid flow from porous formation to the borehole. We also take into account the fluid velocity variation in the wellbore due to multiple pay zones. We present coupled boundary conditions at the interfaces between the wellbore and adjacent formation, the wellbore and pay zone, and the pay zone and adjacent formation. Finally, an alternating direction implicit difference method (ADI) is used to solve the temperature model for the downhole temperature distribution. The comparison of modeled temperature curve with actual temperature log indicates that simulation result is in general quite similar to the actual temperature log. We found that the total production rate, production time, porosity, thickness of pay zones, and geothermal gradient, all have effects on the downhole temperature distribution.展开更多
In this paper, a new forming model of the feed direction burr for drilling process is presented. The feed direction burr formation is experimented and studied. The related theories are analyzed, and the influential ...In this paper, a new forming model of the feed direction burr for drilling process is presented. The feed direction burr formation is experimented and studied. The related theories are analyzed, and the influential factors of the feed direction burrs are pointed out. Furthermore, a certain number of new measures to prevent and decrease the burr in drilling process are advanced.展开更多
Azimuthal electromagnetic(EM)logging while drilling(LWD)has been extensively used in high-angle and horizontal(HA/HZ)wells.However,due to the effects of formation anisotropy,accurate geosteering decision and formation...Azimuthal electromagnetic(EM)logging while drilling(LWD)has been extensively used in high-angle and horizontal(HA/HZ)wells.However,due to the effects of formation anisotropy,accurate geosteering decision and formation evaluations have become increasingly difficult.To quantitatively analyze the effect of anisotropy on tool responses and data processing,this paper investigates the sensitivity of EM LWD measurements to electric anisotropy and inversion accuracy via forward modeling and inversion.First,a sensitivity factor is defined to quantitatively analyze the sensitivity of the magnetic field components and synthetic signals to electric anisotropy.Then,azimuthal EM LWD responses in anisotropic layered formations are simulated,and the sensitivities to formation parameters for compensated and uncompensated tool configurations are comparatively analyzed.Finally,we discuss the effects of the inversion model on bed boundary inversion in anisotropic formations.Numerical simulation and inversion results show that azimuthal EM LWD can be significantly affected by electric anisotropy.Fortunately,by using a symmetrical compensation configuration,the sensitivity of the geosignals to electric anisotropy can be suppressed,and the boundary detection capability can be further enhanced.Anisotropy normally gives rise to separated resistivity curves and abnormal"horns";moreover,complicated nonlinear distortion can also arise in geosignals as the tool approaches a bed boundary.If anisotropy effects are ignored in the inversion process,the estimated bed boundary and formation resistivity are usually unreliable,which may mislead geosteering decisions.展开更多
In order to obtain the distribution rules of in situ stress and mining-induced stress of Beiminghe Iron Mine, the stress relief method by overcoring was used to measure the in situ stress, and the MC type bore-hole st...In order to obtain the distribution rules of in situ stress and mining-induced stress of Beiminghe Iron Mine, the stress relief method by overcoring was used to measure the in situ stress, and the MC type bore-hole stress gauge was adopted to measure the mining-induced stress. In the in situ stress measuring, the technique of improved hollow inclusion cells was adopted, which can realize complete temperature compensation. Based on the measuring results, the distribution model of in situ stress was established and analyzed. The in situ stress measuring result shows that the maximum horizontal stress is 1.75-2.45 times of vertical stress and almost 1.83 times of the minimum horizontal stress in this mineral field. And the mining-induced stress measuring result shows that, according to the magnitude of front abutment pressure the stress region can be separated into stress-relaxed area, stress- concentrated area and initial stress area. At the -50 m mining level of this mine, the range of stress-relaxed area is 0-3 m before mining face; the range of stress-concentrated area is 3-55 m before mining face, and the maximum mining-induced stress is 16.5-17.5 MPa, which is 15-20 m from the mining face. The coefficient of stress concentration is 1.85.展开更多
This paper presents a nonlinear finite element method to resolve the problem of the nonlinear contact between the drillstring and hole wall by using a Multi-directional Contact Gap Element (MCGE) contacting at appropr...This paper presents a nonlinear finite element method to resolve the problem of the nonlinear contact between the drillstring and hole wall by using a Multi-directional Contact Gap Element (MCGE) contacting at appropriate positions in each beam element. The method was successfully applied to the Daqing Oil Field GP1 well. It was shown that the drillstring's contact resistance at any well depth could be obtained by calculations and that as the error in the calculation of the hole top load is below 10%, the calculation result can provide theoretical basis for the design and operation of drillstrings.展开更多
According to the working environment of the polar glacier core drilling, a separatory funnel capillary viscometer which can be used for testing the viscosity of liquids under ultra-low temperature conditions was de- s...According to the working environment of the polar glacier core drilling, a separatory funnel capillary viscometer which can be used for testing the viscosity of liquids under ultra-low temperature conditions was de- signed. This viscometer has a simple structure and it is easy to operate, which can meet the testing requirements of different temperature conditions. The viscosity of the dimethyl silicone oil KF-96L-2.0cs was measured under different temperatures using this designed viseometer, and it is found that the viscometer coefficient K changes linearly with temperature. This testing method has relatively high test accuracy and its relative error is less than 4%, which can be used to test the viscosity of the different liquids in ultra-low temperature conditions.展开更多
Rock drillability reflects the drill bit fragments rock hardly or easily. At present, rock drillability classification indexes have rock single axle compressive strength, point load intensity, fracture stress during c...Rock drillability reflects the drill bit fragments rock hardly or easily. At present, rock drillability classification indexes have rock single axle compressive strength, point load intensity, fracture stress during chiseling, drill speed, chiseling specific work, acoustic parameter, cutting magnitude and so on. Every index reflects rock drillability but isnt overall. It is feasible that using many indexes of fuzzy mathematics method etc. to evaluate rock drillability.展开更多
We present a 3D approach to numerical modeling of the borehole-surface electromagnetic (BSEM) method. The 3D electromagnetic response created by a vertical line current source in a layered medium is modeled using th...We present a 3D approach to numerical modeling of the borehole-surface electromagnetic (BSEM) method. The 3D electromagnetic response created by a vertical line current source in a layered medium is modeled using the 3D integral equation method. The modeling results are consistent with analytical solutions. 3D Born approximation inversion of BSEM data is also conducted for reservoir delineation. The inversion method is verified by a synthetic reservoir model.展开更多
Based on the kinematic features of the Planetary Cutter Head (PCH),this paper deals with specificity of rock breakage with PCH, and the basic laws of this process are presented. The study shows that the axial force ac...Based on the kinematic features of the Planetary Cutter Head (PCH),this paper deals with specificity of rock breakage with PCH, and the basic laws of this process are presented. The study shows that the axial force acting on the drilling face by a single cutter and the uniaxial com-pressive strength of the drilled rock have the marked influence on the penetrating depth of the cutter ,and the breakage of the drilling face with the PCH has two forms,the cyclical and the noncyclical. The major reason why the PCH has higher drilling efficiency than other cutter heads is that the PCH can make more free faces on the drilling face.展开更多
In this study,the cylindrical finite-volume method(FVM)is advanced for the efficient and high-precision simulation of the logging while drilling(LWD)orthogonal azimuth electromagnetic tool(OAEMT)response in a three-di...In this study,the cylindrical finite-volume method(FVM)is advanced for the efficient and high-precision simulation of the logging while drilling(LWD)orthogonal azimuth electromagnetic tool(OAEMT)response in a three-dimensional(3 D)anisotropic formation.To overcome the ill-condition and convergence problems arising from the low induction number,Maxwell’s equations are reformulated into a mixed Helmholtz equation for the coupled potentials in a cylindrical coordinate system.The electrical fi eld continuation method is applied to approximate the perfectly electrical conducting(PEC)boundary condition,to improve the discretization accuracy of the Helmholtz equation on the surface of metal mandrels.On the base,the 3 D FVM on Lebedev’s staggered grids in the cylindrical coordinates is employed to discretize the mixed equations to ensure good conformity with typical well-logging tool geometries.The equivalent conductivity in a non-uniform element is determined by a standardization technique.The direct solver,PARDISO,is applied to efficiently solve the sparse linear equation systems for the multi-transmitter problem.To reduce the number of calls to PARDISO,the whole computational domain is divided into small windows that contain multiple measuring points.The electromagnetic(EM)solutions produced by all the transmitters per window are simultaneously solved because the discrete matrix,relevant to all the transmitters in the same window,is changed.Finally,the 3 D FVM is validated against the numerical mode matching method(NMM),and the characteristics of both the coaxial and coplanar responses of the EM field tool are investigated using the numerical results.展开更多
Based on the action mechanism of linear shaped charge( LSC ), penetration performance of LSC on rock was studied. The optimal standoff and the vertex angle of LSC were studied and determined by lab experiments. Thro...Based on the action mechanism of linear shaped charge( LSC ), penetration performance of LSC on rock was studied. The optimal standoff and the vertex angle of LSC were studied and determined by lab experiments. Through cutting sand-cement grout samples, the spacing interval of boreholes can approach 17.5 times of the bore-hole' s diameter, and the result of the directional expansion of crack is satisfactory. The result of field experiment indicates cutting effect is very good, the ruggedness in fracture plane is less than 50 mm, the rate of half-hole marks is nearly 100 % , and the crack inspection shows that there is no damage in the internal of the cutting part. All these suggest that the orientation fracture blasting with LSC is a good means in directional fracture controlled blasting and is worth popularizing widely.展开更多
One of the most important factors controlling the morphology of the modem Huanghe (Yellow) River delta is consolidation settlement, which is impacted by fast deposition, high water content, and low density of seaflo...One of the most important factors controlling the morphology of the modem Huanghe (Yellow) River delta is consolidation settlement, which is impacted by fast deposition, high water content, and low density of seafloor sediment. Consolidation settlement of the Huanghe River subaqueous delta was studied based on field data, laboratory experiments on 12 drill holes, and the one-dimensional consolidation theory. Results show that vertical sediment characteristics varied greatly in the rapidly forming sedimentary bodies of the modem Huanghe River subaqueous delta. Sediments in the upper parts of drill holes were coarser than those in the deeper parts, and other physical and mechanical properties changed accordingly. On the basis of the one-dimensional consolidation theory and drilling depth, the final consolidation settlement of drill holes was between 0.6 m and 2.8 m, and the mean settlement of unit depth was at 1.5-3.5 cm/m. It takes about 15-20 years for the consolidation degree to reach 90% and the average sedimentation rate within the overlying 50 m strata was at 5 cm/a to 12 cm/a. This study helps to forecast the final consolidation settlement and settlement rate of the modem Huanghe River subaqueous delta, which provides key geotechnical information for marine engineers.展开更多
A new numerical algorithm is presented to simulate the explosion reaction process of mine explosives based on the equation of state, the equation of mass conservation and thermodynamics balance equation of explosion p...A new numerical algorithm is presented to simulate the explosion reaction process of mine explosives based on the equation of state, the equation of mass conservation and thermodynamics balance equation of explosion products. With the affection of reversible reaction of explosion products to explosion reaction equations and thermodynamics parameters considered, the computer program has been developed. The computation values show that computer simulation results are identical with the testing ones.展开更多
The concern on formation damage control of high permeability sandstone reservoir has been growing in oil industry in recent years. The invasion of particles and the filtrate of drilling fluid are proven as one of the ...The concern on formation damage control of high permeability sandstone reservoir has been growing in oil industry in recent years. The invasion of particles and the filtrate of drilling fluid are proven as one of the key factors accounting for reservoir damage. Based on the ideal packing theory, the practical software has been developed to optimize the blending proportion of several bridging agents, and the core flooding tests were conducted to evaluate return permeability of core samples contaminated with different drilling fluids. Experimental results show that the ideal packing approach can reduce the dynamic filtration rate, improve the return permeability and drawdown the breakthrough pressure, indicating that this kind of drilling fluids can meet the demands of formation damage control for high permeability sandstone reservoirs. Some applying procedures for formation damage control are also proposed in this paper.展开更多
基金sponsored by the National Nature Science Foundation of China (Grant No. 40830424).
文摘An improved numerical simulation method is presented to calculate the downhole temperature distribution for multiple pay zones in producing oil wells. Based on hydrodynamics and heat transfer theory, a 2-D temperature field model in cylindrical coordinates is developed. In the model, we considered general heat conduction as well as the heat convection due to fluid flow from porous formation to the borehole. We also take into account the fluid velocity variation in the wellbore due to multiple pay zones. We present coupled boundary conditions at the interfaces between the wellbore and adjacent formation, the wellbore and pay zone, and the pay zone and adjacent formation. Finally, an alternating direction implicit difference method (ADI) is used to solve the temperature model for the downhole temperature distribution. The comparison of modeled temperature curve with actual temperature log indicates that simulation result is in general quite similar to the actual temperature log. We found that the total production rate, production time, porosity, thickness of pay zones, and geothermal gradient, all have effects on the downhole temperature distribution.
文摘In this paper, a new forming model of the feed direction burr for drilling process is presented. The feed direction burr formation is experimented and studied. The related theories are analyzed, and the influential factors of the feed direction burrs are pointed out. Furthermore, a certain number of new measures to prevent and decrease the burr in drilling process are advanced.
基金supported by the National Natural Science Foundation of China(No.41674131,No.41974146,and No.41904109)the Shandong Province Postdoctoral Innovation Projects(sdbh20180025)the Fundamental Research Funds for the Central Universities(No.17CX06041)。
文摘Azimuthal electromagnetic(EM)logging while drilling(LWD)has been extensively used in high-angle and horizontal(HA/HZ)wells.However,due to the effects of formation anisotropy,accurate geosteering decision and formation evaluations have become increasingly difficult.To quantitatively analyze the effect of anisotropy on tool responses and data processing,this paper investigates the sensitivity of EM LWD measurements to electric anisotropy and inversion accuracy via forward modeling and inversion.First,a sensitivity factor is defined to quantitatively analyze the sensitivity of the magnetic field components and synthetic signals to electric anisotropy.Then,azimuthal EM LWD responses in anisotropic layered formations are simulated,and the sensitivities to formation parameters for compensated and uncompensated tool configurations are comparatively analyzed.Finally,we discuss the effects of the inversion model on bed boundary inversion in anisotropic formations.Numerical simulation and inversion results show that azimuthal EM LWD can be significantly affected by electric anisotropy.Fortunately,by using a symmetrical compensation configuration,the sensitivity of the geosignals to electric anisotropy can be suppressed,and the boundary detection capability can be further enhanced.Anisotropy normally gives rise to separated resistivity curves and abnormal"horns";moreover,complicated nonlinear distortion can also arise in geosignals as the tool approaches a bed boundary.If anisotropy effects are ignored in the inversion process,the estimated bed boundary and formation resistivity are usually unreliable,which may mislead geosteering decisions.
基金Projects(10702072, 10632100) supported by the National Nature Science Foundation of China
文摘In order to obtain the distribution rules of in situ stress and mining-induced stress of Beiminghe Iron Mine, the stress relief method by overcoring was used to measure the in situ stress, and the MC type bore-hole stress gauge was adopted to measure the mining-induced stress. In the in situ stress measuring, the technique of improved hollow inclusion cells was adopted, which can realize complete temperature compensation. Based on the measuring results, the distribution model of in situ stress was established and analyzed. The in situ stress measuring result shows that the maximum horizontal stress is 1.75-2.45 times of vertical stress and almost 1.83 times of the minimum horizontal stress in this mineral field. And the mining-induced stress measuring result shows that, according to the magnitude of front abutment pressure the stress region can be separated into stress-relaxed area, stress- concentrated area and initial stress area. At the -50 m mining level of this mine, the range of stress-relaxed area is 0-3 m before mining face; the range of stress-concentrated area is 3-55 m before mining face, and the maximum mining-induced stress is 16.5-17.5 MPa, which is 15-20 m from the mining face. The coefficient of stress concentration is 1.85.
文摘This paper presents a nonlinear finite element method to resolve the problem of the nonlinear contact between the drillstring and hole wall by using a Multi-directional Contact Gap Element (MCGE) contacting at appropriate positions in each beam element. The method was successfully applied to the Daqing Oil Field GP1 well. It was shown that the drillstring's contact resistance at any well depth could be obtained by calculations and that as the error in the calculation of the hole top load is below 10%, the calculation result can provide theoretical basis for the design and operation of drillstrings.
基金Supported Project of the Natural Science Foundations of China(No.41106158)
文摘According to the working environment of the polar glacier core drilling, a separatory funnel capillary viscometer which can be used for testing the viscosity of liquids under ultra-low temperature conditions was de- signed. This viscometer has a simple structure and it is easy to operate, which can meet the testing requirements of different temperature conditions. The viscosity of the dimethyl silicone oil KF-96L-2.0cs was measured under different temperatures using this designed viseometer, and it is found that the viscometer coefficient K changes linearly with temperature. This testing method has relatively high test accuracy and its relative error is less than 4%, which can be used to test the viscosity of the different liquids in ultra-low temperature conditions.
文摘Rock drillability reflects the drill bit fragments rock hardly or easily. At present, rock drillability classification indexes have rock single axle compressive strength, point load intensity, fracture stress during chiseling, drill speed, chiseling specific work, acoustic parameter, cutting magnitude and so on. Every index reflects rock drillability but isnt overall. It is feasible that using many indexes of fuzzy mathematics method etc. to evaluate rock drillability.
文摘We present a 3D approach to numerical modeling of the borehole-surface electromagnetic (BSEM) method. The 3D electromagnetic response created by a vertical line current source in a layered medium is modeled using the 3D integral equation method. The modeling results are consistent with analytical solutions. 3D Born approximation inversion of BSEM data is also conducted for reservoir delineation. The inversion method is verified by a synthetic reservoir model.
文摘Based on the kinematic features of the Planetary Cutter Head (PCH),this paper deals with specificity of rock breakage with PCH, and the basic laws of this process are presented. The study shows that the axial force acting on the drilling face by a single cutter and the uniaxial com-pressive strength of the drilled rock have the marked influence on the penetrating depth of the cutter ,and the breakage of the drilling face with the PCH has two forms,the cyclical and the noncyclical. The major reason why the PCH has higher drilling efficiency than other cutter heads is that the PCH can make more free faces on the drilling face.
基金supported jointly by Strategic Pilot Science and Technology Project of Chinese Academy of Sciences (No. XDA14020102)National key research and development plan (No. 2017YFC0601805)+5 种基金National Natural Science Foundation of China (No. 41574110)Youth Foundation of Hebei Educational Committee (No. QN2018217)Hebei Higher Education Teaching Reform Research and Practice(No. 2018GJJG328)Zhangjiakou science and technology bureau(No. 1821011B)Doctoral Fund of Hebei Institute of Architecture and Civil Engineering (No. B-201606)Academic Team Innovation Ability Improvement Project of Hebei Institute of Architecture and Civil Engineering(TD202011)。
文摘In this study,the cylindrical finite-volume method(FVM)is advanced for the efficient and high-precision simulation of the logging while drilling(LWD)orthogonal azimuth electromagnetic tool(OAEMT)response in a three-dimensional(3 D)anisotropic formation.To overcome the ill-condition and convergence problems arising from the low induction number,Maxwell’s equations are reformulated into a mixed Helmholtz equation for the coupled potentials in a cylindrical coordinate system.The electrical fi eld continuation method is applied to approximate the perfectly electrical conducting(PEC)boundary condition,to improve the discretization accuracy of the Helmholtz equation on the surface of metal mandrels.On the base,the 3 D FVM on Lebedev’s staggered grids in the cylindrical coordinates is employed to discretize the mixed equations to ensure good conformity with typical well-logging tool geometries.The equivalent conductivity in a non-uniform element is determined by a standardization technique.The direct solver,PARDISO,is applied to efficiently solve the sparse linear equation systems for the multi-transmitter problem.To reduce the number of calls to PARDISO,the whole computational domain is divided into small windows that contain multiple measuring points.The electromagnetic(EM)solutions produced by all the transmitters per window are simultaneously solved because the discrete matrix,relevant to all the transmitters in the same window,is changed.Finally,the 3 D FVM is validated against the numerical mode matching method(NMM),and the characteristics of both the coaxial and coplanar responses of the EM field tool are investigated using the numerical results.
文摘Based on the action mechanism of linear shaped charge( LSC ), penetration performance of LSC on rock was studied. The optimal standoff and the vertex angle of LSC were studied and determined by lab experiments. Through cutting sand-cement grout samples, the spacing interval of boreholes can approach 17.5 times of the bore-hole' s diameter, and the result of the directional expansion of crack is satisfactory. The result of field experiment indicates cutting effect is very good, the ruggedness in fracture plane is less than 50 mm, the rate of half-hole marks is nearly 100 % , and the crack inspection shows that there is no damage in the internal of the cutting part. All these suggest that the orientation fracture blasting with LSC is a good means in directional fracture controlled blasting and is worth popularizing widely.
基金Supported by the National Science Foundation for Young Scientists of China(No.41206054)the Joint Fund of the National Natural Science Foundation of China and Shandong Province(No.U1606401)the Special Fund of Chinese Central Government for Basic Scientific Research Operations in Commonweal Research Institutes(No.2015G08)
文摘One of the most important factors controlling the morphology of the modem Huanghe (Yellow) River delta is consolidation settlement, which is impacted by fast deposition, high water content, and low density of seafloor sediment. Consolidation settlement of the Huanghe River subaqueous delta was studied based on field data, laboratory experiments on 12 drill holes, and the one-dimensional consolidation theory. Results show that vertical sediment characteristics varied greatly in the rapidly forming sedimentary bodies of the modem Huanghe River subaqueous delta. Sediments in the upper parts of drill holes were coarser than those in the deeper parts, and other physical and mechanical properties changed accordingly. On the basis of the one-dimensional consolidation theory and drilling depth, the final consolidation settlement of drill holes was between 0.6 m and 2.8 m, and the mean settlement of unit depth was at 1.5-3.5 cm/m. It takes about 15-20 years for the consolidation degree to reach 90% and the average sedimentation rate within the overlying 50 m strata was at 5 cm/a to 12 cm/a. This study helps to forecast the final consolidation settlement and settlement rate of the modem Huanghe River subaqueous delta, which provides key geotechnical information for marine engineers.
文摘A new numerical algorithm is presented to simulate the explosion reaction process of mine explosives based on the equation of state, the equation of mass conservation and thermodynamics balance equation of explosion products. With the affection of reversible reaction of explosion products to explosion reaction equations and thermodynamics parameters considered, the computer program has been developed. The computation values show that computer simulation results are identical with the testing ones.
文摘The concern on formation damage control of high permeability sandstone reservoir has been growing in oil industry in recent years. The invasion of particles and the filtrate of drilling fluid are proven as one of the key factors accounting for reservoir damage. Based on the ideal packing theory, the practical software has been developed to optimize the blending proportion of several bridging agents, and the core flooding tests were conducted to evaluate return permeability of core samples contaminated with different drilling fluids. Experimental results show that the ideal packing approach can reduce the dynamic filtration rate, improve the return permeability and drawdown the breakthrough pressure, indicating that this kind of drilling fluids can meet the demands of formation damage control for high permeability sandstone reservoirs. Some applying procedures for formation damage control are also proposed in this paper.