Esempio n. 1
0
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)
Esempio n. 2
0
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)