Optimal Design and Analysis of Sector-Coupled Energy System in Northeast Japan

: As for research on sector-coupled energy systems, few studies comprehensively deal with energy carriers and energy demand sectors. Moreover, few studies have analyzed energy conversion functions such as Power-to-Gas, Power-to-Heat, and Vehicle-to-Grid on the energy system performance. This study clariﬁes the required renewable resources and costs in the sector-coupled energy system and cost-optimal installed capacity and operation. We formulated an optimization model considering sector coupling and conducted a case study applying the model in the Tohoku region. As a result, due to sector coupling, the total primary energy supply (TPES) is expected to decrease, and system costs are expected to increase from 1.8 to 2.4 times the current level. System costs were minimized when maximizing the use of V2G by electric vehicles and district heating systems (DHS). From the hourly analysis, it becomes clear that the peak cut effect by Power-to-Heat and the peak shift effect by Vehicle-to-Grid result in leveling the output of electrolyzer and fuel synthesizer, which improves the capacity factor reducing capacity addition. Since a large amount of renewable energy is required to realize the designed energy system, it is necessary to reduce the energy demand mainly in the industrial sector. Besides, in order to reduce costs, it is required to utilize electric vehicles by V2G and provide policy support for district heating systems in Japan.


Introduction
The Paris Agreement was adopted in 2015, and it was agreed globally that CO 2 emissions would be reduced to zero by mid-century. It is essential to switch from a conventional fossil-fueled energy system to an innovative energy system with renewable energy resources. Many studies have been performed on integrating large amounts of renewables into energy systems in recent years, and many results have been reported. Lund et al. stated the need for smart energy systems that leverage the synergies of individual sectors to build the sustainable energy system at the lowest cost by consistently capturing energy systems [1]. Connolly et al. presented a scenario for a 100% renewable energy system in Europe in 2050, citing specific steps in terms of scientific and political certainty [2]. Niemi et al. also showed that the share of renewable energy in a region can be increased by considering energy interchange for multiple energy carriers [3]. The studies mentioned above provide suggestions on how to offset the intermittentness of volatile renewable energies and integrate them throughout the energy system. Based on these research results, Thellufsen et al. presented the concepts of "cross-border" and "cross-sector" and verified the effect [4]. As a result, both are important concepts, but it is stated that cross-sectorbased system integration should be prioritized over cross-border-based transmission line enhancements. Syranidou et al. developed the transmission model with variable renewable energy sources (VRES) as main resources and applied it to integrate large-scale variable renewable energy sources in the future European power system [5].
"Cross-sector" is an approach for integrating energy systems by converting electricity to other energy carriers. In recent years, this concept has often been regarded as the term "sector coupling". Sector coupling is defined as the concept that promotes the use of variable renewable energy through system integration between electricity, heat, and transportation.
The coupling of electricity and synthetic fuels using hydrogen as an intermediate is generally called Power-to-Gas (P2G) and Power-to-Liquid (P2L). Various studies have been conducted on evaluating the manufacturing costs of synthetic fuels. Tremel et al. evaluated six types of synthetic fuels from the viewpoints of technology, economy, and acceptability [6]. Gotz et al. summarized the technical and economic barriers of electrolyzers, catalytic methanation, and biological methanation for P2G using renewable energy [7]. Guilera et al. showed that the technical and economic evaluation of SNG could reduce the cost of natural gas to the current cost level in the future [8]. These studies have had significant achievements in research and development by presenting the technical issues of synthetic fuels. However, since the role played by synthetic fuels differs depending on the application area and input conditions of the technology, research considering integration into the energy system is required. For instance, Ikaheimo et al. optimized plant capacity and operation in nine European countries with a model that integrates the power sector with synthetic fuel plants [9]. Varone et al. validated the contribution of P2L and P2G in the German energy system in 2050 in some scenarios of increased renewable energy introduction [10]. Blanco et al. analyzed 80-95% CO 2 reduction scenarios by 2050 (compared to 1990) to identify barriers and drivers of methanation [11]. Nastasi applied a solar Power-to-Gas system to an island energy system [12].
The coupling between electricity and heat is called Power-to-Heat (P2H), and there are many reports on the integration of electrical grid and heating systems. For instance, Bloess et al. simulated P2H and heat storage effects on the power sector by extending the long-term model of power development. As a result, P2H has been shown to increase electricity demand [13]. Obara et al. applied the optimization model that minimizes the system's economic efficiency to Hokkaido. As a result, it became clear that the system used was much cheaper than the storage battery. It was also confirmed that the introduction of renewable energy increased by 21.5% [14]. Ashfaq et al. used the model to examine the possibility of coupling the power grid and heating. The results suggest that current surplus electricity generation is not sufficient to support the heating sector fully, but heat storage can play an important role when renewable energy generation increases by 50% [15].
Vehicle-to-Grid (V2G), the system that utilizes battery electric vehicles as batteries for the electrical grid, is also the coupling between electricity and transportation. Many studies reported the impact of introducing V2G into the power system. For instance, Taljegard et al. simulated how the power supply mix changes with the spread of electric vehicles in Scandinavia and Germany. As a result, it was found that the V2G charging strategy for passenger car EVs smooths the net load curve and almost completely reduces the need for peak power capacity [16]. Child et al. analyzed a 100% renewable energy scenario by using EnergyPLAN. As a result, participation in V2G connectivity promoted a high share of variable RE on a daily and weekly basis [17].
These studies examine individual couplings in the energy system. However, incorporating multiple couplings may cause unanticipated effects because of differences in technical features, costs, and energy qualities. Therefore, research incorporating multiple couplings into the energy system is important for introducing renewable energy into the entire energy system. For instance, Brown et al. comprehensively analyzed cross-border and cross-sector synergies in Europe from electricity, heat, and transportation perspectives. As a result, it was shown that the cost reduction effect of sector coupling becomes weaker as the transmission line is strengthened. However, the energy demands of the industry are not considered in this study [18]. The industry occupies a large share of energy demand and is one of the consumers that has not yet fully introduced renewable energy. Therefore, it is necessary to analyze the amount of renewable energy introduced and the cost of the entire system required by sector coupling in the energy system that also considers industry.  Bailera et al. presented the review of Power-to-X processes as applied to the iron/steel industry [19]. In addition, since sector coupling is a concept for the overall optimization of the energy system, many stakeholders need to cooperate to manage the energy system. It is necessary to consider this cooperation method from the viewpoint of appropriate market management and the priority of technology introduction. However, so far, there has not been sufficient discussion about the optimal operation of each technology in the energy system and the role of individual couplings in Japan.
The purpose of this study is to clarify how much renewable energy is needed and how much it costs in a sector-coupled energy system (including industry) compared to the current situation in this country. We also clarify the cost-optimal installed capacity and operation of each technology. The present study provides useful knowledge for energy planning to realize decarbonization by clarifying the required technologies, their costs, and optimal operation. Among other technologies, sector coupling gives us another solution to lead a zero-carbon society having renewable energy resources. The study focuses on the Tohoku region in Northeast Japan, as it has more abundant renewable energy potential than other regions. Based on actual data on electricity and heat demand, we designed and analyzed a cost-effective energy system to show the benefits of introducing a sector coupling function replacing the conventional energy system. On 26 October 2020, the Prime Minister of Japan announced a carbon-neutral target for Japan by 2050. The newest policy upgraded Japan's previous long-term strategy of reducing emissions by 80% by 2050 to full decarbonization. The approach and results shown here could first contribute to a practical design for zero-carbon Japan toward 2050, making full use of domestic renewable resources.
First, the optimization model for designing sector-coupled energy systems is developed and applied to Japan's particular regional energy system. Then, some scenarios assuming decarbonized energy systems are analyzed using the model. Finally, the optimally designed energy systems' performances are compared, and each technology's optimal operations are explored.
In October 2020, the Government of Japan declared that it would aim to realize a decarbonized society by 2050. According to an IEA report [20], Japan's total primary energy supply (TPES) in 2017 was 18.01 EJ, of which about 93% depended on fossil fuels (including nuclear power). Furthermore, 93% of the total primary energy supply depends on imports from abroad. This is shown in the study on analyzing the energy flow and cash flow diagrams in Japan's energy system [21]. Therefore, expanding the introduction of renewable energy in Japan is a top priority from the environmental, economic, and energy security perspectives.
Furthermore, Japan is the country that emits the fifth largest amount of CO 2 in the world. If Japan manages to reduce most of its emissions, it would significantly impact the global environment and the international community. The implementation of decarbonized society plans by increasing transmission lines based on the cross-border concept is underway by the continued introduction of renewable energy in Japan. On the other hand, there is a description of "sector coupling" in the documents published by the Japanese Government and Ministries, but no detailed vision is drawn, and no research report has been made to support it. This research can be a primary material for examining sector coupling as one direction of Japanese future energy policy.

System Configuration of the Sector-Coupled Energy System
The energy system configuration with sector-coupling in this study is shown in Figure 1. Solar energy, wind energy, hydro energy, run-of-river, geothermal energy, and woody biomass are considered resources. Nuclear energy is not included as no nuclear power plants are in operation in the target area since the Great East Japan Earthquake of 2011. This model consists of five secondary energy carriers: electricity, hydrogen, synthetic fuels, high-temperature heat, and low-temperature heat. The conversion technologies Regarding the inputs and outputs of these conversion technologies, only the conversion of energy is considered, and the input of water to the electrolyzer is not considered in this model. Low-temperature waste heat from 60 to 80 • C is generated by the electrolyzer [22]. This waste heat was difficult to recover in the conventional district heating system (DHS) because of its low temperature. However, based on the concept of 4th generation district heating [23], heat recovery from electrolyzer is possible because the supply temperature can be lowered, and various heat resources can be used. Therefore, in this study, it is assumed that the electrolyzer's waste heat can be recovered. The fuel cell's waste heat can also be retrieved as low-temperature heat. As energy storage technology, stationary battery, hot water tank (HWT), high-pressure tank (HPT), and fuel tanks are handled in the model. If the heat source and heat demand are close to each other, heat storage can be considered, but this model does not reflect the location. Therefore, storage equipment of high-temperature heat is not handled in this research.
Synthetic fuels such as methane, methanol, dimethyl ether (DME), diesel, and gasoline can be produced by reacting hydrogen with CO 2 or CO. In this model, only CO 2 is considered a reactant for fuel synthesis, and they are supplied by direct air capture (DAC), biomass, and gas CHP plants. In the transportation demands, fuel cell vehicles (FCV) and battery electric vehicles (BEV) are taken into account because they do not emit exhaust gases such as CO 2 and NOx when running. The non-manufacturing industry consists of agriculture, forestry, fishery, construction, and mining. Therefore, fuels are mainly consumed by agricultural machines, forestry machines, fishing boats, and construction machines.
The remaining fuels are consumed in the manufacturing industry's heating process. Therefore, only synthetic fuel is supplied to the non-manufacturing industrial and hightemperature heat and synthetic fuels supplied to the industrial manufacturing sector. Local heat pumps are installed in customer's buildings. In the buildings connected to DHS, heat exchangers and radiators are installed.

Optimization Model
This section describes the optimization model for designing sector-coupled energy systems. The proposed optimization model is formulated as a linear programming prob-

Optimization Model
This section describes the optimization model for designing sector-coupled energy systems. The proposed optimization model is formulated as a linear programming problem that minimizes the total annual system cost with each technology's installed capacity and the operations of these technologies as variables. The annual total system cost, which is the objective function, comprises capital, operation and maintenance costs, CO 2 procurement costs for fuel synthesis, costs for transmission between regions, and liquefaction cost of gaseous fuel as shown in Equations (1)- (10).
The first constraints are the hourly balances of electricity. In the electricity sector of region r, it is necessary to balance generation, transmission to other regions, conversion from other carriers, demand, conversion to other carriers, and charging of batteries with a time resolution of 1 h, as shown in Equation (11). The amount of electricity generation is calculated by adding up the product of installed capacity as a decision variable and the standardized output that depends on the time fluctuation of renewable energy resources as a given parameter for each generation technology as shown in Equation (12). The amount of electricity transmission is expressed as the sum of exports to other regions, which are decision variables. Here, the exports to other regions can take a negative value, which is an import (Equation (13)). The amount of conversion from other energy carriers to electricity is the sum of the amount of electricity generated by the fuel cell, biomass, and gas CHP, which are the decision variables (Equation (14)). The total electricity demand includes the amount of electricity required for DAC and methane liquefaction in addition to the local electricity demands (Equation (15)). The local electricity demands are given parameters, and the electricity required for DAC and methane liquefaction is calculated using the CO 2 captured from DAC and the amount of methane liquefied as decision variables. The amount of conversion from electricity to other energy carriers is the sum of the inputs of electricity to electric boilers and heat pumps in the DHS and electrolyzer, which are decision variables (Equation (16)).
Regarding charging, we distinguish between stationary and in-vehicle use. In this optimized model, the charging and discharging of batteries are represented by one variable. Therefore, in Equation (17), there are only two variables related to electricity storage. In the case of a stationary battery, if the variable "chg" is positive, it means charging. In the case of in-vehicle use, if the variable "v2g" is positive, it means discharging.
The amount of electricity stored in the batteries at a specific time is defined by adding the amount of charge to the amount of electricity stored 1 h ago, as shown in Equation (18). The amount of electricity storage is limited by the installed capacity of the stationary battery, which is a decision variable and its depth-of-discharge (DoD), which is assumed to be 80% in this model (Equation (19)). The maximum output of charging or discharging depends on the stationary battery's installed capacity (Equation (20)). The stored electricity in vehicles participating in V2G at time h is defined as shown in Equation (21). The amount of electricity stored in BEVs in the entire area is defined as the value obtained by subtracting the amount of electricity discharged to the grid and the electricity consumption due to running, a given parameter from the amount of electricity stored 1 h ago. The amount of electricity stored in vehicles is restricted by the storage capacity of the battery equipped in the vehicle and state of charge (SOC) (Equation (22)). The discharge to the grid is also limited by the number of inactive BEVs in the area, which is calculated from the given parameters as shown in Equation (24).
− or battery ·Cap battery (r) ≤ chg elect. (r, h) ≤ or battery ·Cap battery (r), Low-temperature heat in DHS also needs to balance generation, conversion from other carriers, demand, and heat storage, as shown in Equation (25). The generation of low-temperature heat is defined as the amount of waste heat recovered from electrolyzers and fuel cells which are decision variables (Equation (26)). The conversion of other energy carriers to low-temperature heat is defined as the sum of the amount of heat supplied by each energy carrier to EB and HP in the DHS, and HOB of biomass and fuels, which are variables multiplied by the respective conversion efficiency (Equation (27)). The energy demand in DHS is defined as a given parameter which is the actual demand divided by the efficiency of the heat exchanger (Equation (28)). The amount of low-temperature heat stored in the hot water tank (HWT) at time h is defined by adding the amount of charge, which is a variable of the amount of low-temperature heat stored 1 h ago, as shown in Equation (30), and is restricted by the installed capacity of HWT as a decision variable (Equation (31)).
For high-temperature heat, heat production by CHP plants, waste heat during fuel synthesis, supply to the manufacturing industry, and heat demand for DAC must be balanced hourly as shown in (32). In the operation of CHP using woody biomass and synthetic methane, the output of electricity and the output of high-temperature heat satisfy (33)- (36).
Hydrogen produced by the electrolyzer is stored in a high-pressure tank (HPT). The amount of hydrogen stored in the HPT is defined by the amount of hydrogen 1 h ago, the amount of hydrogen produced by the electrolyzer, and the amount of hydrogen required for fuel synthesis and FCV driving, which is a given parameter (Equation (37)). As shown in Equation (38), this storage amount should not exceed the installed capacity of HPT as a variable. If the production and consumption of hydrogen are always balanced, HPT is not installed. The electricity consumption of the compressor is not considered in this model. The storage amount of synthetic methane is defined as shown in Equation (39) based on the energy balance of the storage amount 1 h ago, and decision variables (methane production by methanation, consumption at CHP and HOB, supply to the industrial sector, and liquefaction volume) and is restricted by the installed capacity of the existing gas holder, which is a given parameter (Equation (40)). Synthesized methane in a gaseous state can be stored in a liquid state after liquefaction. At this time, electricity corresponding to the amount of liquefaction is consumed. The storage amount of methane in the liquid state is defined by Equation (41) and is restricted by the existing gas holder's installed capacity, which is a given parameter as shown in Equation (42). The storage amount of another synthetic fuel is defined as shown in Equation (43) based on the energy balance of the storage amount 1 h ago and decision variables such as fuel production by fuel synthesis, consumption at HOB, and supply to the industrial sector and is restricted by the installed capacity of existing storage technology, as shown in Equation (44).
The stock of woody biomass starts with the annual potential for use, and the stock decreases with the use of CHP and HOB, as shown in Equations (45)- (47). The storage of CO 2 is defined as in Equation (48) based on the balance of storage 1 h ago, the CO 2 recovery from CHP, and the CO 2 demand for fuel synthesis. The amount of CO 2 stored should not exceed the capacity of the introduced CO 2 tank, as shown in Equation (49).
The fuel consumption of the industrial sector is not captured every hour. In addition, many consumers often have storage facilities. Therefore, in this model, the fuel supplydemand balance in the industrial sector is satisfied in 1 year. Agriculture, forestry and fisheries, and construction are the main consumers in the non-manufacturing industry. Agricultural machinery and construction machinery account for a large proportion, making it difficult to supply energy with high-temperature heat or gas. Therefore, the current city gas is replaced with synthetic methane, gasoline/light oil is replaced with methanol, and LPG is replaced with DME in this model. This is defined by Equation (50). Since fuel consumption in the manufacturing industry is ultimately supplied as process heat, it is considered that most of the fuel is used for combustion, and no distinction is made, as shown in Equation (51).
The installed capacity of generation technology, which is a decision variable, is limited by the introduction potential considering the restrictions on land use, etc., as shown in Equation (52). The outputs from conversion technologies at each time must not exceed the installed capacity (Equation (53)). The transmission line between regions is restricted not to exceed the transmission line's capacity, as shown in Equation (54).

Target Region
In this study, the Tohoku region of Japan is considered a target region because the renewable resources potential is sufficient to give the optimization model feasible solutions. This region is located in the northeastern part of Japan and consists of Aomori, Akita, Iwate, Miyagi, Yamagata, Fukushima, and Niigata prefectures. These seven prefectures are defined as "regions" in the optimization model. The Tohoku region has 79,531 km 2 and a population of 11.28 million. In other words, 9% of the population resides in 21% of Japan's area. The spatial characters of energy supply and demand in the Tohoku region are shown in Figure 2. Figure 2a-c show the estimation results applying methods described later in Sections 2.3.2 and 2.3.3. Energy demand is high in the southern regions such as Miyagi, Niigata, and Fukushima. Renewable energy potential is estimated by considering the capacity factor of each technology. Tohoku region has much potential for offshore wind, particularly in Aomori and Akita. This is because offshore wind power's capacity utilization rate is high in these regions, leading to a reduction in the power generation cost of offshore wind power. The grid is modeled considering the capacity of significant transmission lines [24,25] in the Tohoku region, as shown in Figure 2d. Since the Tohoku region is connected to the other areas by inter-regional transmission lines, it can transmit these regions. However, this study does not consider transmission to the other regions. The transmission cost is assumed to be 4.5 JPY/kWh per link [26].

Estimation of Energy Demands
For estimating the hourly electricity demand in each region, the actual data are referenced [27]. Since these data are the demand for the entire target area, it is prorated by the annual data [28] for each of the seven prefectures. Low-temperature heat demand is defined as the demand for space and water heating in residential and commercial sectors. In this study, the low-temperature heat demand in the target area is estimated based on spatial information analysis. The residential sector's annual heat demand is calculated from households [29] and the annual heat demand per household [30]. That of the commercial sector is estimated using the same method as Fujii et al. [31]. This method estimates annual heat demand from the total floor area [32] and the heat demand intensity of each business.
Furthermore, the hourly heat load is estimated from the above heat demand and the daily heat load patterns [33] for each month. These patterns are categorized by household and business. Some of the estimated heat demand has already been electrified. Therefore, electrified heat demand is deducted from electricity demand. Figure 2c shows the spatial distribution of heat demand in the Tohoku region. When designing power and heat coupling, it is desirable to use existing DHSs to reduce capital costs. However, there are very few DHSs in the target area. In addition, these systems are small-scale systems with a small number of consumers. Therefore, it is necessary to estimate the capital cost of distribution networks if new DHSs are installed in this region. This "distribution capital cost" of DHS is obtained from the heat-demand density and spatial information [34].
Transportation demands are classified into passengers and freights. Each annual mileage is referenced from the statistics on traffic [35]. Annual transportation fuel demand is estimated by multiplying this mileage by each vehicle's tank-to-wheel efficiency. The tank-to-wheel efficiency is summarized in Table 1 [36][37][38]. The hourly energy demand for FCV and BEV is assumed to occur when the energy is charged, not during driving. Energy consumption not related to driving, such as air conditioning in a car, is not considered. The amount of filling per hour is referenced from the database on charging electric vehicles [39].    The fuel demands for industrial use in each region are referenced from statistics on energy demand in each prefecture [28]. In addition to the energy demands mentioned above, there are energy demands for cooking and cooling use in the residential and commercial sectors. These energy demands are partially electrified, but many are still filled with fossil fuels. For exhaustive analysis, it is necessary to consider the electrification of these energy demands. However, it is difficult to estimate its hourly demand accurately, and the proportion of the total consumption is tiny. Therefore, the energy demands that are not electrified are excluded from the analysis. Energy demands for ships, railways, and aircraft in the transport sector are also excluded from the study because it is difficult to distinguish between specific regions.

Estimation of Energy Supply from Renewable Energy
As the upper limit on renewable generation technology has installed capacity in the optimization model, the introduction potential is referenced [40]. The hourly output standardized by the installed capacity from PV and wind turbine is estimated from the sample year's weather data. First, the standardized output of PV is calculated from Equation (55). The amount of solar radiation is referenced from meteorological data at an observation point in each prefecture [41]. Other technical parameters of PV are referenced from a product catalog [42].
The standardized output of wind power generators is estimated by the same method as Obara et al. [14]. The wind speed at the ground level at the observation point is converted to the wind speed at the hub height, and this wind speed is substituted into the wind turbine's output curve. Cut-in wind speed is 2.5 m/s, rated wind speed is 12.5 m/s, and cut-out wind speed is 25.0 m/s. The wind speed at all observation points in each prefecture [41] is referenced for estimating the output of onshore wind power generators. The wind speed at the cape in each prefecture [43] is referenced to calculate offshore wind power generation's output.
The standardized output from geothermal is obtained by dividing the power company's actual data by the current geothermal power generation capacity. As for hydropower generation, it is assumed that the existing power plant is maintained in light of the situation where no further expansion is expected in Japan. The generation from run-of-river is considered to have a constant output at a load factor of 60%. The biomass resource optimized by the model is limited to woody biomass. It is assumed that the woody biomass resource that can be supplied per year is the annual stand growth volume [44]. Table 2 summarizes the model's technologies and its technical and economic parameters. The annual unit cost of technology t is calculated (Equations (56)-(58)), considering the overnight cost, operation and maintenance cost, lifetime, and efficiency, in Table 2. The interest rate I is assumed to be uniformly 2%. In this study, the target synthetic fuels are methane, methanol, and DME. This is because the synthesis of gasoline and diesel oil is costly and has a lower conversion efficiency [45]; therefore, there is no opportunity to introduce it by optimization. In this model, the introduction cost of energy conversion Energies 2021, 14, 2823 12 of 26 equipment in the customer is not considered; therefore, the introduction cost of local heat pump, heat exchanger, FCV, and BEV is ignored. Similarly, the introduction cost of filling facilities that supply energy to FCVs and BEVs is not considered. The efficiencies of solar PV, onshore wind, offshore wind, run-of-river, and geothermal are 1.00 because those resources are not depleting and utilize some of the almost inexhaustible energy. The efficiencies of biomass-CHP and gas-CHP in Table 2 are electrical efficiency. The efficiency of the fuel cell in Table 2 is also electrical efficiency. The thermal efficiency of HOB exceeds 1 because it is calculated using the lower calorific value.

Technical and Economic Parameters
It assumes an air heat source heat pump as a local heat pump in commercial and residential sectors. The local heat pump's coefficient of performance (COP) is affected by the ambient temperature and the supply temperature [54]. In this analysis, the supply temperature is assumed to be 60 • C. The heat exchange efficiency in buildings connected to the DHS is 80%. The value of output ratio of stationary battery is 0.22 [kw/kWh]. It was set from the cost per output [Euro/kW] and the cost per capacity [Euro/kWh]. The vehicle's battery capacity and electricity output for V2G is assumed to be 62 kWh for battery capacity and 6 kW for output [55]. SOC of batteries in vehicles that participate in V2G is between 30% and 90% to maintain the batteries' life cycle within an acceptable range [56].
It is assumed that the fuel cost of woody biomass is 12,750 JPY/t, and the unit calorific value is 15 MJ/kg [46]. The energy consumption for transporting woody biomass is not taken into account. The cost of biomass resources supplied to the construction industry, pulp, and paper industry, wood product manufacturing industry, etc., within current statistics is not evaluated. The heat recovery efficiency from the electrolyzer and fuel synthesis is assumed to be 90%. It is assumed that CO 2 captured by DAC requires 225 kWh/t-CO 2 of electricity and 1500 kWh/t-CO 2 of high-temperature heat [53]. CO 2 capture cost by DAC is assumed to be 2594 JPY/t-CO 2 [53]. The steam turbine is assumed to be the prime mover of the biomass CHP, and the gas turbine is assumed to be the prime mover of the gas CHP. Assuming that the energy required for CO 2 capture from the CHP plant is produced in the plant, it can be expressed by a decrease in energy conversion efficiency. According to [48], it was shown that the power generation efficiency of coal-fired power plants that capture CO 2 is reduced by 23%, compared to the case without capture.
Therefore, it is assumed that the efficiency of the biomass plant is also reduced by 23% from the efficiency without CO 2 capture, which is 46%. For gas CHP plants, efficiency is reported to decrease by 14%. The CO 2 capture efficiency is assumed to be 90%. The CO 2 emission factors for biomass and methane are considered to be 98.7 kg-CO 2 /MJ [57] and 51.1 kg-CO 2 /MJ [58], respectively. The capture costs are also assumed to be 5060 JPY/t-CO 2 and 8140 JPY/t-CO 2 , respectively [48]. It is assumed that synthetic fuels' storage facilities are already in use. In this case, the capital cost of the storage facility is ignored. The storage capacity of methane, methanol, and DME is referenced from [59,60]. Table 3 summarizes the scenarios in this study. There are three scenarios on transportation demands. As with the report by U. Persson et al. [61], heat demand is classified into areas of 120 TJ/km 2 or less, 120-300 TJ/km 2 (dense heat density class), and 300 TJ/km 2 or more (very thick heat density class) according to the heat-demand density. Individual or district heating is provided to each of these areas. There are three scenarios on heat demands. Table 3. Description of scenarios analyzed in this study.

Transportation
All-FCV All transportation is provided by fuel cell vehicles. All-BEV All transportation is provided by battery electric vehicles.
All-BEV-V2G All transportation is provided by battery electric vehicles, and the normal-sized vehicle can be used for vehicle-to-grid.
Heat demand IND All heat demand is met by individual heating. HD300 DHS is introduced in areas where the heat-demand density is over 300 TJ/km 2 HD120 DHS is introduced in areas where the heat-demand density is over 120 TJ/km 2 Figure 3a shows the configuration of the primary energy supply of each scenario based on the results of the optimization calculation. TPES is estimated by Equation (59). TPES in the current energy system is estimated using statistics [28], not using the proposed model. Compared with the current, TPES decreases in almost all scenarios because of the reduction of heat loss at the existing thermal power plants, the improvements of tank-towheel efficiency by a shift to BEV and FCV, and the improvements of energy efficiency to produce heat in commercial and residential sectors. TPES of All-FCV scenarios are larger than that of All-BEV and All-BEV+V2G scenarios because the tank-to-wheel efficiency of FCV is lower. There is no significant difference in TPES between All-BEV and All-BEV-V2G scenarios. TPES of HD300 and HD120 scenarios tend to increase slightly compared to IND scenarios because of the difference in efficiency between the heat exchanger and local HP in consumers. Offshore wind accounts for more than 70% of the total primary energy supply (r, h) + c HOB bio.→lth (r, h) (59) Figure 3b shows the breakdown of the total annual system cost, the optimization model's objective function. Regarding the total annual system cost in the current energy system, the generation cost is estimated from the power generation amount and generation unit price [JPY/kWh] (including capital cost and O&M cost) of each technology, and the fuel cost is estimated from the primary energy supply amount and wholesale unit price of each fuel, not from the proposed model [46,58]. As a result, the current energy system's total annual system cost was 4079.7 billion JPY/year. The total annual system cost increases from 80% to 140% of the current system cost. All-BEV scenarios' annual total system cost was lower than that of all-FCV scenarios because all-FCV scenarios' total primary energy supply is larger. The whole system cost of All-BEV-V2G scenarios is reduced by about 10% compared to all-BEV scenarios. The capital and O&M costs of CHP, electrolyzer, and stationary storage battery are mainly reduced by introducing V2G. Comparing IND scenarios with HD300 and HD120 scenarios, the DHS's introduction reduced the total system cost. The distribution cost in DHS and capital and O&M costs of heat source equipment and heat storage tanks are less than 1% of the entire system cost. Compared to IND scenarios, the cost of electricity generation tends to increase. However, the total system costs are decreased since the cost reduction of stationary storage batteries and electrolyzers exceeds the increase. the fuel cost is estimated from the primary energy supply amount and wholesale unit price of each fuel, not from the proposed model [46,58]. As a result, the current energy system's total annual system cost was 4079.7 billion JPY/year. The total annual system cost increases from 80% to 140% of the current system cost. All-BEV scenarios' annual tota system cost was lower than that of all-FCV scenarios because all-FCV scenarios' total primary energy supply is larger. The whole system cost of All-BEV-V2G scenarios is reduced by about 10% compared to all-BEV scenarios. The capital and O&M costs of CHP, electrolyzer, and stationary storage battery are mainly reduced by introducing V2G. Comparing IND scenarios with HD300 and HD120 scenarios, the DHS's introduction reduced the total system cost. The distribution cost in DHS and capital and O&M costs of heat source equipment and heat storage tanks are less than 1% of the entire system cost. Compared to IND scenarios, the cost of electricity generation tends to increase. However, the total system costs are decreased since the cost reduction of stationary storage batteries and electrolyzers exceeds the increase.

Installed Capacity of Each Technology
This section shows the solved installed capacity, a variable of the total annual system cost. All-BEV_IND, All-BEV_HD120, All-BEV-V2G_IND, and All-BEV-V2G_HD120 scenarios are compared. The installed capacity of each technology is shown in Figure 4. The installed capacity of offshore wind generation in Aomori, Akita, and Niigata prefectures tends to be large. Miyagi and Yamagata prefectures have a large percentage of installed solar PV and onshore wind generation capacity. No generation technology is installed in Iwate. This result is explained by the relationship between the levelized cost of electricity (LCOE) of each generation technology and the transmission costs between regions. A generation technology in a specific region is installed at a lower price than the LCOE in the region and the LCOE in the other region, considering the transmission cost. For example the offshore wind LCOE in Aomori Prefecture is 9.9 JPY/kWh. If this electricity is consumed in Niigata Prefecture, its actual cost would be 24.4 JPY/kWh. In other words, LCOE for different generation technologies in Niigata is cheaper than offshore wind power in Aomori.

Installed Capacity of Each Technology
This section shows the solved installed capacity, a variable of the total annual system cost. All-BEV_IND, All-BEV_HD120, All-BEV-V2G_IND, and All-BEV-V2G_HD120 scenarios are compared. The installed capacity of each technology is shown in Figure 4. The installed capacity of offshore wind generation in Aomori, Akita, and Niigata prefectures tends to be large. Miyagi and Yamagata prefectures have a large percentage of installed solar PV and onshore wind generation capacity. No generation technology is installed in Iwate. This result is explained by the relationship between the levelized cost of electricity (LCOE) of each generation technology and the transmission costs between regions. A generation technology in a specific region is installed at a lower price than the LCOE in the region and the LCOE in the other region, considering the transmission cost. For example, the offshore wind LCOE in Aomori Prefecture is 9.9 JPY/kWh. If this electricity is consumed in Niigata Prefecture, its actual cost would be 24.4 JPY/kWh. In other words, LCOE for different generation technologies in Niigata is cheaper than offshore wind power in Aomori.  The characteristics of the installed capacity of energy conversion technology differ depending on the DHS's utilization. In the All-BEV_IND scenario, many liquefaction facilities are installed in Aomori Prefecture. Since many offshore wind power generation plants are installed in Aomori Prefecture, much flexibility for intermittent power generation is required. The results suggest that P2G is the best option to resolve renewable energy's intermittency. Electrolyzers and fuel synthesizers account for the most installed capacity, particularly in Miyagi, Niigata, and Fukushima, where industrial fuel demand is high. In All-BEV_HD120 and All-BEV-V2G_HD120 scenarios, many EBs are installed in Aomori. This is because heat storage by EBs and HWTs is the cheapest option as flexibility for offshore wind power in Aomori. The electrolyzer and fuel synthesizer's installed capacity tends to be more prominent in the area where energy consumption in the industrial sector is high. In these scenarios, it can be confirmed that the installed capacity of the electrolyzer, fuel synthesizer, stationary battery, and HPT is reduced as compared with the All-BEV_IND scenario. Comparing the All-BEV_HD120 scenario with the All-BEV-V2G_HD120 scenario, the installed capacity of EB, CHP, electrolyzer, and fuel synthesizer are all reduced by the introduction of V2G.
The above results suggest that the use of DHS and V2G can reduce the installed capacity of P2G technologies such as electrolyzer and fuel synthesizers. This is because the flexibility offered by P2H and V2G for electricity networks makes the operations of P2G technologies effective. Figure 5 shows the duration curves of (a) electrolyzer and (b) methane synthesizer. The vertical axis intercept corresponds with the technology's installed capacity in these figures. When comparing All-BEV and All-BEV-V2G scenarios, the installed capacity of P2G technologies in All-BEV-V2G scenarios is reduced. The demand for hydrogen is almost equal between the All-BEV scenario and the All-BEV-V2G scenario, but the slope of the duration curve has become gentler with the introduction of V2G. This means that the electrolyzer's capacity factor (CF) has been improved, and the installed capacity has been reduced. When comparing IND scenarios with HD120 scenarios, it can be confirmed that the installed capacity of P2G technologies in HD120 scenarios has been reduced. The electrolyzer's capacity factor has improved from 0.42 in the All-BEV_IND scenario to 0.55 in the All-BEV_HD120 scenario. The characteristics of the installed capacity of energy conversion technology differ depending on the DHS's utilization. In the All-BEV_IND scenario, many liquefaction facilities are installed in Aomori Prefecture. Since many offshore wind power generation plants are installed in Aomori Prefecture, much flexibility for intermittent power generation is required. The results suggest that P2G is the best option to resolve renewable energy's intermittency. Electrolyzers and fuel synthesizers account for the most installed capacity, particularly in Miyagi, Niigata, and Fukushima, where industrial fuel demand is high. In All-BEV_HD120 and All-BEV-V2G_HD120 scenarios, many EBs are installed in Aomori. This is because heat storage by EBs and HWTs is the cheapest option as flexibility for offshore wind power in Aomori. The electrolyzer and fuel synthesizer's installed capacity tends to be more prominent in the area where energy consumption in the industrial sector is high. In these scenarios, it can be confirmed that the installed capacity of the electrolyzer, fuel synthesizer, stationary battery, and HPT is reduced as compared with the All-BEV_IND scenario. Comparing the All-BEV_HD120 scenario with the All-BEV-V2G_HD120 scenario, the installed capacity of EB, CHP, electrolyzer, and fuel synthesizer are all reduced by the introduction of V2G.
The above results suggest that the use of DHS and V2G can reduce the installed capacity of P2G technologies such as electrolyzer and fuel synthesizers. This is because the flexibility offered by P2H and V2G for electricity networks makes the operations of P2G technologies effective. Figure 5 shows the duration curves of (a) electrolyzer and (b) methane synthesizer. The vertical axis intercept corresponds with the technology's installed capacity in these figures. When comparing All-BEV and All-BEV-V2G scenarios, the installed capacity of P2G technologies in All-BEV-V2G scenarios is reduced. The demand for hydrogen is almost equal between the All-BEV scenario and the All-BEV-V2G scenario, but the slope of the duration curve has become gentler with the introduction of V2G. This means that the electrolyzer's capacity factor (CF) has been improved, and the installed capacity has been reduced. When comparing IND scenarios with HD120 scenarios, it can be confirmed that the installed capacity of P2G technologies in HD120 scenarios has been reduced. The electrolyzer's capacity factor has improved from 0.42 in the All-BEV_IND scenario to 0.55 in the All-BEV_HD120 scenario. Furthermore, it has improved from 0.50 in the All-BEV-V2G_IND scenario to 0.63 in the All-BEV-V2G_HD120 scenario. As in the electrolyzer, the slope of the duration curves Furthermore, it has improved from 0.50 in the All-BEV-V2G_IND scenario to 0.63 in the All-BEV-V2G_HD120 scenario. As in the electrolyzer, the slope of the duration curves of the methane synthesizer becomes gentle due to the introduction of DHS and V2G.
In other words, the capacity factor has been improved, and DHS and V2G have reduced the installed capacity. Figure 6 shows the hourly electricity balances in each scenario. The electrolyzer and methane synthesizer's hourly operations are also shown in Figure 7. In the electricity balances of the All-BEV_IND scenario, fluctuations due to renewable energy can be adjusted by a stationary battery and the electrolyzer. As a result, the electrolyzer operation is selected as the primary option to offset the fluctuations. Therefore, the electrolyzer operation greatly depends on the fluctuation of renewable energy and is not stable, as shown in Figure 7a. Most of the produced hydrogen is used for methane synthesis, and the output of the synthesizer also fluctuates significantly, as shown in Figure 7b. In the electricity balances of the All-BEV_HD120 scenario, renewable energy fluctuations are adjusted by operating the stationary battery, electrolyzer, and electric boiler. Electric boilers operate at the peak of renewable energy supply. The operation of electrolyzer and methane synthesizer depends on renewable energy fluctuation. However, the outputs are partially suppressed by the peak cut effect of P2H, as shown in Figure 7a,b.

Hourly Energy Balances in the Sector-Coupled Energy System
In the electricity balances of the All-BEV-V2G_IND scenario, fluctuations are adjusted by the operating electrolyzer and the use of V2G. The peak shift effect by V2G is working most of the time. Therefore, the electrolyzer and methane synthesizer operation in this scenario is more stable than that of the All-BEV_IND scenario, as shown in Figure 7a,b. In the electricity balance of the All-BEV-V2G_HD120 scenario, fluctuations of renewable energy are adjusted by operating the electrolyzer, the electric boiler, and V2G. The peak shift effect by V2G is working most of the time. Electric boilers are partially in operation at the peak of renewable energy supplies. The electrolyzer operation is leveled by the peak shift effect of V2G and the peak cut effect of P2H. The electrolyzer operation keeps an almost constant output, showing the peak cut and peak shift effects of P2H and V2G. Therefore, the output of the fuel synthesizer is also stable. of the methane synthesizer becomes gentle due to the introduction of DHS and V2G. In other words, the capacity factor has been improved, and DHS and V2G have reduced the installed capacity. Figure 6 shows the hourly electricity balances in each scenario. The electrolyzer and methane synthesizer's hourly operations are also shown in Figure 7. In the electricity balances of the All-BEV_IND scenario, fluctuations due to renewable energy can be adjusted by a stationary battery and the electrolyzer. As a result, the electrolyzer operation is selected as the primary option to offset the fluctuations. Therefore, the electrolyzer operation greatly depends on the fluctuation of renewable energy and is not stable, as shown in Figure 7a. Most of the produced hydrogen is used for methane synthesis, and the output of the synthesizer also fluctuates significantly, as shown in Figure 7b. In the electricity balances of the All-BEV_HD120 scenario, renewable energy fluctuations are adjusted by operating the stationary battery, electrolyzer, and electric boiler. Electric boilers operate at the peak of renewable energy supply. The operation of electrolyzer and methane synthesizer depends on renewable energy fluctuation. However, the outputs are partially suppressed by the peak cut effect of P2H, as shown in Figure 7a,b.

Hourly Energy Balances in the Sector-Coupled Energy System
In the electricity balances of the All-BEV-V2G_IND scenario, fluctuations are adjusted by the operating electrolyzer and the use of V2G. The peak shift effect by V2G is working most of the time. Therefore, the electrolyzer and methane synthesizer operation in this scenario is more stable than that of the All-BEV_IND scenario, as shown in Figure  7a,b. In the electricity balance of the All-BEV-V2G_HD120 scenario, fluctuations of renewable energy are adjusted by operating the electrolyzer, the electric boiler, and V2G. The peak shift effect by V2G is working most of the time. Electric boilers are partially in operation at the peak of renewable energy supplies. The electrolyzer operation is leveled by the peak shift effect of V2G and the peak cut effect of P2H. The electrolyzer operation keeps an almost constant output, showing the peak cut and peak shift effects of P2H and V2G. Therefore, the output of the fuel synthesizer is also stable. Feb. 13  Previous discussions suggest that the introduction of DHS and V2G could reduce the installed capacity of electrolyzer and fuel synthesizers, contributing to lower total annual system costs. When this is analyzed in detail, it becomes clear that the peak cut effect by P2H and the peak shift effect by V2G result in the leveling of the output of the electrolyzer and fuel synthesizer, which improves the capacity factor and reduces the installed capacity. In other words, this result suggests that the introduction of DHS and V2G not only enhances the flexibility of the electrical grid but also reduces the production cost of synthetic fuels for industry and transportation.

Comparison with the Current Energy System
The sector-coupled energy system's characteristics are explained by comparing the current energy system with the energy flow diagrams. Figure 8 shows the energy flow diagram of (a) the current energy system and (b) the energy system considering the sector  Hourly electricity input to electrolyzer Feb. 13 -19 All-BEV_IND All-BEV_HD120 All-BEV-V2G_IND All-BEV-V2G_HD120 Feb. 13 -19 All-BEV_IND All-BEV_HD120 All-BEV-V2G_IND All-BEV-V2G_HD120 Previous discussions suggest that the introduction of DHS and V2G could reduce the installed capacity of electrolyzer and fuel synthesizers, contributing to lower total annual system costs. When this is analyzed in detail, it becomes clear that the peak cut effect by P2H and the peak shift effect by V2G result in the leveling of the output of the electrolyzer and fuel synthesizer, which improves the capacity factor and reduces the installed capacity. In other words, this result suggests that the introduction of DHS and V2G not only enhances the flexibility of the electrical grid but also reduces the production cost of synthetic fuels for industry and transportation. Previous discussions suggest that the introduction of DHS and V2G could reduce the installed capacity of electrolyzer and fuel synthesizers, contributing to lower total annual system costs. When this is analyzed in detail, it becomes clear that the peak cut effect by P2H and the peak shift effect by V2G result in the leveling of the output of the electrolyzer and fuel synthesizer, which improves the capacity factor and reduces the installed capacity. In other words, this result suggests that the introduction of DHS and V2G not only enhances the flexibility of the electrical grid but also reduces the production cost of synthetic fuels for industry and transportation.

Comparison with the Current Energy System
The sector-coupled energy system's characteristics are explained by comparing the current energy system with the energy flow diagrams. Figure 8 shows the energy flow diagram of (a) the current energy system and (b) the energy system considering the sector  Feb. 13 -19 All-BEV_IND All-BEV_HD120 All-BEV-V2G_IND All-BEV-V2G_HD120 Feb. 13 -19 All-BEV_IND All-BEV_HD120 All-BEV-V2G_IND All-BEV-V2G_HD120

Comparison with the Current Energy System
The sector-coupled energy system's characteristics are explained by comparing the current energy system with the energy flow diagrams. Figure 8 shows the energy flow diagram of (a) the current energy system and (b) the energy system considering the sector coupling designed by optimization. The system shown in (b) is based on a scenario in which all vehicles in the transport sector are electric vehicles, of which passenger cars are used for V2G, and a district heating system is introduced into areas where the heat-demand density is 120 TJ/km 2 or more. The current energy flow diagram is based on statistical data estimates [27,28,30].
In the current energy system in the Tohoku region, TPES is 1495 PJ/year. The breakdown is 11.4 PJ for solar, 6.1 PJ for onshore wind, 37.5 PJ for hydro, 2.3 PJ for geothermal, 76.2 PJ for biomass, 343.5 PJ for coal, 335.1 PJ for natural gas, and 682.7 PJ for crude oil. Therefore, fossil fuels make up 91% of the total. The primary energy input to electricity generation is 589.9 PJ, with renewable energy accounting for 12.4% and fossil fuels accounting for 87.6%. In addition, 293.9 PJ, which accounts for 50% of the primary input energy, is lost at thermal power plants. The final energy consumption is 1270.7 PJ (electricity, 23.3%; biomass, 4.7%; and fossil fuels, 71.0%). The energy system's rejected energy is 293.9 PJ, and the system energy efficiency is 80.3%. The Tohoku region's current amount of CO 2 emissions was estimated to be 95.0 Mt-CO 2 /year.
In the All-BEV-V2G_HD120 scenario, which is an example of a sector-coupled energy system, TPES is 1323.5 PJ/year. This is approximately 11.5% less than the current situation. The breakdown is 72.4 PJ for solar, 148.0 PJ for onshore wind, 938.1 PJ for offshore wind, 46.6 PJ for hydro, and 118.8 PJ for biomass. In the electricity generation sector, 1297.0 PJ of primary energy is supplied, generating 1237.9 PJ, approximately 4.2 times the current electricity generation. The electricity generation sector's rejected energy is 29.3 PJ (about 10% of the current level).
Final energy consumption decreases by 17.2% from the current level to 1052.7 PJ. This is due to the shift from internal combustion engine vehicles to electric vehicles and heatdemand electrification. The energy system's rejected energy is 300.7 PJ, and the system's energy efficiency is 77%. This is due to the energy loss during hydrogen production for fuel synthesis. Approximately 64% of the total electricity generation is consumed by the electrolyzer's operation, producing 489.2 PJ of hydrogen. The conversion loss is 224.1 PJ, and the waste heat recovery in the district heating system is 75.7 PJ. The energy demand of the CO 2 capture process by DAC is 140.4 PJ, accounting for about 12% of the total final energy consumption. In this study, it is assumed that synthetic fuels cover the fuel demand of the industrial sector, meaning a large amount of energy is lost when producing hydrogen for fuel synthesis. To improve energy system efficiency, reduction and electrification of fuel demand in the industrial sector are necessary. coupling designed by optimization. The system shown in (b) is based on a scenario in which all vehicles in the transport sector are electric vehicles, of which passenger cars are used for V2G, and a district heating system is introduced into areas where the heat-demand density is 120 TJ/km 2 or more. The current energy flow diagram is based on statistical data estimates [27,28,30].
In the current energy system in the Tohoku region, TPES is 1495 PJ/year. The breakdown is 11.4 PJ for solar, 6.1 PJ for onshore wind, 37.5 PJ for hydro, 2.3 PJ for geothermal, 76.2 PJ for biomass, 343.5 PJ for coal, 335.1 PJ for natural gas, and 682.7 PJ for crude oil . Therefore, fossil fuels make up 91% of the total. The primary energy input to electricity generation is 589.9 PJ, with renewable energy accounting for 12.4% and fossil fuels accounting for 87.6%. In addition, 293.9 PJ, which accounts for 50% of the primary input energy, is lost at thermal power plants. The final energy consumption is 1270.7 PJ (electricity, 23.3%; biomass, 4.7%; and fossil fuels, 71.0%). The energy system's rejected energy is 293.9 PJ, and the system energy efficiency is 80.3%. The Tohoku region's current amount of CO2 emissions was estimated to be 95.0 Mt-CO2/year.
In the All-BEV-V2G_HD120 scenario, which is an example of a sector-coupled energy system, TPES is 1323.5 PJ/year. This is approximately 11.5% less than the current situation. The breakdown is 72.4 PJ for solar, 148.0 PJ for onshore wind, 938.1 PJ for offshore wind, 46.6 PJ for hydro, and 118.8 PJ for biomass. In the electricity generation sector, 1297.0 PJ of primary energy is supplied, generating 1237.9 PJ, approximately 4.2 times the current electricity generation. The electricity generation sector's rejected energy is 29.3 PJ (about 10% of the current level).
Final energy consumption decreases by 17.2% from the current level to 1052.7 PJ. This is due to the shift from internal combustion engine vehicles to electric vehicles and heatdemand electrification. The energy system's rejected energy is 300.7 PJ, and the system's energy efficiency is 77%. This is due to the energy loss during hydrogen production for fuel synthesis. Approximately 64% of the total electricity generation is consumed by the electrolyzer's operation, producing 489.2 PJ of hydrogen. The conversion loss is 224.1 PJ, and the waste heat recovery in the district heating system is 75.7 PJ. The energy demand of the CO2 capture process by DAC is 140.4 PJ, accounting for about 12% of the total final energy consumption. In this study, it is assumed that synthetic fuels cover the fuel demand of the industrial sector, meaning a large amount of energy is lost when producing hydrogen for fuel synthesis. To improve energy system efficiency, reduction and electrification of fuel demand in the industrial sector are necessary. (a)

Toward Implementation of the Sector-Coupled Energy System in Japan
In a sector-coupled energy system, the total primary energy supply decreases, but the amount of generation may be required to be about four times the current amount. In the Tohoku region, the share of renewable resources in the current generation is 12.4%. In order to realize this system, it is necessary to use 34 times as many renewable resources as is currently introduced. From the analysis of energy flow, it is confirmed that in the optimization model proposed in this study, the synthetic fuels supplied to the industrial sector have a considerable influence on the amount of power generation. The higher the share of energy consumption of the industrial sector in final energy consumption, the more power generation required. In the Tohoku region, the energy consumption of the industrial sector accounts for 35.4% of total final energy consumption, but in Japan, the energy consumption of the industrial sector accounts for 46.1% of total final energy consumption. In addition, the share of renewable resources in the current generation is 16% (2017). If this model is applied to Japan, simple calculations would not be possible, but it can be expected that it would be necessary to introduce more than 35 times the current amount of renewable energy. To verify the effects of sector coupling in Japan as a whole, an energy system model covering all regions is needed.
On the other hand, this study does not assume changes in energy demand. Energy demand is greatly affected by the economy and vital statistics. The Japanese population is already declining, and energy demand is also expected to decline. In addition, countermeasures against COVID-19 in 2020 have significantly changed our lifestyles, significantly impacting energy demand. If energy demand decreases, less renewable energy would be required; therefore, efforts to transform lifestyles on the demand side into sustainable ones are also important.
In addition, in sector-coupled energy systems, the total system cost may be more than double the current cost. To reduce the cost as much as possible, it is appropriate to promote electric vehicles in the transportation sector because fuel cell vehicles are inferior to electric vehicles in terms of overall energy efficiency. In addition, there is a possibility that

Toward Implementation of the Sector-Coupled Energy System in Japan
In a sector-coupled energy system, the total primary energy supply decreases, but the amount of generation may be required to be about four times the current amount. In the Tohoku region, the share of renewable resources in the current generation is 12.4%. In order to realize this system, it is necessary to use 34 times as many renewable resources as is currently introduced. From the analysis of energy flow, it is confirmed that in the optimization model proposed in this study, the synthetic fuels supplied to the industrial sector have a considerable influence on the amount of power generation. The higher the share of energy consumption of the industrial sector in final energy consumption, the more power generation required. In the Tohoku region, the energy consumption of the industrial sector accounts for 35.4% of total final energy consumption, but in Japan, the energy consumption of the industrial sector accounts for 46.1% of total final energy consumption. In addition, the share of renewable resources in the current generation is 16% (2017). If this model is applied to Japan, simple calculations would not be possible, but it can be expected that it would be necessary to introduce more than 35 times the current amount of renewable energy. To verify the effects of sector coupling in Japan as a whole, an energy system model covering all regions is needed.
On the other hand, this study does not assume changes in energy demand. Energy demand is greatly affected by the economy and vital statistics. The Japanese population is already declining, and energy demand is also expected to decline. In addition, countermeasures against COVID-19 in 2020 have significantly changed our lifestyles, significantly impacting energy demand. If energy demand decreases, less renewable energy would be required; therefore, efforts to transform lifestyles on the demand side into sustainable ones are also important.
In addition, in sector-coupled energy systems, the total system cost may be more than double the current cost. To reduce the cost as much as possible, it is appropriate to promote electric vehicles in the transportation sector because fuel cell vehicles are inferior to electric vehicles in terms of overall energy efficiency. In addition, there is a possibility that the cost can be further reduced by utilizing the electric vehicle as a flexible carrier by V2G. It was confirmed that the district heating system also reduces the total system cost due to its flexibility through P2H. Japan is in a situation where it is difficult to secure the flexibility to the grid based on the "cross-border" concept that accommodates electricity between countries. Therefore, securing flexibility based on the "cross-sector" concept is the key to the mass introduction of renewable resources into the energy system. However, the penetration rate of district heating systems in Japan is low, and it is believed that they will not spread significantly in the future. Therefore, policy support is needed to promote district heating systems toward decarbonized society.
From the results of this study, the following suggestions were obtained regarding the roles required for each technology: First, regarding power generation technology, when there are few options for flexibility, it tends to be difficult to introduce PV generation. This is due to the characteristic that PV generation peaks during the day. The mass introduction of PV generation increases storage batteries and conversion technologies to cope with this peak, and the cost may increase accordingly. The results show that it is desirable to introduce a large amount of offshore wind power in areas with low LCOE in the target areas of this study. In the energy system where all couplings were taken into consideration, electric boiler, electric heat pump, and hot water tank in the district heating system have the effect of improving the capacity factor of electrolyzer in addition to the function as the flexibility to the electrical grid. Therefore, it was suggested that it might reduce the production cost of synthetic fuel. V2G also showed behavior that improves the capacity factor of the electrolyzer. In this study, P2G became a bottleneck for cost reduction. Therefore, when synthetic fuels are needed in areas where the share of industrial energy is high, it is better to avoid using P2G as the flexibility to the electrical grid. It is necessary to improve the capacity factor by making the output of the electrolyzer as stable as possible.

Possible Improvements
To consider reduction and electrification of fuel demands and direct use of hydrogen in the industrial sector, it is necessary to update the optimization model that includes the consumer's replacement cost of equipment. In the industrial sector, the chemical industry uses naphtha derived from crude oil distillation. Therefore, it is necessary to consider raw materials for chemical products not derived from crude oil. For example, there is synthesis from electrolyzed hydrogen and CO 2 . There is a need for a model that considers reducing energy consumption due to improved energy efficiency through insulation and equipment replacement in the commercial and residential sectors. It is also desirable to consider waste heat from waste incineration plants in the district heating system. In the transportation sector, it is necessary to develop models that consider the recycling of batteries in BEV, the deterioration of batteries due to the use of V2G, the cost of power supply facilities to BEV, and the energy consumption of railways, ships, and aircraft. Moreover, it is necessary to evaluate the impact on V2G usage due to the decrease in the number of vehicles by car-sharing and ridesharing. In addition, the use of pumped storage power generation should be considered as flexible for variable renewable energy.
Moreover, in this study, it was not possible to analyze the technology that utilizes solar thermal. There was a lack of data on solar thermal in Japan, and reliable data could not be obtained; therefore, the model could not handle it. In the future, it will be necessary to prepare reliable data on solar heat and incorporate it into the optimization model. As with solar thermal, there is a lack of data on a load of refrigeration and cooling, meaning it is necessary to first prepare basic data and incorporate it into the optimization model. Furthermore, there is a need to develop a model that considers cost reductions due to scale factors. Although this leads to a non-linear model, it is possible to estimate the introduction costs of electrolyzer and fuel synthesizer more precisely. In this study, we assumed that only electric power is the cross-border target. However, there is a possibility that a more efficient system can be designed by considering the cross-border of synthetic fuels.
In addition, in this model, it is a prerequisite that the introduction potential of renewable energy greatly exceeds the demand, which gives a feasible solution. For example, in Japan, it would be difficult to apply this model in areas where energy demand exceeds the introduction potential of renewable energy, such as around Tokyo and Osaka. Therefore, it is necessary to make it possible for the model to handle centralized power sources such as hydrogen power generation, ammonia power generation, and nuclear power generation. It is also necessary to extend this model to all regions of Japan by considering inter-regional transmission, reinforcement of transmission lines, and fuel transportation.

Conclusions
In this study, an optimization model of the energy system considering sector coupling has been developed and applied to the Tohoku region to reveal the required technologies and their costs and optimal operation in the sector-coupled energy system. As a result, the following findings are obtained.

•
The total primary energy supply will decrease, but the amount of generation may be required to be about four times the current amount. In Japan, it can be expected that it would be necessary to introduce more than 35 times the current amount of renewable energy to realize decarbonization by renewable resources. The annual total system cost of the designed energy systems increases from 1.8 to 2.4 times the current level. • The annual system cost of the scenario introducing district heating systems and using electric vehicles for V2G is minimal. As a result of the analysis of hourly energy balance, it becomes clear that the peak cut effect by P2H and the peak shift effect by V2G result in the leveling of the electrolyzer, and fuel synthesizer's output improves the capacity factor and reduces the introduction capacity. It is necessary to utilize electric vehicles by V2G and provide policy support for district heating systems to reduce costs. • From the analysis of energy flow, it is confirmed that in the optimization model proposed in this study, the synthetic fuels supplied to the industrial sector have a considerable influence on the amount of power generation.

•
There is a concern that the capacity of storage batteries and energy conversion technologies will increase to cope with the peak of PV. In areas where the industrial energy share is high, P2H and V2G may help reduce the cost of synthetic fuels. In addition, P2G should be mainly aimed at producing synthetic fuels, not flexibility.

•
To reduce the amount of renewable energy required, it is necessary to reduce the energy demand in the industrial sector and electrify the energy demand. To reduce energy demand in the industrial sector, it is important to transform lifestyles that were not considered in this study into sustainable ones. • For future research, it is important to incorporate more technologies, improve the model's accuracy with reliable data, and model the optimization on the demand side. It is also required to apply the model to all regions of Japan.  Total costs for methane liquefaction (million JPY/year) C total Total annual system costs (million JPY/year) C trans.
Total transmission costs (million JPY/year) CC bio. , CC gas , CC DAC The total amount of captured CO 2 (t-CO 2 /year) pes bioCHP bio.
Primary energy supply to CHP plants by woody biomass (GW) pes Woody biomass stock (GWh) Sto CO 2 The amount of CO 2 stored in tanks (t-CO 2 ) V2GCap Total allowed maximum vehicle-to-grid outputs (GW) Parameters A Photovoltaic panel area (m 2 /kW) c overnight t Overnight capital cost (million JPY/GW) Cap trs.
Transmission line capacity (GW) cap vbattery The storage capacity of a battery in a vehicle (GWh/vehicle) cd n CO 2 demands for fuels synthesis (t-CO 2 /GW hydrogen ) ce f bio. , ce f CH 4 The carbon emission factor of woody biomass and methane (t-CO 2 /GWh) CRF t Capital recovery factor (-) dmd manu f . Annual fuels demand in the manufacturing industry (GWh) dmd nmanu f . Annual fuels demand in the non-manufacturing industry (GWh) dmd transport. elect.
Electricity consumption by BEV driving (GW) dmd elect.