Rapport Final - Groupe 2
Rapport Final - Groupe 2
Introduction
The objective of our industrial project is to establish a viable process for the
preparation of methanol, using green hydrogen obtained by electrolysis of water and
capturing collected carbon dioxide as feedstock, and producing methanol by a catalytic
reaction.
At the start of the project, we consulted a number of literature to determine
information including the appropriate conditions for the chemical reaction to produce
methanol, the appropriate choice of catalyst, and the amount of feedstock required.
Based on the results of this preliminary investigation, we carried out simulations using
PROII software to model the methanol preparation process in detail and to simulate
parameters such as flow rate, temperature, and pressure for each unit in the process.
From the simulation results, we then were required to carry out a preliminary estimate
of the equipment size. Using the equations and assumptions provided, we calculated the
dimensions of the reactor, pumps, compressors, and other equipment. However, we also
realized that these results were only to be used as a reference and that the actual design
would require further refinement and cooperation with the equipment supplier.
Finally, we carried out a CAPEX and OPEX economic analysis based on the previously
obtained simulation and sizing data to assess the economic viability of the project.
Methods
Preliminary research:
To carry out the project, we first carried out a literature search and had an initial insight
into the methanol preparation process. This revealed that favorable operating conditions for
the reaction are a temperature of 200-300°C and a pressure of 50-100 bar. A heterogeneous
catalyst of the Cu/ZnO/Al2O3 type is also required. To achieve the reaction conditions,
equipment such as heat exchangers, compressors and separation tanks had to be used.
Simulation:
We carried out simulations of the process to determine the data required to size our
equipment. We chose the SRK and Lee Kesler methods as the appropriate thermodynamic
models to bring the reactants to our desired temperature and pressure (200 degrees and
85barg respectively) in the early stages of the simulated process through multiple
pressurisation and cool-down of the filler.
In order to achieve energy savings and less waste of feedstock, we included a recycling
loop in the simulation and thermally integrated a suitable heat exchanger. Through several
simulations we ended up with a methanol product that met the product requirements.
PFD and PID:
Results
Final Product
We obtained at the end of our analysis the final product is composed of 63 % methanol and
9.23.10^-7 ppmwt of. Therefore, the methanol produced respects the specification required.
Costs analysis
Fixed Capital
Equipements ISBL ($) Investment ($) Share (%)
Vessels 100.461,94 401.847,76 4,05%
Groupe 1.9 - Academic Year 2022/2023 5
Reactor and catalyst 329.661,90 329.661,90 3,32%
Pumps 7.153,80 21.461,40 0,22%
Compressors 2.860.843,08 8.582.529,24 86,53%
Heat Exchangers 194.430,36 583.291,07 5,88%
Total: $3.492.551,08 $9.918.791,37 100,00%
As we can see, compressors are the most expensive equipment. The number of compressors
must then be optimized.
CO2 emissions
CO2 emissions due to the production of H2 by electrolysis represent 74% of the global
emissions. However, the total footprint is less than the classic methanol process (
2. INTRODUCTION
Based on bibliographic research, for the synthesis reaction of methanol, the ratio of
reactants hydrogen to carbon dioxide is approximately 3:1 and that the production is
accompanied by the production of carbon monoxide as a by-product. The following equation
shows the main reaction equation (1) and the side reaction equation (2):
CO2 + 3H2 = CH3OH + H2O ΔrH = -49,51 kJ/mol (1)
CO2 + H2 = CO + H2O ΔrH = 41,17 kJ/mol (2)
Moreover, we also understand from the literature consulted that favorable operating
conditions include reaction temperatures in the range of 200-300°C and reaction pressures in
the range of 50-100 bar. Since the primary reaction is an exothermic reaction, low
temperatures are favorable for the reaction equilibrium to advance. This is because the lower
ambient temperature allows the system to dissipate heat more efficiently and avoid
overheating; the sum of the gaseous reactant coefficients of the main reaction is greater than
the sum of the gaseous product coefficients, and according to Le Chatelier's Principle, high
pressure drives the reaction equilibrium in the direction of the products. A relatively low
temperature and high-pressure reaction environment is therefore our first choice.
To make the reaction progress more easily, copper/zinc oxide compounds are often
added as catalysts, and sometimes additives such as Al, Zr, Cr, and Ga are also added to
improve the reaction performance.
For the collection of reaction materials, hydrogen is obtained by electrolysis of water
molecules in an electrolytic cell using renewable energy sources such as solar photovoltaic or
wind energy, which decomposes the water molecules and forms hydrogen and oxygen. For
carbon dioxide, we use captured carbon dioxide. The most common technologies for
capturing carbon dioxide are chemical absorption, physical separation, membrane separation,
etc.
Methanol also has a great variety of industrial applications: as a versatile compound,
methanol can be used to produce chemical derivatives such as formaldehyde, acetic acid,
methyl tertiary butyl ether (MTBE), and dimethyl ether (DME). In addition, methanol also
serves as a fuel substitute, for example by blending with gasoline to produce DME and/or
biodiesel.
2.3. PROCESS PRINCIPLE
First of all, it is essential to select the production process as the basis of the project, as
this will be related to the selection of the type and amount of raw materials, the
determination of the reaction conditions, and a number of other subsequent processes.
According to the documentation provided, the process chosen is a carbon dioxide
hydrogenation process, which is depicted in the following simplified diagram:
3. PROCESS CALCULATIONS
3.1. PROCESS SIMULATION
For the simulation of processes, we used PROII, a simulation software that allows us
to simulate and analyze various process operations and process parameters. PROII's
simulation data also allows us to build up mass and material balances and design equipment
sizes using.
The main objective of the simulation is to optimize the process to improve the
efficiency and performance of the equipment by varying the simulation parameters so that
the methanol product meets the specifications. In addition to this, we are required to design
all process equipment based on the simulation data obtained (e.g., temperature, pressure,
flow rate) and to estimate the energy consumption generated and the overall cost of building
the process.
Standardised, uniform units can ensure the accuracy of calculations, improve the
efficiency of simulations, facilitate comparison of results and exchange of data, as well as
ensuring consistency and reproducibility of modelling.
Figure 4. O2 Feedstock
By setting a battery limit condition, the boundaries of the system design can be
clarified, helping to ensure that important components are not overlooked or missed by our
group during the project. Additionally, the battery limit condition ensures consistency and
control throughout the engineering design process. Once it has been set, each group must
design the system within this condition, ensuring coordination and consistency between the
various components.
According to the requirements for the design of the simulated compressor, we needed
to comply with the basic condition that the compressor outlet temperature should not exceed
140 degrees Celsius. In order to meet this limit, we have placed separate heat exchangers
before the compressors to cool down the heat in time. In addition to this, as all our fluids are
gaseous in the pre-process stage, each heat exchanger has to be set at a pressure drop of 0.3
bar according to the Process Design Criteria.
3.1.6. Reactor
To focus on the process simulation, design and optimization, reactors will be
considered as a Gibbs reactor. According to our research, this should be because the Gibbs
reactor allows the use of different reaction kinetic models to describe the rate of chemical
reactions. Furthermore, it can consider the transfer of substances and energy between the
different phases (assuming that some of our products are liquid) to better describe and
optimise the reaction process in terms of transfer phenomena. The following table shows the
basic conditions of the reactor:
Figure 11. Relation between Methanol Product Mass and Recycle Rate
3.1.8. Vessels
In our simulations, we have used the vessel for separating liquids and gases several
times. They are shown in the following diagrams:
The equipment sizing is divided into six parts, each one of which corresponds to a
different equipment category: vessels, piping, reactor, heat exchanger, compressor, and pump.
Each equipment category has its own characteristic dimensions that are interesting for the
calculation of the final price but also for the proper functioning of the process. Together with
the equations that were given, the process simulation and the PFD the sizing process is
described below.
In each equipment category, a design temperature and pressure were used. They are
obtained by different relations. Regarding the design temperature, the following relation is
used.
𝑇𝐷 = 𝑇𝑀𝐶𝑂 + 15℃ (3)
Equation 3. Design Temperature
Regarding the design pressure there are two different relations: one for the equipment
(vessels, piping, reactor, heat exchanger and compressor) and the other one for the discharge
of pumps.
In this process, there are six vessels in total represented in the simulation, one of them
(V01 as its PROII name) is only a simulation tool, so there is no need to represent it neither in
the PFD nor in the sizing since it will not be used in the real process plant. Therefore, five
vessels were sized, four of them are vertical vessels and one of them is horizontal, the last one
is not represented in the PRO II Simulation, since it is only a real-life tool for the process
functioning.
The vessel sizing starts with the specification of the design temperature and pressure,
then followed by the identification of the inlets flows in each vessel. One of the vertical vessels
(C-002) had two inlets, so the parameters were considered as the ratio of the parameters of
each inlet. Other simulation parameters such as liquid fraction, volumetric flow, liquid and
vapor density were also identified to its future utilization in the calculations. In the following
tables, these parameters are presented for each equipment. Also, it was identified the vessels
that needed to have Mash Pad or not, since this is a key information for some sizing
calculations.
Subsequently, it is important to establish a critical vapor velocity for each vessel which
will guide all the other size parameters. This critical vapor velocity is calculated by the
following equation. But, depending on the presence or absence of Mash Pad this velocity will
be multiplied for a security factor of 2.2, if there is Mash Pad, or 0.8, if there is not. Except for
the horizontal vessel, it is necessary to multiply its critical vapor velocity by a factor of 170%.
Moreover, another difference between the horizontal vessel that differs from the vertical
vessel is that its final diameter must be increased by a factor of security that corresponds to a
HHLL and LLL marge.
𝜌𝐿 −𝜌𝐺
𝑣𝑐 = 0,048√ (4)
𝜌𝐺
The volumetric flow is multiplied also by a security factor of 1.1 and after that it is used
for calculating the cross-sectional area of the Vessel with the following equation to obtain an
initial value of the area and consequently an initial value of the diameter.
𝑚3
𝑚 𝑄( )
𝑠
𝑣𝑐 ( 𝑠 ) = 2
(5)
𝐴 (𝑚 )
Equation 5. Critical vapor velocity related to the vessel flow and area.
𝜋𝐷 2
𝐴= (6)
4
Once we obtain the initial diameter, it is possible to obtain the first estimation of the
equipment height, knowing that the maximum height is equal to four times the diameter of
the equipment and the minimum height is equal to two times the diameter of the equipment.
Other parameters calculated in this part are the corrosion allowance, the maximum
allowable stress value and the equipment thickness, since they are parameters that depend
on the type of flowing process, the design temperature and the diameter of the vessels. The
following table gathers these parameters.
Table 2. Corrosion allowance, the maximum allowable stress value and equipment thickness.
To calculate the final diameter of each vessel, it is crucial to consider the pressure drop
in the lines connected with them but also the requested liquid level control. Therefore, it is
important to find a new diameter using the initial diameter and incorporating those other
pressure and liquid level constraints. An interactive process is done to find the adequate
diameter that can respect the conditions of maximum and minimum height (between two and
four times the diameter of the equipment).
This interactive process consists of choosing an initial diameter (for most vessels we
choose 1m of diameter), which will be used together with the found value of nozzle diameter
(calculated in the piping section) to calculate the small heights level control (h1, to h7) that
3.2.2. PIPING
The Piping Sizing is important for the Vessel Sizing as well as the PID representation.
All the Piping appearing in the PID and those that were considered as the inlets of the vessels
were sized, therefore, fifteen lines should be sized but since some of them were referred to
the same flows in the simulation, only 9 lines were sized. The lines were divided into liquid,
mixed and vapor lines, each line type has its own categories of classification that also changes
the calculation that needed to be done. In the following table, the lines names together with
its classifications are resumed.
For the liquid lines, the sizing starts with the choice of the velocity, with which the
cross-sectional area is found and subsequently the line diameter can be obtained. For this
reason, the liquid density and volumetric flow in each line is needed and they can be obtained
from the PROII Simulation. The velocity and the diameter are associated with each other.
Thus, it is possible to relate the value of the velocity and the diameter until they meet
the conditions of maximum velocity and corresponding diameter present in " Process Design
Criteria". Once the velocities are chosen, the pressure drop calculation takes place and the
diameter chosen will be readjusted until an acceptable value of pressure drop is reached. For
the pressure drop in liquid lines the equations used were those indicated for single phase
lines. The following equations were used and a value of 0,0018" was used for the absolute
roughness of the carbon steel pipe.
𝐿 𝐺2
∆𝑃 = 𝐹𝐷 . 𝐷 . 2𝜌 (7)
1
8 12
12
1
𝐹𝐷 = 8. [(𝑅𝑒) + 3 ] (8)
(𝐴+𝐵)2
37530 16
𝐵=( ) (10)
𝑅𝑒
Although the pressure drop does not correspond exactly to the provided normal and
maximum pressure drop, the results were considered acceptable since they are in the same
order of magnitude as the ones requested. This could also be due to the diameter chosen
having to be standardized into integers and the inlet flows being relatively small. The results
found for the liquid lines were gathered in the following table.
For the mixed lines, the vapor fraction was obtained from the PROII Simulation to
calculate the pressure drop mixed. Besides the condition of the pressure drop for mixed lines
also must respect a condition between the value of the density and velocity of the fluid, as
shown by the following relations:
2
5000 𝑃𝑎 < 𝜌𝑚 𝑣𝑚 < 15000 𝑃𝑎 (11)
As for the classifications of the liquid lines, the classification Partial Condenser Outlet
(General liquid) was chosen for the outlet lines of heat exchangers and the category of Mixed
Phase Condensates was chosen for the rest of the lines. The calculations of the pressure drop
for the mixed line were based on the given model of pressure drop for two phase flow. This
model is based on the following equations.
𝜌 𝜌
𝐶2 = √𝜌 𝐿 + √ 𝜌𝐺 (13)
𝐺 𝐿
𝜌𝐿 .𝜌𝐺
𝜌𝐻 = 𝑥.𝜌 (15)
𝐿 +(1−𝑥).𝜌𝐺
32.𝑌
𝐶1 = 2 + (16)
1+1,005664.𝐶𝐹0,8
3
𝑌 = [1 − [0,16. (2,5 + log(𝑋))2 ]] (17)
𝜌𝑉 𝜇 0,2
𝑋 = max [0,00316; . ( 𝐿) ] (18)
𝜌𝐿 𝜇𝑉
(300−𝐺)2
𝐶𝐹 = 𝐺 , if 𝐶𝐹 < 300, then 𝐶𝐹 = 300 + (19)
40
The priority established was to obey the maximum velocity in the line and then with
the resulting diameter calculate the pressure drop until a reasonable value was found. Lastly,
the velocity was adapted to follow the condition of integer value to the diameter size. The
final pressure drops were considered high in comparison to the references, but this can be
justified by the fact that the flow rating was very low in these lines, so it was necessary to
establish a low velocity to a minimum pressure drop and a coherent piping diameter. The
following table resumes the final values for the chosen velocities and pressures drops.
For the vapor lines, the lines were categorized as normal gas lines since they were not
related to any compressor, columns, reboilers or steam lines. The vapor lines were classified
as pump suction bubble point, because the line presents gases dissolved in it, and pump
discharge (with a pressure smaller than 50 bar), because of the location and the pressure of
the line. Once the gas lines do not have pressure drop restrictions, the only calculation made
The adequate velocity was chosen to obey the condition above and to obtain an
integer value of diameter. In the table below, the chosen velocities and diameters are
represented.
3.2.3. REACTOR
Above all, it is necessary to obtain the volumetric feed of the reactor PROII Simulation
since it is the only parameter needed for the calculation. This value was equal to the
volumetric inlet flow of the line INLET_FIN. With this inlet flow, it is possible to calculate all
the parameters of interest, such as the diameter, length, and reactor volume. The equations
used are described below.
𝑄
𝑉 = 𝐺𝐻𝑆𝑉 (21)
Equation 21. Reactor volume and GSHV as being Gas Hourly Space Velocity (1/h).
The final values are resumed in the following table. We used 3 as the factor between
the length and the diameter, which could be chosen between 2 and 4. Although the final value
of 2 m³ for the whole reactor is not acceptable and relatively small for the parameters of an
industry, it is possible to justify this value with the small volumetric inlet in comparison to the
required GHSV given.
There are six heat exchangers in the PROII Simulation but only five of them are really
concerned with industrial installation, in the PFD and PID, thus only five heat exchangers were
sized.
Firstly, each heat exchanger is identified with its utility fluid, which was chosen
whether the heat exchanger was a condenser or a reboiler and whether the working
temperature was close to those of the utility’s fluids. For the equipments E-001, E-002, E-003,
and E-005 cold water (CW) was used as the cold fluid and for the E-004 medium pressure
steam (MP Steam) was used as the hot fluid.
The minimum approach temperatures were given, and they were inserted into the
PROII Simulation to obtain more information about the heat exchanger parameters, and they
were also used to obtain the outlet temperature of the cold side utility. With this information,
it was possible to obtain the correction factor, Ft, the logarithmic temperature difference,
LMTD, and the energy flow of each heat exchanger.
Next, the heat exchangers were classified according to their service and utilities fluids
to obtain the heat transfer coefficient, U. The E-001 and E-002 were classified as Inert Vent
Gas Chiller, the E-003 was classified as Feed Reactor/Effluent reactor, the E-004 was classified
as Effluent Reactor / Saturated steam and the E-005 was classified as Effluent reactor/CW.
Once we obtain all the parameters, we can calculate the area of heat transfer by the following
equation.
𝑄 = 𝑈. 𝐴. 𝑀𝐷𝑇 (22)
3.2.5. COMPRESSOR
There are three compressors in total in the present process and all of them are used
in the real-life industrial plant. The first information needed is the presence or absence of
recycling flow in the compressor inlet. Afterward, the compressor inlet line must be identified,
and its temperature, shaft work, pressure, and volumetric flow are obtained from the PROII
Simulation and noted. If there is a recycling flow in the current equipment, then the
volumetric flow will be multiplied by a factor of 5%.
With all the updated values of the volumetric flow and the following equation, the
efficiency of each compressor is obtained.
𝑄 𝑄 2
𝑉 𝑉
𝜂 = 0,751 + 0,01985. ln (8500) − 0,02. (ln (8500)) (23)
The final parameters for each compressor are represented in the following table.
3.2.6. PUMP
In the current process, only two pumps were needed in the final line of the methanol
product but only one pump was sized since the two of them have the same characteristics.
For the sizing of the pump, some parameters of the PROII Simulation are needed, such as:
volumetric flow, inlet line vapor pressure, pump pressure drop, efficiency, differential head,
and density of the inlet flow. Other parameters were obtained from the vessels and piping
sizing, namely the design temperature and pressure, the LLL of the vessel C-005, the pressure
drop in the suction and the difference of pressure in the suction and discharge lines. Still, the
water density was obtained from the literature as equal to 997 kg/m³ [XX]. The parameters of
control equipment, pressure drop, pump suction and discharge height were given. All the
values used for the calculations are resumed in the following table:
The volumetric flow must be multiplied for a correction factor of 1.1. Other parameters
such as the liquid minimum heigh (∆𝐻𝑠 ) are obtained from the C-005 vessel sizing and
corresponds to the LLL plus 1 m as suggested. Also, in some calculations it is considered the
relative density, that is equal to the liquid density divided by the water density.
Two different methods were considered in the pump sizing, one method using only
the equations and relations given (Calculation 1) and the other using also the notions of
Bernoulli pressure drop (Calculation 2). Each method will be explained in the sequence.
3.2.6.1. Calculation 1
In "Calculation 1", all the parameters available – such as difference of pressure in the
pump, volumetric flow and pump efficacity - in the PRO II Simulation were used. With them,
it is possible to calculate the power with the equation 24. This value is only a basis value to
find the nominal power, that is found by following the instructions in "Process Design Criteria".
Since it is smaller than 22 kW, a factor of 1,25 is multiplied to it to find the nominal power.
With the nominal power it is possible to find the performance, by the values on the table X.
Primary to the NPSH calculation, it is necessary to calculate the ℎ𝑠𝑚𝑖𝑛 and the ∆𝐻𝑠
with the following equations.
(25)
(26)
Equation 26. Height suction variation.
(27)
Equation 27. Pressure drop suction.
Then the ℎ𝑠𝑚𝑖𝑛 is calculated. Finally, the NPSH can be obtained with the equation
below. For the 𝑃𝑠𝑣𝑚𝑖𝑛 , is ℎ𝑠𝑣𝑚𝑖𝑛 used.
(28)
Equation 28. Pump NPSH.
In this calculation procedure, the pump performance was very low, smaller than 70%,
since the nominal power was also very low. Although, the NPSH obtained had a value bigger
than 3 m, that means It presented no cavitation. These results are coherent with the PROII
Simulation values, from where the parameters were obtained. It is important to notice that
the pump sized was not crucially important to the process in terms of power availability
because of the great pressure that the vessel C-005 and the lines in the suction were already
submitted.
3.2.6.2. Calculation 2
In "Calculation 2", the pressure drop of every control equipment was considered, as
well as the efficiency obtained from the correlation shown in figure above and the pressure
drop at the pump discharge. The pressure calculation was based in Bernoulli's model, as the
following equation 29 shows, although the velocity term was neglected.
1 1
𝑃1 + 2 . 𝜌. 𝑉12 + 𝜌. 𝑔. ℎ1 = 𝑃2 + 2 . 𝜌. 𝑉22 + 𝜌. 𝑔. ℎ2 (29)
For the power calculation the same equation was used, but the pressure drop, and the
efficiency was calculated in a different manner. Starting with the pressure drop calculation,
first, it was obtained separately, the values of the suction and the discharge pressure. The
equations 30 to 31 were used.
The efficiency was found by the values of volumetric flow and the pump head, using
figure above, the value found was 51.7 %. Afterwards, the ∆𝐻𝑠 is calculated with the same
equation as in "Calculation 1" as well as ℎ𝑠𝑚𝑖𝑛 . Lastly, the NPSH is calculated by the same
equation.
The nominal power obtained in this calculation method is significantly greater than
that found in the first calculation method which results in higher performance as well (85%).
However, in this case the results are not coherent with the PRO II Simulation results.
Moreover, the NPSH obtained was a negative number, which has no physical meaning. That
problem was due to the value of vapor pressure being greater than the pressure in the C-005
Vessel.
The final parameters for the pump are represented in the following table.
4. PROCESS ANALYSIS
4.1. HEALTH, SAFETY AND ENVIRONMENT
4.1.1. HEALTH
Methanol is a colorless, volatile liquid that can be toxic when ingested, inhaled or in
contact with the skin. It is therefore essential to implement personal protection measures,
such as appropriate protective equipment (gloves, masks, safety glasses), as well as safe
handling procedures to avoid direct contact. The actions that can be taken for the health are:
Workers training: Workers must be informed of the health risks associated with methanol,
including exposure by inhalation, skin contact, or accidental ingestion. They must also be
trained in good methanol handling practices, including the use of personal protective
equipment (PPE) such as gloves, safety glasses and protective suits.
Exposure control measures: Procedures should be in place to minimize worker exposure to
methanol, such as adequate ventilation of work areas, and proper storage and handling of the
product. implementation of containment and vapor capture systems, as well as air quality
monitoring.
Personal hygiene: Adequate sanitary facilities must be available for workers, including hand-
washing facilities and showers.
The table below shows the causes, consequences and actions that can be taken for the safety
of the personnel.
- Regular inspection of
equipment, piping and valves to
identify and repair leaks.
- Use of chemical-
resistant materials and seals
quality seals.
Exposure to Health problems - Use of adequate
methanol vapors ventilation to remove vapors of
methanol.
- Use of personal
protective equipment (PPE)
such as breathing masks, safety
glasses and overalls.
- Training workers on
the health risks associated with
methanol exposure and
preventive measures.
Fires and Injuries, loss of - Implementation of
explosions life fire detection and suppression.
- Worker training on
evacuation procedures and
measures to combat fires.
- Regular maintenance
of fire protection equipment.
Chemical Spills, toxic -In-depth assessment
unwanted reactions discharges of potential chemical reactions
and product incompatibilities
chemicals used.
- Proper storage and
handling of chemicals to avoid
spills.
- Training workers in
reaction management
procedures chemicals.
Isolation valves were also placed to isolate each device of the process. This ensures safety for
the employee in case of maintenance. Furthermore, we connected another pump in parallel
to the main one in case of maintenance needs.
6.1.3 ENVIRONMENT
Air and water pollution Deterioration of air - Use of emission control systems to reduce
and water quality atmospheric emissions of pollutants
Comparison between CO2 emissions from normal methanol synthesis and e-methanol
synthesis of our process:
CO2 emissions from other processes (source : Carbon footprint evaluation of coal-to-methanol chain
with the hierarchical attribution management and life cycle assessment Zhen Qin a,b , Guofu Zhai c , Xiaomei
Wu a,b , Yunsong Yu a , Zaoxiao Zhang)
It can be noticed that the CO2 produced from our process is less than the basic coal-
to-coal chain process. Furthermore, one solution is to install a secure CO2 capture and storage
system, then reinject it into the process. This reduces the impact on the environment.
The calculation of the cost of the flashes was estimated based on the thickness. To
calculate it, we use the following relation:
∅
𝑃×( +𝑐)
2
𝑒 = 𝑆×𝐸 −0,6×𝑃 + 𝑐 (32)
Equation 32. Thickness
With: P [MPa]: Design pressure for vessel; E [-]: E=1; ∅ [mm]: diameter of each vessel;
c [mm]: corrosion allowance; S [MPa]: Maximum constraint.
Maximum allowable stress values, S, exist for several typical materials as a function of
design temperature. As the vessels were at a design temperature between 55 and 222 °C,
according to Process Design Criteria (figure 1 - Working limits for steels in hydrogen service to
avoid high temperature hydrogen attack), the appropriate material for the vessels is Carbon
steel. The choice of Carbon steel as a material has a positive impact on the economic analysis,
since between Cr-Mo steel and Stainless steel it presents the most advantageous cost-benefit.
Next, the thickness was used to find a relationship that would provide the volume
value of each vessel. Knowing the density of the carbon steel and the volume of the
equipment it was then possible to calculate the mass. Finally, the carbon steel price ratio given
in the PDC was used to calculate the price of the equipment and the installation price.
Fixed Capital
Equipment Name ISBL ($) Investment ($)
Vessel C-001 2.421,00 9.683,99
Vessel C-002 34.775,61 139.102,44
Vessel C-003 42.264,52 169.058,09
4.2.1.2. Reactor
The reactor cost is estimated considering the reference cost of 150€/m3. The volume
of the reactor was estimated at 2.06 m3 during sizing, which results in a cost of about €
308,5M. Considering the currency conversion of July 1th 2023, €1 = $1.07, the cost of the
reactor is $329M, which represents 2.52% of the total project investment.
4.2.1.3. Pump
To estimate the value of the pumps, the following relation was used:
𝑃 0,7
𝐶 = 800 × (1000) (33)
Equation 33. Pump cost
With: C: Equipment cost in k$; P: Shaft Power in kW
4.2.1.4. Compressor
To estimate the value of the compressor, the reference cost of €0.5 million/ MW and
the Power (MW) is used. Currency conversion from euros to dollars resulted in a total
equipment cost of $2,9M, resulting in a fixed capital investment of $8,6M. The total
expenditure for compressors represents the highest expenditure for the project at 86.5%.
𝐴 0,8
𝐶 = 150 (1000) (34)
Equation 34. Heat exchangers costs
With: C: Heat exchanger cost in k$; A: Area in m2;
When the design pressure is higher than 50 barg, the following factors were considered:
Thus, for the 5 heat exchangers, the total cost is $194M, resulting in a final cost of
$538,3M, which is 5.9% of the total project cost.
Fixed Capital
Equipment Name ISBL ($) Investment ($)
Heat Exchanger E-001 9.588,53 28.765,59
Heat Exchanger E-002 25.864,02 77.592,07
Heat Exchanger E-003 12.981,63 38.944,90
Heat Exchanger E-004 27.262,65 81.787,94
Heat Exchanger E-005 118.733,52 356.200,56
Total: $ 194.430,36 $ 583.291,07
Table 16. Heat exchangers costs
Finally, the total cost of the e-methanol production plant was estimated at $9,9M. The
description of the values of each equipment was made in the following table:
Groupe 1.9 - Academic Year 2022/2023 42
Fixed Capital Share
Equipements ISBL ($) Investment ($) (%)
Vessels 100.461,94 401.847,76 4,05%
Reactor and catalyst 329.661,90 329.661,90 3,32%
Pumps 7.153,80 21.461,40 0,22%
Compressors 2.860.843,08 8.582.529,24 86,53%
Heat Exchangers 194.430,36 583.291,07 5,88%
Total: $3.492.551,08 $9.918.791,37 100,00%
Table 17. CAPEX results
To better illustrate the data below, we have also created a graph illustrating the
distribution of costs:
4.2.2. OPEX
Thus, a value of $159M was obtained for the utilities, which represents 38,6% of
the overall OPEX value.
• Raw materials (feed):
To calculate the amount spent on feed, the following expression was used:
With this, we had a value of $247M which represents the largest portion of the OPEX
expenses, 60%.
4.2.2.2. Fixed cost
For the fixed cost, expenses with labor, supervision, quality control, maintenance
labor, maintenance material, and operating supplies were considered.
• Operating labor:
To determine the labor expenses, it is necessary to determine the amount of salary and
the amount of employees needed to operate in the plant for this the following equation
was used:
0,5
𝑁𝑂𝐿 = (6,29 + 31,7𝑃2 + 0,23 × 𝑁𝑛𝑝 ) (38)
Equation 38. Number of operators
With: NOL: number of operators per team; P: number of process steps (unit operation);
Nnp: number of steps of type: compression, heating/cooling, mixing, separation, and
reaction.
In the synthesis of e-methanol, there is one unit operation and three compression stages,
five heating stages, two mixing stages, four separation stages, and two reaction stages, for
Thus, the total salary cost for the project is $659M, which represents 0,16% of the OPEX
expenses.
• Other expenses:
The following formulas were used to determine the other expenses included in the fixed
costs of OPEX:
𝑂𝑝𝑒𝑟𝑎𝑡𝑖𝑛𝑔 𝑆𝑢𝑝𝑒𝑟𝑣𝑖𝑠𝑖𝑜𝑛 = 0,20 × 𝑂𝑝𝑒𝑟𝑎𝑡𝑖𝑛𝑔 𝐿𝑎𝑏𝑜𝑟 𝐶𝑜𝑠𝑡 (40)
𝑄𝑢𝑎𝑙𝑖𝑡𝑦 𝐶𝑜𝑛𝑡𝑟𝑜𝑙 = 0,20 × 𝑂𝑝𝑒𝑟𝑎𝑡𝑖𝑛𝑔 𝐿𝑎𝑏𝑜𝑟 𝐶𝑜𝑠𝑡 (41)
𝑀𝑎𝑖𝑛𝑡𝑒𝑛𝑎𝑛𝑐𝑒 𝐿𝑎𝑏𝑜𝑟 = 0,027 × 𝐹𝑖𝑥𝑒𝑑 𝐶𝑎𝑝𝑖𝑡𝑎𝑙 𝐶𝑜𝑠𝑡 (42)
𝑀𝑎𝑖𝑛𝑡𝑒𝑛𝑎𝑛𝑐𝑒 𝑀𝑎𝑡𝑒𝑟𝑖𝑎𝑙 = 0,018 × 𝐹𝑖𝑥𝑒𝑑 𝐶𝑎𝑝𝑖𝑡𝑎𝑙 𝐶𝑜𝑠𝑡 (43)
𝑂𝑝𝑒𝑟𝑎𝑡𝑖𝑛𝑔 𝑆𝑢𝑝𝑝𝑙𝑖𝑒𝑠 = 0,0075 × 𝐹𝑖𝑥𝑒𝑑 𝐶𝑎𝑝𝑖𝑡𝑎𝑙 𝐶𝑜𝑠𝑡 (44)
Through the analysis it is possible to verify that the highest operating costs are
associated with the purchase of materials, which reflect in a share of 60% of total costs.
Finally, it was possible to establish that the plant has annual operating costs of $413M.
𝐶𝐴𝑃𝐸𝑋𝑁𝐸𝑊 +𝑂𝑃𝐸𝑋
𝑃𝐶 = (46)
𝑃𝑅𝑂𝐷𝑈𝐶𝑇𝐼𝑂𝑁
Equation 46. Production cost
Production Cost
CAPEX ($) 9.918.791,37
fa (10%/year, 20 years) 0,12
CAPEXnew ($) 1.165.057,51
OPEX ($) 413.004.322,60
Methanol Product (t/h) 22,80
Annual Product (t/year) 182.360,00
Production cost ($/tMeOH) $2.271,16
Table 19. Production cost
Thus, the total expenditure for production is $2,273 per ton of e-methanol. This value
means that the product must be sold at this price to ensure that the factory does not make
losses or profits after 20 years.
• Pump
o Backup pump
Stopping the production for maintenance can result in a significant loss in revenue.
Therefore, a backup pump was installed and will be used when the other is having
maintenance. For each pump, a system of drainage and isolation valves is added to ease
maintenance.
o Filters
A filter is installed at the pump inlet to obstruct and retain any unwanted substances
such as residues, particles, and erosions from the piping thus protecting the pump from
damage.
o Safety system
A safety system (HSL) is activated if the pump speed is too low, by sending an electric
signal to the pump motor.
o Sensors
The suction and discharged pressure can be monitored by a pressure gauge (PG) which
is a local indicator. This will help the operator to identify any unusual behavior on the pump.
We find that our methanol is calculated to be priced at approximately $2200 per ton
compared to commercial methanol at approximately $400 per ton. This indicates that there is
still room for improvement in our process and indirectly proves that our calculations contain
a large margin of error.
5.2. Conclusion
The goal of the project was to study and simulate the process of industrial methanol
production. Through an in-depth study of methanol synthesis reactions, operating conditions,
and other aspects, we were able to understand and optimise the key factors in the methanol
production process. The following is a conclusion of the project:
In the early stages of the project, we conducted extensive literature research to
understand the concept, preparation process and application areas of e-methanol. We
learned that methanol is a versatile compound that can be used not only as a raw material for
the production of chemical derivatives but also as a fuel substitute. This provided the basis for
defining our project objectives and approach.
During the simulation, we used PROII software for process simulation and optimisation.
This software provides powerful modelling tools and calculation capabilities that can help us
to analyse and optimise key parameters in the process. PROII simulations enabled us to predict
and evaluate the product quality and capacity under different operating conditions and to
rationalise our process with the aim of increasing process performance, improving yields, and
saving energy.
In the middle of the project, we studied and tried to make PFD and PID, which helped
us to gain an insight into the basic principles and procedures of industrial processes. Learning
about the role and interrelationship of different equipment and instrumentation in the
process has also helped us to develop an overall understanding of the engineering field.
6. BIBLIOGRAPHIC REFERENCES
XI. YU, Minli, WANG, Ke et VREDENBURG, Harrie, 2021. Insights into low-
carbon hydrogen production methods: Green, blue and aqua hydrogen.
International Journal of Hydrogen Energy. [en ligne]. juin 2021. Vol. 46,
n° 41, pp. 21261‑21273. [Consulté le 5 février 2023]. DOI
10.1016/j.ijhydene.2021.04.016.
XII. LEI, Hong, Zhaoyin HOU et Jianwei XIE (jan. 2016). « Hydrogenation of
CO2 to CH3OH over CuO/ZnO/Al2O3 catalysts prepared via a solvent-
free routine ». In : Fuel 164.Supplement C, p. 191-198
XIII. STEINBERG, Meyer. United States Patent Steinberg <<
ELECTROLYTIC SYNTHESIS OF METHANOL FROM CO >>.
Huntington Station, 1976, N.Y.
XIV. HOSEINY, S., ZARE, Z., MIRVAKILI, A., SETOODEH, P. et
RAHIMPOUR, M.R., 2016. Simulation–based optimization of operating
parameters for methanol synthesis process: Application of response
surface methodology for statistical analysis. Journal of Natural Gas
Science and Engineering. [en ligne]. août 2016. Vol. 34, pp. 439448.
[Consulté le 7 mars 2023]. DOI 10.1016/j.jngse.2016.06.075.
XV. Dieterich, Vincent & Buttler, Alexander & Hanel, Andreas & Hartmut,
Spliethoff & Fendt, Sebastian. (2020). Power-to-liquid via synthesis of
methanol, DME or Fischer–Tropsch-fuels: a review. Energy &
Environmental Science. 13. 10.1039/D0EE01187H.
XVI. Hank, Christoph, et al. ‘Economics & Carbon Dioxide Avoidance Cost of
Methanol Production Based on Renewable Hydrogen and Recycled
Carbon Dioxide – Power-to-Methanol’. Sustainable Energy & Fuels, vol.
2, no. 6, May 2018, pp. 1244–61. pubs.rsc.org,
https://doi.org/10.1039/C8SE00032H.
PDS
Numéro doc. Client : Page
CLIENT
Revision
SITE
BALLON
UNITE C-005
1 REPERE
2 SERVICE
3 SERVICE CALCUL PROCEDE
4 Pression 17.5 barg 19.3 barg
5
PI T E
6 (1)
7 Température 49.5°C 64.5°C
8
LT
9
10 (1)
11 VIROLE/FONDS INTERNES (2)
12 Materiaux carbon steel carbon steel RD
Eliminateur
14
19 Calorifuge no
H
20 Revêtement
21 Elévation LT inf. au-dessus du sol = 2.42m mini
22 Notes F EXT.
23
24
25
26
28 REP Nb Ø SERVICE N
LT
29 A 1 6" Alimentation
30
31
F
32 D 1 2" Purge
33 E 1 2" Event
34 F 1 4" Sortie fond
35 H 1 8" Trou d'homme
36 L - - Ligne d'équilibre
37 N - - Nettoyage
38 S - - Soutirage
39 T 1 2" Sortie tête
40 V - - Entrée vapeur
41 RD 1 - Organe de sécurité A
r=
42 LG 1 - Niveau à glace
43 LT 4 - Transmetteur de niveau
44
45
46 TI 1 - Prise de température
47 PI 1 - Prise de pression (1) Pour les services concernés, préciser aussi les pressions partielles d'H2 et d'H2S
48 ainsi que les conditions de dépressurisation (T,P)
49 (2) Seulement en cas de demande Procédé specifique, specifier les internes et les
50 materiaux / surepaisseur de corrosion correspondant.