def pressure_vary_time(): params = Params() depth = params.params['well']['casing1']['depth'] p = params.params # params.set_production_rate(200) oil_temp = OilTemp(p) oil_temp.load_params() oil_temp.run() # oil_temp.plot() print(oil_temp.params['etc']['t'], oil_temp.temps_in_C[-1]) annular_temp = AnnularTemp(p, oil_temp.temps_in_K, oil_temp.zindex) annular_temp.run() # annular_temp.plot() t = np.concatenate((np.arange(0, 100, 10), np.arange(100, 600, 50))) # t = np.arange(0, 100, 1 / 24) t = [0, 1, 5, 20, 100, 200] t = np.arange(0, 25, 1) r_c = [] r_b = [] for i in t: params.set_time_day(i / 24) ot = OilTemp(params.params) # print(params.params['thermal']['W']) ot.load_params() ot.run() at = AnnularTemp(params.params, ot.temps_in_K, ot.zindex) at.run() p_b = Pressure(params.params, at.temps_B_in_C, at.zindex_B) p_c = Pressure(params.params, at.temps_C_in_C, at.zindex_C) # r r_b.append(p_b.pressure_delta) r_c.append(p_c.pressure_delta) # fig, ax = plt.subplots(1, 1) # # ax.set_yticks(np.arange(50, 100, 10)) # ax.plot(t, r) # fig.show() # export_to_excel(t, r_b, r_c)
def temp_dis_vary_time(): params = Params() depth = params.params['well']['casing1']['depth'] p = params.params # params.set_production_rate(200) oil_temp = OilTemp(p) oil_temp.load_params() oil_temp.run() # oil_temp.plot() print(oil_temp.params['etc']['t'], oil_temp.temps_in_C[-1]) annular_temp = AnnularTemp(p, oil_temp.temps_in_K, oil_temp.zindex) annular_temp.run() # annular_temp.plot() t = np.concatenate((np.arange(0, 100, 10), np.arange(100, 600, 50))) # t = np.arange(0, 100, 1 / 24) t = [0, 1, 5, 20, 100, 200] r = [] zindex = None for i in t: params.set_time_day(i) ot = OilTemp(params.params) # print(params.params['thermal']['W']) ot.load_params() ot.run() at = AnnularTemp(params.params, ot.temps_in_K, ot.zindex) at.run() # r r.append(ot.temps_in_C) zindex = ot.zindex # fig, ax = plt.subplots(1, 1) # # ax.set_yticks(np.arange(50, 100, 10)) # ax.plot(t, r) # fig.show() # export_to_excel_time(depth - zindex, r)