ABSTRACT
Crude oil stabilization is one of the important processes for crude oil treatment, which can reduce the evaporation loss of crude oil and improve the process safety of crude oil storage and transportation. The difference of operating conditions affects the product benefit and energy consumption of crude oil stabilization. In order to improve the overall economic efficiency, this paper used the data interaction function of MATLAB and HYSYS, and the temperature and pressure of the stabilized tower were optimized by using the particle swarm optimization (PSO) algorithm. Meanwhile, combined with the distributed energy system, the mixed-integer linear programming model was developed, and the design strategy of the energy supply equipment was determined by taking the minimum annual equipment cost and fuel cost as the objective function, so as to reduce energy costs and carbon emissions. The method was validated in an oil field in China. The results showed that the optimal operating conditions after optimization were 0.1MPa and 130°C. The energy supply equipment was selected as an electric heater, a CHP engine, and a waste heat boiler. The energy consumption cost at this time is CNY/a, which is 17.65% less than the actual operating cost of the site. At the same time, carbon emissions fell by 7.8%. This method can guide the operation of the oilfield site to a certain extent.
Highlights
A mixed-integer linear programming model for distributed energy systems is established.
Apply distributed energy system to crude oil positive pressure flash stabilization process.
Using particle swarm algorithm based on MATLAB and HYSYS connection function.
The validity of the method is verified by the data of an oilfield in China.
The proposed method reduces energy consumption by 17.65%.
Nomenclature
Sets and indices | = | |
= | All facilities: COG, CHP, GB, EH, WHB | |
= | Energy types: electricity, oil, gas | |
Continuous Parameters | = | |
= | Light hydrocarbon price [CNY/t] | |
= | Natural gas price [CNY/m3] | |
= | Crude oil price [CNY/t] | |
= | The yield of crude oil when evaporation losses are not considered [t/a] | |
= | The evaporation loss rate of crude oil [%] | |
= | Electricity price [CNY/(kWh)] | |
= | The calorific value of crude oil [J/(kg·℃)] | |
= | The density of crude oil [kg/m3] | |
= | The calorific value of natural gas [J/m3] | |
= | The heating furnace efficiency [%] | |
= | The cooler efficiency [%] | |
= | Fixed cost of facility [CNY] | |
= | Linear cost of equipment installation [CNY/kW] | |
= | Price of energy types [CNY/(kWh)] | |
= | Working hours of the year [hour] | |
= | Discount rate [-] | |
= | The service life of facility [years] | |
= | Conversion efficiency of facility [-] | |
= | Electricity demand [kW] | |
= | Heating demand [kW] | |
= | Maximum capacity for facilities [kW] | |
Positive continuous variables | = | |
= | Light hydrocarbon production [t/a] | |
= | Natural gas production [m3/a] | |
= | Stabilised crude oil production [t/a] | |
= | Lifting pump power [kW] | |
= | Lifting pump power [kW] | |
= | The heating furnace boost temperature [℃] | |
= | The cooler lowering temperature [℃] | |
= | Capacity of facility [kW] | |
= | Consumption of energy type [kW] | |
= | Energy input to the facility [kW] | |
= | Grid input power [kW] | |
Binary variables | = | |
= | If the facility is selected = 1. Otherwise, = 0 | |
Acronyms | = | |
PSO | = | Particle Swarm Optimization |
MILP | = | Mixed Integer Linear Programming |
DES | = | Distributed Energy System |
COG | = | Crude Oil-fired Generator |
CHP | = | CHP engine |
GB | = | GAS Boiler |
EH | = | Electric Heater |
WHB | = | Waste Heat Boiler |
Acknowledgements
This work was funded by the National Natural Science Foundation of China (52202405),and Science Foundation of Zhejiang Ocean University (11025092122).
Disclosure statement
No potential conflict of interest was reported by the author(s).