Gas Load Research Articles

A cross-sectoral analysis of energy shortages in Pakistan: based on supply-driven input-output model

The impacts of energy shortages characterized by regular blackouts, natural gas, and electricity load shedding in Pakistan affected each economic sector, causing an energy-induced crisis and ecological sustainability issues. This study was conducted to reveal the benefits of renewable energy and describe the economic losses associated with electricity unavailability using supply-driven input-output as a price model across 34 sectors. The results revealed that exogenous shocks in electricity prices are responsible for bringing significant fluctuations across the business cycle in the country. Similarly, the overall output of Pakistan’s economy will decrease by 24.89 rupees due to a 1-kilowatt-hour reduction in electricity supply. Moreover, both forward and backward linkages of Pakistan’s economy revealed that higher electricity allocation coefficients pose significant output impacts on most sectors. We conclude that indirect output impacts require due consideration to avoid the underestimation problem due to total electricity shortages. It is recommended that the government provide a social and legal framework to boost the environmental sustainability and economic activities in the textile, oil refining, production of cement, and fertilizer sectors for sustainable economic growth.

Research on energy management of integrated energy system coupled with organic Rankine cycle and power to gas

Integrated energy system could offer multiple products, such as heat and power. However, the supplied energy and thermoelectrical ratio may not conform with the demand and or thermoelectrical ratio of end user. In order to solve the energy mismatch between source side and load side an integrated energy system coupled with organic Rankine cycle and power to gas was constructed, and energy management strategies based on organic Rankine cycle and power to gas were proposed. Considering the economic, energy and environmental performance of the system, the evaluation index and multi-objective optimization model were constructed, respectively. Then, a case study of an industrial park including factory, dormitory and commercial buildings was employed and the hourly heating, electrical and gas loads were simulated. Finally, the capacity configuration optimization, system performance of integrated energy system under four strategies: without organic Rankine cycle and power to gas, with power to gas, with organic Rankine cycle and with organic Rankine cycle and power to gas were compared. Simultaneously, the suitable integrated energy system configuration between different thermoelectrical ratio range of end user was analyzed. Results showed that when organic Rankine cycle and power to gas were added to integrated energy system, the output range of thermoelectrical ratio could be widened. When thermoelectrical ratio was less than 0.42, the comprehensive performance of integrated energy system coupled with organic Rankine cycle was the best. When thermoelectrical ratio was greater than 3.89, the comprehensive performance of integrated energy system coupled with power to gas was the best. Integrated energy system coupled with organic Rankine cycle and power to gas presented better performances if the thermoelectrical ratio was between 0.43 and 3.21.

Natural gas demand response strategy considering user satisfaction and load volatility under dynamic pricing

Under the background that natural gas is continuously consumed as a non-renewable energy source and the growing demand for natural gas in human society, the contradiction of supply and demand between natural gas utilization and social development has become increasingly prominent. Demand response (DR) in the system helps to address this issue. Currently, there is a lack of diverse DR for natural gas systems. However, DR not only maintains the balance between supply and demand in the gas system, but also relieve peak shaving pressure in natural gas systems and instruct users to save and reduce energy. Therefore, based on the price elasticity theory and user satisfaction, this study focuses on the purpose of natural gas users and gas suppliers to participate in DR, and constructs a natural gas DR model from the perspective of dynamic pricing. An effective natural gas DR strategy is proposed by classifying peak, plain and valley of natural gas load by K-means algorithm, and solving the DR model by multi-objective particle swarm optimization (MOPSO). The DR strategy proposed in this study is evaluated in terms of both non-heating and heating periods through actual examples. The results show that the DR strategy proposed in this study reduces the daily load demand of natural gas, reduces the maximum peak-to-valley difference of daily load, slows down the daily load volatility, improves user satisfaction, maximizes the interests of users and gas suppliers, and effectively balances the system demand.

Collaborative optimization method of the regional integrated energy system for energy cascade utilization

Energy is an important resource, and its main reserves are obviously insufficient compared with consumption. The energy structure is mainly dominated by fossil energy. The pollution emissions from energy utilization are becoming more and more serious. The energy problem has seriously affected the development of new urbanization. This paper studies the optimal scheduling problem of typical urban integrated energy systems participating in different demand response methods. Based on the traditional power demand response model, the demand response model is constructed according to the equivalent thermal parameters. Its model can optimize the electricity and gas load curves and realize the balance of electricity, heat, and gas in an economical and environmentally friendly way. It can also reduce the total cost of the economic operation of the system and reduce the overall carbon emissions of the system. Through the optimization and regulation of multi-type “source-storage-load”, multi-time scale, and multi-objective joint for the distributed thermal power, the matching of thermal and electric hybrid energy and the peak and valley of the power grid can be realized.

Natural Gas Load forecasting based on the fusion of VMD-SE and CNN-BiLSTM-Attention

The rational dispatch and steady supply of natural gas resources can be effectively supported by accurate forecasting of the natural gas load. The accuracy of natural gas load forecasting is, however, low due to the fundamental characteristics of the non-stationary and non-linear natural gas load. Based on this, this paper proposes a new natural gas load forecasting model based on the fusion of VMD-SE and CNN-BiLSTM-Attention. First, the average instantaneous frequency method is used to determine the optimal K value of VMD, and then the sample entropy method is introduced to reconstruct the sub-sequence of VMD decomposition, which effectively reduces the computational overhead of the algorithm. Then, using CNN to mine the level information of multi-dimensional time series variables, using BiLSTM to mine the time series feature information, and introducing the Attention mechanism to give corresponding weight to different feature factors, highlighting the influence of key factors, the anticipated values of each reconstructed sub-series are superimposed to provide the final predicted value of the natural gas load. The model proposed in this paper has lower error than other forecasting models, can effectively improve the precision of natural gas load forecasting, and has a better forecasting effect, according to the results of establishing various forecasting models using the natural gas load data of one region of China as an example.

Investigation of methane mass transfer and sorption in Marcellus shale under variable net-stress

Natural gas in shale exists as free and adsorbed gas, subject to prevailing pore pressures and stress conditions. Accordingly, to accurately estimate/predict the shale gas recovery potential, a central requirement is to represent gas transport and sorption behavior under varying stress conditions. The objective of this work is to facilitate the interpretation of laboratory-scale experiments, at relevant conditions, in an attempt to bridge the gap in scales between laboratory- and field-scale observations.We have conducted a series of high-pressure experiments on a full-diameter core sample from the Marcellus shale. These include gas loading (pressure-decay) and depletion (production) experiments with pure methane (CH4) at variable stress conditions to characterize transport and sorption behavior under reservoir-relevant conditions. We have formulated and applied a novel integral model for mass transfer and storage in multi-porosity shale systems that allows us to effectively investigate transport and sorption phenomena: We delineate gas transport by interpreting helium (He) pressure-decay experiments and demonstrate how to use the information gained to calculate the relevant transport coefficients of CH4 and other gases. A separate measurement of the CH4 sorption isotherm on a smaller sample (a shale cube) was interpreted and combined with the transport description to predict the production behavior of CH4 from the experiments with the full-diameter core.Our experiments demonstrate that the representation of sorption hysteresis is crucial for predicting and guiding shale gas production: At the end of both gas production experiments, approximately 20% of the initial gas in place remained in the core. Without accounting for sorption hysteresis, our modeling demonstrates that the CH4 production could be overestimated by 10%. We demonstrate that our integral, triple-porosity model provides an effective approach for the interpretation and prediction of gas transport and sorption behavior during loading and production experiments on shale cores under variable net-stress conditions.In summary, our work combines measurements and modeling of mass transfer and sorption in shales at different scales to validate a characterization approach that facilitates an improved understanding of shale gas production. Furthermore, the triple-porosity model utilized in our work defines a potential pathway for the translation of laboratory-scale experimentation to larger-scale applications.

An integrated demand response dispatch strategy for low-carbon energy supply park considering electricity–hydrogen–carbon coordination

Driven by the goal of ‘carbon peak, carbon neutrality’, an integrated demand response strategy for integrated electricity–hydrogen energy systems is proposed for low-carbon energy supply parks considering the multi-level and multi-energy characteristics of campus-based microgrids. Firstly, considering the spatial and temporal complementary nature of wind and photovoltaic generation and energy utilization, the energy flow framework of the park is built based on the electricity and hydrogen energy carriers. Clean energy is employed as the main energy supply, and power, heat, cooling, and gas loads are considered energy consumption. Secondly, the operation mechanism of coupled hydrogen storage, hydrogen fuel cell, and carbon capture equipment is analyzed in the two-stage power-to-gas conversion process. Thirdly, considering the operating costs and environmental costs of the park, an integrated demand response dispatch model is constructed for the coupled electricity–hydrogen–carbon system while satisfying the system equipment constraints, network constraints and energy balance constraints of the park system. Finally, Case study in an energy supply park system is implemented. The dispatch results of the integrated demand response with customer participation in the conventional, electricity–hydrogen and electricity–hydrogen–carbon modes are compared to verify the effectiveness of the proposed strategy in renewable accommodation, environmental protection, and economic benefits.

Robust synthesis of amine-based natural gas sweetening plants

The natural gas processing and refining industry uses acid gas removal units to improve the quality of their sour gas feed to a certain level. Various disturbances such as sour gas feed conditions or solvent quality can influence the operation of such plants. A framework is developed using linear control tools which can be used to assess the flexibility of the operating condition of a sweetening plant against such disturbances. A quantitative measure of the flexibility and controllability of a design is provided as a systematic tool for the designer to evaluate different process structures or conditions. The framework can also be used to calculate the maximum tolerable disturbances by the process. The framework is illustrated in a plant that removes acid gas from a raw gas stream using MDEA as the solvent. The parameter σ1 as the maximum singular value of the steady state transfer matrix of the plant is considered as the sensitivity of the outlet H2S and CO2 purities against external disturbances where a higher value means more sensitivity. The results show that higher solvent concentrations are preferred considering the robustness of the absorption column against the sour gas feed disturbances and the uncertainties of the solvent properties. For example, at the solvent flow rate of 40,000 kg/h, the value of σ1 is 2.817, 0.736 and 0.148 at the wt% of 20%, 30% and 40%, respectively. Also, it is shown that for a plant designed with the sweet gas H2S purity of 2 ppm, when a solvent with 20 wt% and flow of 88897.3 kg/h is used, the feed quality disturbance of 0.59 kg/h and 505.19 kg/h in the H2S and CO2 mass flows relative to their nominal can be tolerated before reaching the maximum allowable value of 3 ppm. Such calculations at the 45% point with solvent flow of 44015.8 kg/h lead to the maximum tolerable disturbance of 1.79 kg/h and 967.36 kg/h in the H2S and CO2 flows. The rationale behind this behavior can be well understood if the plot of acid gas loading of rich amine is considered. A lower acid gas loading is observed at higher solvent wt% which means that more free solvent molecules are available in the solvent stream to absorb new acid gases. Accordingly, it has a better capability of confronting external disturbances. Also, the effect of the number of stages of the absorption column on the maximum tolerable disturbance is studied. It is concluded that the robustness of the process will decrease against both types of disturbances as the number of stages increases and lower disturbances can be handled without H2S outlet constraint violation.

Interdisciplinary design optimization of compressor blades combining low- and high-fidelity models

Multidisciplinary design optimization has great potential to support the turbomachinery development process by improving designs at reduced time and cost. As part of the industrial compressor design process, we seek for a rotor blade geometry that minimizes stresses without impairing the aerodynamic performance. However, the presence of structural mechanics, aerodynamics, and their interdisciplinary coupling poses challenges concerning computational effort and organizational integration. In order to reduce both computation times and the required exchange between disciplinary design teams, we propose an inter- instead of multidisciplinary design optimization approach tailored to the studied optimization problem. This involves a distinction between main and side discipline. The main discipline, structural mechanics, is computed by accurate high-fidelity finite element models. The side discipline, aerodynamics, is represented by efficient low-fidelity models, using Kriging and proper-orthogonal decomposition to approximate constraints and the gas load field as coupling variable. The proposed approach is shown to yield a valid blade design with reasonable computational effort for training the aerodynamic low-fidelity models and significantly reduced optimization times compared to a high-fidelity multidisciplinary design optimization. Especially for expensive side disciplines like aerodynamics, the multi-fidelity interdisciplinary design optimization has the potential to consider the effects of all involved disciplines at little additional cost and organizational complexity, while keeping the focus on the main discipline.

Evaluation of PVC-Type Insulation Foam Material for Cryogenic Applications

With the International Maritime Organization (IMO) reinforcing environmental regulations on the shipbuilding industry, the demand for fuels, such as liquefied natural gas (LNG) and liquefied petroleum gas (LPG), has soared. Therefore, the demand for a Liquefied Gas Carrier for such LNG and LPG also increases. Recently, CCS carrier volume has been increasing, and damage to the lower CCS panel has occurred. To withstand liquefied gas loads, the CCSs should be fabricated using a material with improved mechanical strength and thermal performance compared with the conventional material. This study proposes a polyvinyl chloride (PVC)-type foam as an alternative to commercial polyurethane foam (PUF). The former material functions as both insulation and a support structure primarily for the LNG-carrier CCS. To investigate the effectiveness of the PVC-type foam for a low-temperature liquefied gas storage system, various cryogenic tests, namely tensile, compressive, impact, and thermal conductivity, are conducted. The results illustrate that the PVC-type foam proves stronger than PUF in mechanical performance (compressive, impact) across all temperatures. In the tensile test, there are reductions in strength with PVC-type foam but it meets CCS requirements. Therefore, it can serve as insulation and improve the overall CCS mechanical strength against increased loads under cryogenic temperatures. Additionally, PVC-type foam can serve as an alternative to other materials in various cryogenic applications.

Multi-objective optimization of integrated energy systems based on demand response

ABSTRACT The implementation of demand response (DR) in the operation optimization of the integrated energy system has become one of the research hotspots worldwide with the rapid development of the energy network. Integrated gas-electricity system (IGES) operators can allocate different energy carriers more rationally and decrease energy cost and exergy input considering the DR scheme. However, DR was mainly applied to electric power systems in the past, while natural gas systems were less considered. In this paper, based on the theory of price-elasticity, and simultaneously with the consideration of the DR for end-user (EDU) of electricity and natural gas, the optimization model of the IGES is established, which includes two objective functions, the minimum operating cost of the integrated system and the maximum EDU satisfaction. The constraints include compressors and gas storages in the natural gas system, coal and gas generator units in the electric power system, Pow to gas units, and EDU gas/electric load changes. Subsequently, the proposed model is employed for an integrated 7-node natural gas system and 6-node power system with the sensitivity analysis for operating cost. The results demonstrate that the implementation of DR brings specific economic benefits to the operation of the entire system, the total operating cost is reduced by 16.6%, and the natural gas load is reduced by 89.9%. The implementation of DR balances the 24-hour load, reduces the fluctuation of the EDU load, increases the efficiency of the load shifting, and alleviates the instability of upstream production and transportation.

Performance of biocover in controlling methane emissions from landfill: A decade of full-scale investigation

Methane (CH4) emissions generated from waste management facilities represent a serious global warming concern. The objective of this study was to monitor a fully instrumented biocover capable of abating fugitive CH4 emissions from a closed landfill located at Kudjape, Saaremaa Island, Estonia. This investigation documented the alteration in the emission of CH4, carbon dioxide, and temperature fluctuation for a monitoring period of 10 years after the implementation of biocover. The fine fraction reclaimed from landfill mining along with natural mineral soil and sewage sludge compost in a proportion of 3:1:1 was used as a biocover substrate. CH4 emissions were monitored at the surface and from three locations at a regular depth of 25 cm up to 2 m inside the biocover layer. The measurements recorded from three gas monitoring wells present different scenarios of a typical landfill due to the heterogenous nature of waste and non-homogeneous distribution of landfill gas load to the methane-oxidation layer, namely, a) desired situation where a permanent reduction in CH4 concentration took place over time; b) hot-spot with an extremely high concentration of CH4 albeit consistent decrease over the years; c) a location with minimal initial concentration of CH4. Multivariate analysis showed that the CH4 concentration dynamics and parameters reflecting the CH4 oxidation process activity were different in biocover across the landfill. Based on data analysis results, the CH4 oxidation process stabilizes in biocover in 5–6 years after establishing the cover.

Research on optimal scheduling of integrated energy system in low-carbon parks based on demand response

Abstract The comprehensive energy system of the park is one of the effective ways to solve the problems of low utilization efficiency of comprehensive energy and difficulty in absorbing renewable energy. By coordinating the output of each unit to optimize the scheduling of the system, the operating cost of the system can be reduced to a certain extent, and the space for new energy generation to be connected to the Internet can be increased. Aiming at the characteristics of electric, heat and gas load demand of the integrated energy system of industrial low-carbon parks and the overall needs of low-carbon development, a research on the optimal scheduling of integrated energy systems in low-carbon parks based on demand response is proposed. By analyzing the structure and components of the integrated energy system of the low-carbon park, different types of demand response are modeled. On this basis, based on the comprehensive demand response of electricity and heat, the optimal dispatch model of the integrated energy system of the low-carbon park is constructed, and the model is solved to realize the optimal dispatch of the integrated energy system of the low-carbon park. The experimental resu lts show that the proposed method has lower operating cost of the integrated energy system, better economy of the integrated energy system and can effectively improve the power supply reliability and energy saving rate of the integrated energy system. The average system load rate of the proposed method is up to 98.7%, the average comprehensive energy utilization rate is up to 97.9% and the system operation cost is only 10343.1 yuan.

Numerical Study on Key Factors Affecting Synchronization of Multiple Spark-Gap Switches in Single-Stage LTD

Gas switch jitter in fast linear transformer driver (LTD) has a great influence on output current waveform. Aiming at explaining the phenomenon that the breakdown jitter increases when gas switches are synchronously triggered in a single-stage LTD cavity, key factors that affect gas switch breakdown jitter were analyzed in this article. Based on Monte Carlo analysis, a novel LTD circuit model with switch breakdown delay and jitter function was developed, and its accuracy was verified by experimental results. Key factors that result in breakdown jitter increasing were analyzed by changing the gas switch self-jitter, breakdown sequence, and loop inductance. It is found that the gas switch breakdown process can influence each other when triggered in a single-stage LTD cavity. The simulation results showed that the gas switch self-jitter and load inductance have a great influence on the branches-jitter, where the branches-jitter decreased from 12.01 to 5.81 ns when the load inductance decreased from 5 to 4 nH. The simulation results gave an in-depth understanding of the branches-jitter increasing in a single-stage LTD cavity.

Molecular Dynamics Simulation of Pore-Size Effects on Gas Adsorption Kinetics in Zeolites

AbstractStrong gas-mineral interactions or slow adsorption kinetics require a molecular-level understanding of both adsorption and diffusion for these interactions to be properly described in transport models. In this combined molecular simulation and experimental study, noble gas adsorption and mobility is investigated in two naturally abundant zeolites whose pores are similar in size (clinoptilolite) and greater than (mordenite) the gas diameters. Simulated adsorption isotherms obtained from grand canonical Monte Carlo simulations indicate that both zeolites can accommodate even the largest gas (Rn). However, gas mobility in clinoptilolite is significantly hindered at pore-limiting window sites, as seen from molecular dynamics simulations in both bulk and slab zeolite models. Experimental gas adsorption isotherms for clinoptilolite confirm the presence of a kinetic barrier to Xe uptake, resulting in the unusual property of reverse Kr/Xe selectivity. Finally, a kinetic model is used to fit the simulated gas loading profiles, allowing a comparison of trends in gas diffusivity in the zeolite pores.

Determination of Parameters Affecting the Estimation of Iceberg Draft

Recent offshore oil and gas loading facilities developed in the Arctic area have led to a considerable awareness of the iceberg draft approximation, where deep keel icebergs may gouge the ocean floor, and these submarine infrastructures would be damaged in the shallower waters. Developing reliable solutions to estimate the iceberg draft requires a profound understanding of the problem’s dominant parameters. As such, the dimensionless groups of the parameters affecting the iceberg draft estimation were determined for the first time in the present study. Using the dimensionless groups recognized and the linear regression (LR) analysis, nine LR models (i.e., LR 1 to LR 9) were developed and then validated using a comprehensive dataset, which has been constructed in this study. A sensitivity analysis distinguished the premium LR models and important dimensionless groups. The best LR model, as a function of all dimensionless parameters, was able to estimate the iceberg draft with the highest level of precision and correlation along with the lowest degree of complexity. The ratio of iceberg length to iceberg height as the “iceberg length ratio” and the ratio of iceberg width to iceberg height as the “iceberg width ratio” was detected as the important dimensionless groups in the estimation of the iceberg draft. An uncertainty analysis demonstrated that the best LR model was biased towards underestimating the iceberg drafts. The premium LR model outperformed the previous empirical models. Ultimately, a set of LR-based relationships were derived for estimating the iceberg drafts for practical engineering applications, e.g., the early stages of the iceberg management projects.

Morphologies, kinetics and structures of methane hydrate in the system containing tetrahydrofuran and cyclopentane

The morphologies of methane (CH4) hydrate formation in the system containing different mole ratios of tetrahydrofuran (THF) and cyclopentane (CP) at two different temperatures are systematically studied in this work. The results showed that for all cases conducted at 274.2 K, a hydrate film quickly covers the entire interface between the oil (CP) and the THF solution phase after gas loading, and then a hydrate layer forms, thereby significantly decreasing the following gas-liquids mass transfer efficiency. Raman spectrum shows this hydrate layer is caused mainly by pure THF hydrate. Interestingly, as the temperature increase to 278.2 K, the hydrate layer breaks. This is because, after hydrate layer formation, hydrate still grows at the contact interface of the metal, gas, and liquid phase, and then moves upward. In this process, liquids are gradually transported to the reaction front to form hydrate by capillary force and the hydrate layer’s downward curvature gradually increases as the consumption of liquid until rupture. Besides, typical reaction kinetics shows that there are three stages of hydrate formation in the THF + CP + CH4 system at 278.2 K. (1) hydrate layer formation (pure THF hydrate formation); (2) CP-CH4 hydrate formation; (3) large amount of THF-CH4 hydrate formation (after hydrate layer rupture). Further, this work also confirmed that after the hydrate layer rupture, the CH4 occupancy in 512 cages of sII at the same temperature and promoter’s mole fractions is independent of the type of promoters (CP, THF, and CP + THF) used in the system.

Constraint-based interactive approach for equilibrium of interdependent gas and electricity markets

Locational marginal gas price (LMGP) has been adopted in the modeling of interactive markets operated by electricity system operators (ESOs) and gas system operators (GSOs) because of its capability for preserving privacy in the coordination of electricity and gas systems. In this paper we look at the current price-demand iterative (PDI) algorithm approach to solving the LMGP-based interactive model and find that the PDI algorithm may have some difficulty achieving convergence. To address this convergence issue, a new interactive solution approach is proposed in this paper to find the equilibrium between an ESO and a GSO. First, a critical load interval (CLI) model of the gas system is proposed in order to identify new LMGPs as well as their preconditions, which can characterize the variation relationship of gas loads within several hours. Second, a solution technique is proposed to set the equality constraints from the CLI model and LMGPs as interactive information to search for the equilibrium point. The case study indicates that the proposed interactive approach has a better performance on convergence than the PDI algorithm with fewer convergence numbers when solving the interaction problem of GSOs and ESOs.

Optimal Operation of Integrated Energy Systems Based on Multi-Energy Coupling With Power-to-Gas Technology

Abstract Natural gas and renewable energy sources make up an increasingly large proportion of power generation due to their being environmentally friendly. Because of the intermittent and fluctuating output of renewable energy sources, the emerging technology of power-to-gas (P2G) units is of great significance to alleviate. This paper focuses on the coordinated optimization of the combined gas and electricity network (CGEN) system with bidirectional energy conversion. Gas-fired power generation (GPG) plants and P2G are considered as linkages between the two networks. A unified CGEN mathematical model is established with the minimum operation cost as the objective function. The injection and production process of underground natural gas storage (UNGS) is also taken into consideration that is not available in other literature. Variables such as the output of P2G and gas-fired power generation plants, the supply of UNGS, and the wind curtailment are optimized correspondingly. The improved 24-node power grid (PG) and Belgium 20-node natural gas pipeline network (NGPN) are adopted to test the validity and capability of the proposed model, and then we compare the results CGEN model with the single Belgium NGPN to analyze the impact on the NGPN. Moreover, we adopt a coupled operation system of a 6-node power system and a 7-node NGPN to further analyze the influence of P2G on the CGEN. The results show that compared with the system without P2G, the total operating cost of the system is reduced by 9.39%, the natural gas load shedding is reduced by 26.1%, and the wind curtailment is reduced from 326 MWh to zero.

Optimal Operation of Energy Hub Systems under Resiliency Response Options

The economic and resilient operation of power systems has always been one of the main priorities of energy systems. In spite of improvements in various fields of energy systems, especially power systems, the issue of resilience has become more important. For this purpose, this paper proposes a multiobjective optimization model to improve the economic performance of energy hub systems and improve the resilience of electrical consumers. Also, consumer welfare, which is a function of the energy not supplied index, is maximized over a 24-hour period by considering extreme weather conditions. The ε-constraint method is applied to solve the proposed model by transforming the multiobjective optimization problem into several single-objective optimization problems. The max-min fuzzy method is also used to select the optimal solution among the Pareto solutions. A sample hub system is made up of electrical, thermal, and gas loads, electrical and thermal energy sources, and storage systems employed as a test system. A group of actions is applied to improve the resilience of the system, which may be affected by outages caused by storms under the resilience response program (RRP). The results proved the efficiency of the proposed RRP in improving economics and resilience.