[vc_empty_space][vc_empty_space]
Pressure transient modeling in geothermal reservoir by using Picard-Mclaurin iteration
Singarimbun A.a, Fujimitsu Y.b, Djamal M.a, Andajani R.D.a
a Earth Physics and Complex Systems Research Group, Bandung Institute of Technology, Bandung, 40132, Indonesia
b Laboratory of Geothermics, Department of Earth Resources Engineering, Kyushu University, Nishi-ku, Fukuoka, 819-0395, Japan
[vc_row][vc_column][vc_row_inner][vc_column_inner][vc_separator css=”.vc_custom_1624529070653{padding-top: 30px !important;padding-bottom: 30px !important;}”][/vc_column_inner][/vc_row_inner][vc_row_inner layout=”boxed”][vc_column_inner width=”3/4″ css=”.vc_custom_1624695412187{border-right-width: 1px !important;border-right-color: #dddddd !important;border-right-style: solid !important;border-radius: 1px !important;}”][vc_empty_space][megatron_heading title=”Abstract” size=”size-sm” text_align=”text-left”][vc_column_text]© 2014 Trans Tech Publications, Switzerland.Exploitation of geothermal resource results the decrease of fluid pressure in geothermal reservoir. In production process, the analysis of reservoir condition is made by observing the pressure state. When the fluid is pumped into the reservoir, the value of pressure varies with time which is called as transient state. The plot of pressure versus time will create a curve with a certain slope. From the graph of the pressure the reservoir condition can be analyzed. The solution of the pressure transient is identified as exponential integral equation Ei(x). When the input to the function is really small for example at x < 0.01, the equation will form an asymptotic curve. Analytical solution involves logarithm natural and Euler constant (γ). In this paper we try to approach the solution of exponential integral equation by numerical integration. The objective of this study is to make a numerical model of the pressure change in a geothermal reservoir and to compare the result between numerical method and analytic. There are two methods used in this study, first is Picard- McLaurin iteration to solve the ordinary differential equations (ODE), and the second is trapezoidal integration to calculate the function of Ei(x). The modeling shows that the result of the calculation with the numerical method matched with the analytic with the range of error between 0.0008 to 4.5 % for drawdown test and 0.19 to 7.7 % for buildup test.[/vc_column_text][vc_empty_space][vc_separator css=".vc_custom_1624528584150{padding-top: 25px !important;padding-bottom: 25px !important;}"][vc_empty_space][megatron_heading title="Author keywords" size="size-sm" text_align="text-left"][vc_column_text]Exponential integrals,Geothermal resources,Numerical integrations,Ordinary differential equation (ODE),Picard-McLaurin iteration,Pressure transient,Pressure-transient model,Reservoir conditions[/vc_column_text][vc_empty_space][vc_separator css=".vc_custom_1624528584150{padding-top: 25px !important;padding-bottom: 25px !important;}"][vc_empty_space][megatron_heading title="Indexed keywords" size="size-sm" text_align="text-left"][vc_column_text]Picard-McLaurin iteration,Pressure transient,Trapezoid integration,Well testing[/vc_column_text][vc_empty_space][vc_separator css=".vc_custom_1624528584150{padding-top: 25px !important;padding-bottom: 25px !important;}"][vc_empty_space][megatron_heading title="Funding details" size="size-sm" text_align="text-left"][vc_column_text][/vc_column_text][vc_empty_space][vc_separator css=".vc_custom_1624528584150{padding-top: 25px !important;padding-bottom: 25px !important;}"][vc_empty_space][megatron_heading title="DOI" size="size-sm" text_align="text-left"][vc_column_text]https://doi.org/10.4028/www.scientific.net/AMR.1025-1026.959[/vc_column_text][/vc_column_inner][vc_column_inner width=”1/4″][vc_column_text]Widget Plumx[/vc_column_text][/vc_column_inner][/vc_row_inner][/vc_column][/vc_row][vc_row][vc_column][vc_separator css=”.vc_custom_1624528584150{padding-top: 25px !important;padding-bottom: 25px !important;}”][/vc_column][/vc_row]