Patent Description:
Conventionally, there has been known a control technique in which an oxygen storage state in a three-way catalyst is grasped by a three-way catalyst provided in an exhaust pipe of an internal combustion engine and exhaust gas sensors assembled on the upstream and downstream sides of the three-way catalyst, and an air-fuel ratio of an air-fuel mixture is corrected according to a result thereof. In this control technique, a rich correction of the air-fuel ratio is determined on the basis of the oxygen storage state captured by the time integration of the amount of air taken into the engine and the product of the air-fuel ratio of the exhaust gas and the difference value between the theoretical air-fuel ratio. Further, feedback correction of the air-fuel ratio control is performed by detecting the presence or absence of oxygen released on the downstream side of the three-way catalyst by the exhaust gas sensor assembled on the downstream side of the three-way catalyst (for example, PTL <NUM>).

An internal combustion engine control device according to preamble of claim <NUM> is known from <CIT> Al.

Incidentally, there is a tendency that the frequency of motoring operation accompanied by fuel cut and engine stop due to idling stop increases due to the demand for reduction in fuel consumption of the internal combustion engine. Therefore, there is a problem that the oxygen storage state and temperature in the three-way catalyst cannot be accurately grasped, the purification efficiency of the three-way catalyst is deteriorated, and emission performance is deteriorated.

In the method described in PTL <NUM> in which the presence or absence of oxygen released on the downstream side of the three-way catalyst is detected using an exhaust gas sensor and feedback correction of the air-fuel ratio control is performed, the oxygen storage state of the catalyst reaches the lower limit value or the upper limit value of the target control range at the time when the presence or absence of oxygen is detected. Therefore, there is a problem that the catalyst purification efficiency cannot be kept high, and emission performance deteriorates.

In addition, there is a method in which the amount of oxygen stored in the catalyst during the motoring operation accompanied by fuel cut is estimated at the timing when the oxygen sensor output on the downstream side of the catalyst changes, and rich correction is performed at the time of restarting the firing operation on the basis of the fuel increase period according to the amount of oxygen. In this method, there is a problem that a temperature state of the catalyst and a transient change in an exhaust gas flow rate cannot be considered, an appropriate rich correction cannot be performed, and emission performance deteriorates.

The present invention has been made in view of the above circumstances, and an object thereof is to keep catalyst purification efficiency high and to prevent deterioration of emission performance. Solution to Problem.

In order to solve the above problems, an internal combustion engine control device according to an aspect of the present invention controls an internal combustion engine including a catalyst disposed in an exhaust pipe, a first exhaust gas sensor disposed on an upstream side of the catalyst, and a second exhaust gas sensor disposed on a downstream side of the catalyst. The internal combustion engine control device includes an oxygen storage ratio calculation unit that calculates an oxygen storage ratio of the catalyst based on a catalytic reaction model having at least a detection value of the first exhaust gas sensor as an input, a statistical model calculation unit that predicts a catalyst downstream exhaust gas concentration by using a statistical model having the oxygen storage ratio as an input and the catalyst downstream exhaust gas concentration as an output, and an air-fuel ratio correction amount calculation unit that calculates an air-fuel ratio correction amount of an air-fuel mixture of the internal combustion engine based on a future catalyst downstream exhaust gas concentration calculated by the statistical model calculation unit. The catalytic reaction model is defined by a reaction rate between at least oxygen on the upstream side of the catalyst, carbon monoxide on the upstream side of the catalyst, and a metal carried in the catalyst. The oxygen storage ratio calculation unit includes an oxygen concentration calculation unit that calculates an oxygen concentration on the upstream side of the catalyst based on a detection value of a first exhaust gas sensor on the upstream side of the catalyst; and a carbon monoxide concentration calculation unit that calculates a carbon monoxide concentration on the upstream side of the catalyst based on the detection value of the first exhaust gas sensor on the upstream side of the catalyst.

According to at least one aspect of the present invention, it is possible to keep catalyst purification efficiency high and prevent deterioration of emission performance.

Objects, configurations, and effects besides the above description will be apparent through the explanation on the following embodiments.

Hereinafter, examples of modes for carrying out the present invention (hereinafter, described as "embodiment") will be described with reference to the accompanying drawings. In the present specification and the accompanying drawings, components having substantially the same function or configuration are denoted by the same reference numerals, and redundant description is omitted.

First, a configuration example of an entire engine system to be controlled by an internal combustion engine control device according to an embodiment of the present invention will be described.

<FIG> is a schematic configuration diagram of an entire engine system to be controlled by an internal combustion engine control device according to an embodiment of the present invention. The engine system includes an internal combustion engine <NUM>, a flow rate sensor <NUM>, a turbocharger <NUM>, an air bypass valve <NUM>, an intercooler <NUM>, a supercharging temperature sensor <NUM>, a throttle valve <NUM>, an intake manifold <NUM>, a supercharging pressure sensor <NUM>, a flow enhancing valve <NUM>, an intake valve <NUM>, an exhaust valve <NUM>, a fuel injection valve <NUM>, an ignition plug <NUM>, a knock sensor <NUM>, and a crank angle sensor <NUM>. The engine system further includes a wastegate valve <NUM>, an air-fuel ratio sensor <NUM>, an exhaust purification catalyst <NUM>, an oxygen sensor <NUM>, an exhaust gas recirculation (EGR) pipe <NUM>, an EGR cooler <NUM>, an EGR valve <NUM>, a temperature sensor <NUM>, a differential pressure sensor <NUM>, and an electronic control unit (ECU) <NUM>.

An intake flow path and an exhaust flow path communicate with each other via the internal combustion engine <NUM>. The flow rate sensor <NUM> and an intake temperature sensor (not illustrated) built in the flow rate sensor <NUM> are assembled to the intake flow path. The turbocharger <NUM> includes a compressor 3a and a turbine 3b. The compressor 3a is connected to the intake flow path, and the turbine 3b is connected to the exhaust flow path. The turbine 3b of the turbocharger <NUM> converts the energy of the exhaust gas from the internal combustion engine <NUM> into the rotational energy of turbine blades. The compressor 3a of the turbocharger <NUM> compresses the intake air flowing in from the intake flow path by rotation of compressor blades connected to the turbine blades.

The intercooler <NUM> is provided on the downstream side of the compressor 3a of the turbocharger <NUM>, and cools the intake air temperature of the intake air that has been adiabatically compressed and increased by the compressor 3a. The supercharging temperature sensor <NUM> is assembled on the downstream side of the intercooler <NUM> and measures the temperature (supercharging temperature) of the intake air cooled by the intercooler <NUM>.

The throttle valve <NUM> is provided on the downstream side of the intercooler <NUM>, and narrows the intake flow path to control an intake air amount flowing into a cylinder of the internal combustion engine <NUM>. The throttle valve <NUM> is configured by an electronically controlled butterfly valve capable of controlling a valve opening independently of an accelerator pedal depression amount by a driver. The intake manifold <NUM> to which the supercharging pressure sensor <NUM> is assembled communicates with the downstream side of the throttle valve <NUM>.

Note that the intake manifold <NUM> provided on the downstream side of the throttle valve <NUM> and the intercooler <NUM> may be integrated. In this case, since the volume from the on the downstream side of the compressor 3a to the cylinder can be reduced, it is possible to improve the responsiveness and controllability of acceleration/deceleration.

The flow enhancing valve <NUM> is disposed on the downstream side of the intake manifold <NUM>, and enhances the turbulence of the flow inside the cylinder by generating a drift in the intake air sucked into the cylinder. When the exhaust gas recirculation combustion described later is performed, the turbulent flow combustion is promoted and stabilized by closing a flow enhancing valve (not illustrated).

The internal combustion engine <NUM> is provided with an intake valve <NUM> and an exhaust valve <NUM>. Each of the intake valve <NUM> and the exhaust valve <NUM> has a variable valve mechanism for continuously changing a phase of valve opening and closing. An intake valve position sensor <NUM> and an exhaust valve position sensor <NUM> for detecting the on-off phase of the valve are respectively assembled to the variable valve mechanisms of the intake valve <NUM> and the exhaust valve <NUM>. A cylinder of the internal combustion engine <NUM> includes a direct fuel injection valve <NUM> that directly injects fuel into the cylinder. The fuel injection valve <NUM> may be of a port injection type that injects fuel into the intake port.

In the cylinder of the internal combustion engine <NUM>, the ignition plug <NUM> that exposes an electrode portion in the cylinder and ignites a combustible air-fuel mixture by a spark is assembled. The knock sensor <NUM> is provided in the cylinder block, and detects the presence or absence of knock by detecting cylinder block vibration caused by combustion pressure vibration generated in the combustion chamber. The crank angle sensor <NUM> is assembled to a crankshaft, and outputs a signal corresponding to a rotation angle of the crankshaft to the ECU <NUM> described later as a signal indicating a rotation speed.

The air-fuel ratio sensor <NUM> is provided on the downstream side of the turbine 3b of the turbocharger <NUM>, and outputs a signal indicating the detected exhaust gas composition, that is, the air-fuel ratio to the ECU <NUM>. The exhaust purification catalyst <NUM> is, for example, a three-way catalyst, is provided on the downstream side of the air-fuel ratio sensor <NUM>, and purifies harmful exhaust gas components such as carbon monoxide, nitrogen compounds, and unburned hydrocarbon in the exhaust gas by a catalytic reaction. In general, platinum and rhodium, or a mixture of platinum and rhodium and palladium are used as the catalyst substance. The oxygen sensor <NUM> is provided on the downstream side of the exhaust purification catalyst <NUM>, and detects the presence or absence of oxygen contained in the exhaust gas after purification by the exhaust purification catalyst <NUM>. Hereinafter, the exhaust purification catalyst <NUM> may be abbreviated as a "catalyst <NUM>".

The turbocharger <NUM> includes the air bypass valve <NUM> and the wastegate valve <NUM>. The air bypass valve <NUM> is disposed on a bypass flow path connecting the upstream and the on the downstream side of the compressor 3a in order to prevent the pressure from excessively increasing from the on the downstream side of the compressor 3a to the on the upstream side of the throttle valve <NUM>. When the throttle valve <NUM> is rapidly closed in the supercharged state, the air bypass valve <NUM> is opened according to the control of the ECU <NUM>, whereby the compressed intake air on the downstream side of the compressor 3a flows back to the on the upstream side of the compressor 3a through the bypass flow path. As a result, by immediately decreasing the supercharging pressure, a phenomenon called surging can be prevented, and damage of the compressor 3a can be appropriately prevented.

The wastegate valve <NUM> is disposed on a bypass flow path connecting upstream and on the downstream side of the turbine 3b. The wastegate valve <NUM> is an electric valve whose valve opening degree can be freely controlled with respect to supercharging pressure by the control of the ECU <NUM>. When the opening degree of the wastegate valve <NUM> is adjusted by the ECU <NUM> based on the supercharging pressure detected by the supercharging pressure sensor <NUM>, a part of the exhaust gas passes through the bypass flow path, so that the work given to the turbine 3b by the exhaust gas can be reduced. As a result, the supercharging pressure can be held at the target pressure.

The EGR pipe <NUM> communicates an exhaust flow path on the downstream side of the exhaust purification catalyst <NUM> and an intake flow path on the upstream side of the compressor 3a, splits the exhaust gas from the on the downstream side of the exhaust purification catalyst <NUM>, and returns the exhaust gas to the on the upstream side of the compressor 3a. The EGR cooler <NUM> provided in the EGR pipe <NUM> cools the exhaust gas. The EGR valve <NUM> is provided on the downstream side of the EGR cooler <NUM> and controls the flow rate of the exhaust gas. The EGR pipe <NUM> is provided with the temperature sensor <NUM> that detects the temperature of the exhaust gas flowing on the upstream side of the EGR valve <NUM> and the differential pressure sensor <NUM> that detects the differential pressure between upstream and on the downstream side of the EGR valve <NUM>.

The ECU <NUM> is an example of an internal combustion engine control device, and controls each component of the engine system and executes various data processing. The above-described various sensors and various actuators are communicably connected to the ECU <NUM>. The ECU <NUM> controls the operations of actuators such as the throttle valve <NUM>, the fuel injection valve <NUM>, the intake valve <NUM>, the exhaust valve <NUM>, and the EGR valve <NUM>. In addition, the ECU <NUM> detects the operating state of the internal combustion engine <NUM> based on signals input from various sensors, and causes the ignition plug <NUM> to ignite at a timing determined according to the operating state.

<FIG> is a block diagram illustrating a hardware configuration example of the ECU <NUM>. The electronic control unit (ECU) <NUM> includes a control unit <NUM>, a storage unit <NUM>, and an input/output interface <NUM> connected to one another via a system bus. The control unit <NUM> includes a central processing unit (CPU) 31a, a read only memory (ROM) 31b, and a random access memory (RAM) 31c. The CPU 31a executes the control program stored in the ROM 31b to implement each function of the ECU <NUM>.

The input/output interface <NUM> is an interface that communicates signals and data with each sensor and each actuator. The ECU <NUM> includes an analog/digital (A/D) converter (not illustrated) that processes input/output signals of each sensor, a driver circuit, and the like. The input/output interface <NUM> may also serve as an A/D converter. Although a central processing unit (CPU) is used as the processor, another processor such as a micro processing unit (MPU) may be used. In addition, the control program may be stored in the storage unit <NUM> which is an auxiliary storage device including a semiconductor memory or the like.

Hereinafter, a control method of the internal combustion engine <NUM> that realizes a fuel-efficient operation by the throttle system, the wastegate system, and the EGR system included in the internal combustion engine <NUM> will be described.

<FIG> is a diagram for explaining a method of controlling the throttle valve <NUM> and the wastegate valve <NUM> and an operation area where EGR (recirculated exhaust gas) is introduced in an operation area defined by the rotation speed and the filling efficiency of the internal combustion engine <NUM>. The filling efficiency is the ratio of the mass of air sucked into the cylinder in one cycle to the mass of air in the standard state corresponding to the cylinder volume. The upper side of <FIG> is a graph when cold EGR is not introduced, and the lower side of <FIG> is a graph when cold EGR is introduced. In each graph of <FIG>, the horizontal axis represents the rotation speed, and the vertical axis represents the filling efficiency.

As illustrated in the upper graph of <FIG>, the operation area of the internal combustion engine <NUM> is roughly divided into a non-supercharging range and a supercharging range. In the non-supercharging range, the filling efficiency is controlled by the throttle valve <NUM>, and the wastegate valve <NUM> is fully opened. In the supercharging range, the throttle valve <NUM> is fully opened, and the supercharging pressure is controlled by the wastegate valve <NUM> to control the filling efficiency. As described above, by switching the means (throttle valve <NUM> and wastegate valve <NUM>) for adjusting torque between the non-supercharging range and the supercharging range, the pump loss generated in the internal combustion engine <NUM> can be reduced, and the fuel-efficient operation can be realized. The thin broken line represents a fresh air equal-flow rate line.

Further, the internal combustion engine <NUM> described in this embodiment is equipped with an EGR system. In a region (a portion surrounded by a thick broken line in the lower graph of <FIG>) from a relatively high load condition of the non-supercharging range to the supercharging region of the internal combustion engine <NUM>, the EGR cooled by the EGR cooler <NUM> is returned to the cylinder, so that the gas sucked into the cylinder is diluted by the EGR as the inert gas, and it is possible to suppress the irregular combustion called knock that is likely to occur under the high load condition. This EGR system is called a cold EGR. Since knocking can be suppressed in this manner, the ignition timing can be appropriately advanced, and a fuel-efficient operation can be realized.

In the engine system of the embodiment described here, the wastegate valve <NUM> is fully opened in the non-supercharging range. However, the present invention is not limited thereto, and a system of fully closing or holding in an intermediate state can also be adopted. As a result, it is possible to increase the rotational speed of the turbocharger <NUM> and improve the acceleration performance while reducing the pump loss reduction effect.

<FIG> is a diagram for explaining a configuration of a post-processing system for purifying exhaust gas of the internal combustion engine <NUM>. As the exhaust purification catalyst <NUM>, a three-way catalyst is used. For the purpose of maintaining the purification efficiency of the three-way catalyst at an optimum point, sensors for detecting exhaust gas compositions are provided upstream and on the downstream side of the three-way catalyst, respectively. Three post-processing systems A, B, and C can be configured according to the type of the sensor described above.

In <FIG>, the upstream side of the catalyst is referred to as "front", and the downstream side of the catalyst is referred to as "rear".

The post-processing system A includes an air-fuel ratio sensor (in the drawing, the front air-fuel ratio sensor) on the upstream side of the catalyst <NUM> and an oxygen sensor (in the drawing, a rear oxygen sensor) on the downstream side of the catalyst. According to this configuration, the air-fuel ratio of the exhaust gas flowing into the three-way catalyst can be measured, and the presence or absence of oxygen contained in the exhaust gas after catalyst purification can be detected.

The post-processing system B includes an oxygen sensor (in the drawing, the front oxygen sensor) on the upstream side of the catalyst <NUM> and an oxygen sensor on the downstream side of the catalyst. According to this configuration, the presence or absence of oxygen contained in the exhaust gas flowing into the three-way catalyst can be detected, and the presence or absence of oxygen contained in the exhaust gas after purification of the catalyst can be detected.

The post-processing system C includes an air-fuel ratio sensor on the upstream side of the catalyst <NUM> and an air-fuel ratio sensor on the downstream side of the catalyst. According to this configuration, the air-fuel ratio of the exhaust gas flowing into the three-way catalyst can be measured, and the air-fuel ratio of the exhaust gas after the catalyst purification can be measured. The post-processing systems A to C described above are selected based on exhaust regulations and cost requirements. The engine system described in this embodiment is based on the post-processing system A.

However, the present invention is not limited to the configuration of the post-processing system A, and it is possible to achieve a similar effect without significantly changing the configuration of the invention by switching the input and output of a control model (see <FIG>) to be described later as described above.

<FIG> is a diagram for explaining the relationship between the equivalent ratio (= theoretical air-fuel ratio/air-fuel ratio) of the exhaust gas and the output of the air-fuel ratio sensor. In <FIG>, the horizontal axis represents the equivalent ratio, and the vertical axis represents the air-fuel ratio sensor output. The air-fuel ratio sensor output tends to decrease as the equivalent ratio increases (as the equivalent ratio becomes richer). The air-fuel ratio sensor is characterized in that the equivalent ratio can be accurately detected in a wide range from the lean state to the rich state for the exhaust gas by acquiring the information indicating the relationship of <FIG> in advance.

<FIG> is a diagram for explaining the relationship between the equivalent ratio of the exhaust gas and the output of the oxygen sensor. In <FIG>, the horizontal axis represents the equivalent ratio, and the vertical axis represents the oxygen sensor output. The oxygen sensor output is represented by an electromotive force associated with a concentration difference between the oxygen concentration contained in the exhaust gas and the oxygen concentration in the air. In the lean condition, almost the minimum electromotive force is shown, and in the rich condition, the maximum electromotive force is shown. Therefore, the output has a characteristic of rapidly changing at the theoretical air-fuel ratio (equivalent ratio <NUM>). A characteristic is that the exhaust gas equivalent ratio can be held in the vicinity of the theoretical air-fuel ratio by capturing the change timing of the oxygen sensor output and feeding back the change timing to the air-fuel ratio control.

<FIG> is a diagram for explaining the tendency of the chemical species concentration of the exhaust gas with respect to the equivalent ratio. The horizontal axis of the upper graph of <FIG> represents the equivalent ratio, the vertical axis represents the concentration [%], the horizontal axis of the lower graph of <FIG> represents the equivalent ratio, and the vertical axis represents the concentration [ppm].

As shown in the upper graph of <FIG>, the combustion gas composition of the hydrocarbon-based fuel shows a tendency that carbon monoxide (CO) and hydrogen (H<NUM>) increase on the rich side and oxygen (O<NUM>) increases on the lean side with the theoretical air-fuel ratio (equivalent ratio <NUM>) as a boundary (upper graph of <FIG>).

On the other hand, as shown in the lower graph of <FIG>, the nitrogen oxide (NOx) shows a maximum value on a slightly lean side of the theoretical air-fuel ratio, and tends to decrease on a lean side and a rich side of the maximum value. Unburned hydrocarbon (HC) is a component that is discharged without reaching combustion, and there is no clear tendency with respect to the equivalent ratio, but when being excessively lean or rich, HC discharged without reaching normal combustion tends to increase.

Further, in a theoretical air-fuel ratio condition in which fuel and air (oxygen) are supplied without excess or deficiency, since a certain amount of CO or NOx is discharged without reaching water (H<NUM>O) or carbon dioxide (CO<NUM>) in the high-temperature combustion gas, it is necessary to appropriately purify the exhaust gas by the post-processing system.

Here, a main reaction process of the three-way catalyst (ceria system) used in the post-processing system will be described.

The reaction process of the three-way catalyst mainly includes an oxidation reaction, an NOx reduction reaction, and an oxygen storage/release reaction. In the oxidation reaction, CO, H<NUM>, and HC generated under rich or high temperature conditions react with oxygen to generate harmless CO<NUM> and H<NUM>O. The unburned hydrocarbon (HC) contains components such as methane, propane, ethylene, and butane, and the reactions proceed at different rates (Reaction Expressions (<NUM>) to (<NUM>)). In the NOx reduction reaction, CO and NO are mainly reacted to produce harmless CO<NUM> and N<NUM> (Reaction Expression (<NUM>)). In the oxygen storage/release reaction, storage/release of oxygen and each oxidation/reduction reaction of HC, CO, and NO proceed via Ce (cerium) which is a catalyst material (Reaction Expressions (<NUM>) to (<NUM>)). Each reaction expression is expressed in the form of an elementary reaction expression.

H<NUM> + O<NUM> ⇒ H<NUM><NUM>.

CnHm + O<NUM> ⇒ CO<NUM> + H<NUM>O.

CO + NO ⇒ CO<NUM> + N<NUM>.

CeO<NUM> + CO ⇒ Ce<NUM>O<NUM> + CO<NUM>.

CnHm + CeO<NUM> ⇒ CeO<NUM>O<NUM> + CO + H<NUM>O.

Ce<NUM>O<NUM> + O<NUM> ⇒ CeO<NUM>.

Ce<NUM>O<NUM> + NO ⇒ CeO<NUM> + N<NUM>.

As described above, harmless CO<NUM> and H<NUM>O are generated by the reaction of cerium dioxide (CeO<NUM>) with CO and HC, and harmless N<NUM> is generated by the reaction of dicerium trioxide (Ce<NUM>O<NUM>) with NO. At this time, the oxygen storage ratio of the three-way catalyst is defined by the balance between CeO<NUM> and Ce<NUM>O<NUM> generated simultaneously. That is, when all of Ce<NUM>O<NUM> in the catalyst becomes CeO<NUM>, a reaction with NO cannot be performed, and NO cannot be purified. As described above, in order to appropriately maintain the purification efficiency of the three-way catalyst, it is necessary to maintain the balance between CeO<NUM> and Ce<NUM>O<NUM>, that is, the oxygen storage ratio at a predetermined value. All the above-described reaction processes strongly depend on the catalyst temperature, and it is necessary to appropriately manage the catalyst temperature so that the temperature becomes equal to or higher than the activation temperature early after the engine is started.

In the engine system shown in this embodiment, the ceria-based three-way catalyst is used, but the present invention is not limited thereto. Even in a catalyst using another material exhibiting a similar effect, the same effect can be exhibited without changing the configuration of the invention by adjusting the constant of the control model. In addition, in the catalytic reaction, an aqueous gas shift reaction or the like may be used in addition to the reaction mechanism described above, but these reaction mechanisms can also be handled by adjusting the constant of the control model.

<FIG> is a diagram for explaining the tendency of the purification efficiency of the three-way catalyst with respect to the exhaust gas equivalent ratio at the catalyst activation temperature or higher. In <FIG>, the horizontal axis represents the equivalent ratio, and the vertical axis represents the catalyst purification efficiency [%].

The purification efficiency characteristic of the three-way catalyst changes with a theoretical air-fuel ratio (equivalent ratio <NUM>) as a boundary. Under the lean condition, while the purification efficiency of CO and HC is maintained at approximately <NUM>% or more, the purification efficiency of NOx decreases as the equivalent ratio decreases. Under the rich condition, the purification efficiency of HC and CO tends to decrease as the equivalent ratio decreases. In the vicinity of the theoretical air-fuel ratio, the purification efficiency of <NUM>% or more can be achieved for any of NOx, HC, and CO, and this point is referred to as a three-way point. In the three-way catalyst, control is performed to keep the purification efficiency high by keeping the equivalent ratio near the theoretical air-fuel ratio (control target) including the three-way point.

<FIG> is a diagram for explaining the output behavior of the oxygen sensor <NUM> installed on the downstream side of the catalyst <NUM> when the air-fuel ratio is varied stepwise over time to the lean side or the rich side with the equivalent ratio <NUM> as a center. In each graph of <FIG>, the horizontal axis represents time, and the vertical axis represents an equivalent ratio (upper graph) and a rear oxygen sensor output (lower graph).

Even when the equivalent ratio is set to the theoretical air-fuel ratio, a very small amount of oxygen is discharged on the downstream side of the catalyst, so that the intermediate state is maintained. Then, when the equivalent ratio is changed stepwise over time to the lean side, the rear oxygen sensor output changes to the minimum value side after the delay period d1. On the other hand, when the air-fuel ratio is changed from the lean side to the rich side, the rear oxygen sensor output exhibits hysteresis that changes with a larger delay period d2.

As described above, the characteristic of the oxygen sensor is that the delay time tends to be different between the change from lean to rich and the change from rich to lean. This is because the reaction rates of CeO<NUM> and Ce<NUM>O<NUM> described above are different. Since the reaction rate depends on the catalyst temperature, the hysteresis described above also changes depending on the catalyst temperature.

<FIG> is a diagram for explaining temporal changes in the rear oxygen sensor output and the catalyst downstream NOx concentration when the firing operation is performed again in the theoretical air-fuel ratio state after the motoring operation (fuel cut) period from the state where the control is performed at the theoretical air-fuel ratio. In each graph of <FIG>, the horizontal axis represents time, and the vertical axis represents an equivalent ratio (upper graph), a rear oxygen sensor output (middle graph), and a catalyst downstream NOx concentration (lower graph). In each graph, the solid line indicates the behavior during the stoichiometric control, and the broken line indicates the behavior after the rich correction.

When the fuel cut-off is performed, O<NUM> sucked into the cylinder of the internal combustion engine <NUM> increases. When the firing operation using the theoretical air-fuel ratio is started again after the fuel cut, the rear oxygen sensor output increases with a delay from the minimum value as indicated by the solid line as described in <FIG>. Then, in a delay period until the rear oxygen sensor output returns to the maximum value, the catalyst downstream NOx concentration increases in a spike manner, and a large amount of NOx is discharged.

On the other hand, when starting the firing operation again after the fuel cut, if the theoretical air-fuel ratio control is performed after performing the rich correction to change the equivalent ratio indicated by the broken line to the rich side, NOx emission on the downstream side of the catalyst can be prevented. In other words, NOx emission can be prevented by sucking rich gas (large amount of fuel and small amount of O<NUM>) into the cylinder of the internal combustion engine <NUM> before the rear oxygen sensor output changes from the minimum value.

As understood from <FIG>, since the oxygen sensor <NUM> detects the oxygen state of the exhaust gas on the downstream side of the catalyst <NUM>, the internal state of the catalyst has already changed to the maximum or minimum state of the oxygen storage state at the time when the oxygen sensor <NUM> reacts. That is, in the conventional method of performing feedback control after the reaction of the rear oxygen sensor, since the control timing is too late for the catalyst <NUM>, emission deterioration cannot be appropriately prevented. Therefore, in the air-fuel ratio control of the internal combustion engine <NUM>, it can be said that it is necessary to perform rich correction control for an appropriate period in consideration of a catalyst state that cannot be directly observed from the outside.

<FIG> is a diagram for explaining the relationship between the deterioration degree of the three-way catalyst and the oxygen storage capacity. In <FIG>, the horizontal axis represents the degree of catalyst deterioration, and the vertical axis represents the oxygen storage capacity. The catalyst deterioration refers to a state in which the catalytic action is reduced due to a thermal influence or an influence of poisoning by sulfur contained in the fuel. As shown in <FIG>, the oxygen storage capacity tends to decrease as the catalyst deterioration progresses. Hereinafter, the influence of the change in oxygen storage capacity on the purification action of the three-way catalyst will be described.

<FIG> is a diagram for explaining a comparison result of the output behavior of the oxygen sensor <NUM> installed on the downstream side of the catalyst when the air-fuel ratio is varied stepwise over time to a lean side or a rich side with an equivalent ratio of <NUM> as a center in the new catalyst and the deterioration catalyst. In each graph of <FIG>, the horizontal axis represents time, and the vertical axis represents an equivalent ratio (upper graph), an oxygen storage ratio (middle graph), and a rear oxygen sensor output (lower graph). In each graph, a solid line indicates the behavior of the deterioration catalyst, and a broken line indicates the behavior of the new catalyst.

In the deterioration catalyst, the delay of the rear oxygen sensor output behavior with respect to the change between the lean and rich air-fuel ratios is reduced as compared with the new catalyst (broken line portion in the lower graph of <FIG>). This can be explained by the temporal transition of the oxygen storage ratio of the catalyst. That is, since the oxygen storage capacity of the catalyst decreases due to deterioration, and the oxygen storage ratio reaches the maximum value or the minimum value more quickly (solid line portion in the lower graph of <FIG>), the oxygen release behavior toward the on the downstream side of the catalyst is accelerated, and the delay in the rear oxygen sensor output behavior decreases. Therefore, the rich correction period after recovery from the fuel cut described in <FIG> needs to be set in consideration of the deterioration state of the catalyst.

<FIG> is a diagram for explaining the relationship between the oxygen storage ratio and the NOx purification efficiency. In <FIG>, the horizontal axis represents the oxygen storage ratio, and the vertical axis represents the NOx purification efficiency. When the oxygen storage amount is <NUM>% relative to the oxygen storage capacity, the oxygen storage ratio is <NUM>.

When the oxygen storage ratio exceeds a predetermined value, the NOx purification efficiency is significantly deteriorated. This is because, as described above, Ce<NUM>O<NUM> in the catalyst is important for NOx purification, but when all of Ce<NUM>O<NUM> reacts to change to CeO<NUM>, a reaction between Ce<NUM>O<NUM> and NO does not occur, and NO cannot be purified. Therefore, as described in <FIG>, in order to keep the catalyst purification efficiency high, it is necessary not only to keep the exhaust gas air-fuel ratio at the catalyst inlet at the three-way point but also to appropriately correct and control the exhaust gas air-fuel ratio at the catalyst inlet so that the oxygen storage ratio falls within a predetermined control range.

<FIG> is a diagram for explaining a block diagram illustrating a function of a control model for correcting and controlling the air-fuel ratio in consideration of the catalyst state. This control model is realized by the ECU <NUM>. The ECU <NUM> includes a catalyst boundary condition calculation unit <NUM>, a catalyst temperature calculation unit <NUM>, an oxygen storage ratio calculation unit <NUM>, a statistical model calculation unit <NUM>, a system identification unit <NUM>, an air-fuel ratio correction amount calculation unit <NUM>, a target air-fuel ratio calculation unit <NUM>, and a deterioration degree determination unit <NUM>.

The catalyst boundary condition calculation unit <NUM> calculates a catalyst boundary condition based on the operating state of the internal combustion engine <NUM> and the like. The catalyst boundary conditions include a catalyst upstream exhaust gas temperature (hereinafter referred to as "exhaust gas temperature"), a catalyst upstream exhaust gas composition (exhaust gas air-fuel ratio), a catalyst inflow exhaust gas flow rate (hereinafter referred to as" exhaust gas flow rate "), atmospheric conditions (outdoor temperature and atmospheric pressure), and a vehicle speed. The calculation result of the catalyst boundary condition is input to the catalyst temperature calculation unit <NUM> and the oxygen storage ratio calculation unit <NUM>.

The catalyst temperature calculation unit <NUM> calculates the catalyst temperature based on the catalyst boundary condition calculated by the catalyst boundary condition calculation unit <NUM>. By providing the catalyst temperature calculation unit <NUM>, the influence of the catalyst temperature change important in the catalytic reaction can be reflected in the air-fuel ratio control.

The oxygen storage ratio calculation unit <NUM> calculates the oxygen storage ratio based on the catalyst boundary condition, the catalyst temperature, the detection value of the air-fuel ratio sensor <NUM> (front air-fuel ratio sensor) on the upstream side of the catalyst, and the oxygen storage capacity defined by the catalyst deterioration degree. At this time, the relationship shown in <FIG> in which the oxygen storage capacity decreases with respect to the progress of catalyst deterioration is considered. Although it is necessary to consider the above-described reaction mechanism for the calculation of the oxygen storage ratio, on-board calculation as a control model is impossible, and thus a simple reaction mechanism is substituted. A specific model configuration will be described later. The calculation result of the oxygen storage ratio is input to the statistical model calculation unit <NUM> and the system identification unit <NUM>. The model handled by the catalyst temperature calculation unit <NUM> and the oxygen storage ratio calculation unit <NUM> is a catalytic reaction model and can be said to be a physical model.

As described above, in the internal combustion engine control device (ECU <NUM>) of this embodiment, the catalytic reaction model (oxygen storage ratio calculation unit <NUM>) receives at least the detection value of the first exhaust gas sensor (for example, the air-fuel ratio sensor <NUM>), the catalyst temperature, and the exhaust gas flow rate on the upstream side of the catalyst <NUM> in the exhaust pipe as inputs.

The catalytic reaction model (Reaction Expressions (<NUM>) to (<NUM>)) of this embodiment is defined by at least a reaction rate between oxygen on the upstream side of the catalyst <NUM>, carbon monoxide on the upstream side of the catalyst <NUM>, and a metal (platinum, rhodium, etc.) carried in the catalyst <NUM>.

As described above, by accurately calculating the temporal change of the oxygen storage ratio in consideration of the influence of the exhaust gas flow rate, the catalyst temperature, and the catalyst upstream exhaust gas air-fuel ratio, the air-fuel ratio correction control for maintaining the catalyst purification efficiency high in this embodiment can be realized.

The statistical model calculation unit <NUM> calculates a statistical model <NUM> having the calculation result of the oxygen storage ratio as an input and the detection value of the oxygen sensor <NUM> (rear oxygen sensor) on the downstream side of the catalyst <NUM> as an output. A tuning parameter is set in the statistical model <NUM>, and is always updated to the latest state by a system identification algorithm described later. With the statistical model <NUM>, the behavior of the future rear oxygen sensor detection value can be predictively calculated.

When an identification permission flag is in the ON state, the system identification unit <NUM> updates the tuning parameter set to the statistical model <NUM> of the statistical model calculation unit <NUM> by the system identification algorithm based on the calculated oxygen storage ratio and the rear oxygen sensor detection value. With the system identification algorithm, the catalyst characteristics and the sensor characteristics that change due to deterioration can be sequentially reflected in the tuning parameter, and the calculation accuracy of the statistical model <NUM> can be kept high. The processing result (tuning parameter and the like) of the system identification unit <NUM> is input to the deterioration degree determination unit <NUM>.

The air-fuel ratio correction amount calculation unit <NUM> calculates the air-fuel ratio correction amount based on the calculation result of the statistical model <NUM>. By predictively calculating the future behavior of the rear oxygen sensor detection value by the statistical model <NUM>, the air-fuel ratio correction amount can be calculated before the rear oxygen sensor actually changes. Therefore, as compared with the feedback control based on the conventional rear oxygen sensor detection value, the oxygen storage ratio in the catalyst can be more suitably maintained within a predetermined range.

The target air-fuel ratio calculation unit <NUM> corrects the air-fuel ratio determined based on the operating state with the calculation result of the air-fuel ratio correction amount calculated by the air-fuel ratio correction amount calculation unit <NUM>, and calculates the target air-fuel ratio. The control unit <NUM> controls the air-fuel ratio based on the target air-fuel ratio calculated by the target air-fuel ratio calculation unit <NUM>, whereby the catalyst purification efficiency can be held high to improve emission performance, and deterioration of emission performance can be prevented by adapting to characteristic changes such as catalyst deterioration. In <FIG>, the target air-fuel ratio calculation unit <NUM> acquires the operating state and determines the air-fuel ratio, but a calculation unit different from the target air-fuel ratio calculation unit <NUM> may determine the air-fuel ratio based on the operating state.

The deterioration degree determination unit <NUM> determines the deterioration degree of the catalyst <NUM> (hereinafter, referred to as "catalyst deterioration degree") based on the processing result of the system identification unit <NUM>. The calculation result of the catalyst deterioration degree is input to the oxygen storage ratio calculation unit <NUM>.

<FIG> is a diagram for explaining calculation contents of the catalyst boundary condition in the catalyst boundary condition calculation unit <NUM> illustrated in <FIG>. The catalyst boundary condition calculation unit <NUM> includes an exhaust gas flow rate calculation unit <NUM> and an exhaust gas temperature calculation unit <NUM>.

The exhaust gas flow rate calculation unit <NUM> calculates the exhaust gas flow rate based on the rotational speed, the load (rotational torque), and the state flag of the internal combustion engine <NUM>. Similarly, the exhaust gas temperature calculation unit <NUM> calculates the exhaust gas temperature based on the rotational speed, the load, and the state flag of the internal combustion engine <NUM>. Note that a temperature sensor may be provided upstream or on the downstream side of the catalyst <NUM>, and the exhaust gas temperature may be detected by the temperature sensor.

The state flag is information for determining a control state such as fuel cut or ignition retard. In the cold start mode, catalyst temperature rise control is performed by ignition retard, intake air increase, and the like. In an internal combustion engine having a sailing stop control function or an internal combustion engine for a hybrid vehicle, a motoring state due to fuel cut frequently occurs, and in this case, air corresponding to atmospheric temperature passes through the catalyst. As described above, since the catalyst upstream temperature is affected by various control states, the catalyst state can be calculated with high accuracy by considering these state transitions.

<FIG> is a diagram for explaining calculation contents of the catalyst temperature in the catalyst temperature calculation unit <NUM> shown in <FIG>. The catalyst temperature calculation unit <NUM> includes an exhaust gas/catalyst heat transfer amount calculation unit <NUM>, a catalyst/atmosphere heat transfer amount calculation unit <NUM>, a catalyst reaction generation heat amount calculation unit <NUM>, a catalyst heat balance calculation unit <NUM>, and a catalyst temperature calculation unit <NUM>.

The exhaust gas/catalyst heat transfer amount calculation unit <NUM> calculates the heat transfer amount between the exhaust gas and the catalyst <NUM> based on the exhaust gas flow rate and the exhaust gas temperature calculated by the catalyst boundary condition calculation unit <NUM> and the past catalyst temperature value. The heat transfer amount can be accurately calculated by considering the exhaust gas flow rate and the exhaust gas temperature, which are main factors of the turbulent heat transfer.

The catalyst/atmosphere heat transfer amount calculation unit <NUM> calculates the heat transfer amount between the catalyst <NUM> and the atmosphere based on the vehicle speed, the outdoor temperature and the atmospheric pressure, and the past catalyst temperature value. Regarding heat transfer around the catalyst <NUM>, a forced convection phenomenon is dominant when the vehicle speed is large, and a natural convection phenomenon is dominant when the vehicle speed is small or in a stopped state. By considering such a change in the heat transfer mechanism based on the vehicle speed information, the heat transfer amount around the catalyst <NUM> can be accurately calculated.

The catalyst reaction generation heat amount calculation unit <NUM> calculates the heat amount generated by the catalytic reaction based on the exhaust gas flow rate calculated by the catalyst boundary condition calculation unit <NUM>, the past catalyst temperature value, and the upstream exhaust gas concentration calculated by an exhaust gas concentration calculation unit <NUM>. The inside of the catalyst <NUM> is affected by heat generation due to the oxidation reaction of the unburned gas, endotherm due to the NOx reduction reaction, and the like, and these reaction rates strongly depend on the catalyst temperature. The catalyst reaction generation heat amount calculation unit <NUM> can accurately calculate a catalytic reaction generation heat amount in the catalyst <NUM> in consideration of these.

The catalyst heat balance calculation unit <NUM> calculates a balance of each heat transfer amount calculated by the exhaust gas/catalyst heat transfer amount calculation unit <NUM> and the catalyst/atmosphere heat transfer amount calculation unit <NUM> and a heat amount calculated by the catalyst reaction generation heat amount calculation unit <NUM>.

The catalyst temperature calculation unit <NUM> calculates the temporal change of the catalyst temperature in consideration of catalyst specifications (volume, aperture ratio, surface area, mass, specific heat, and the like) based on the calculation result of the catalyst heat balance by the catalyst heat balance calculation unit <NUM>. The calculation result of the catalyst temperature is input to the exhaust gas/catalyst heat transfer amount calculation unit <NUM> and the catalyst/atmosphere heat transfer amount calculation unit <NUM>. With such a configuration, even in an engine system in which the starting and stopping of the internal combustion engine <NUM> are repeated, the catalyst temperature important for estimating the catalyst state can be calculated with high accuracy.

Here, a configuration of an oxygen storage ratio calculation model that is assumed to be operated online (in real time) in the ECU <NUM> will be described.

The reaction of the three-way catalyst is described by the eight reaction mechanisms of Reaction Expressions (<NUM>) to (<NUM>) described above. However, it is not realistic to perform the online calculation of the reaction expression from the viewpoint of calculation load, and thus, in the engine system of this embodiment, a simplified reaction model assuming control implementation is used. Among exhaust gas compositions, HC and NOx concentrations are sufficiently smaller than those of other chemical species, and thus only CO and O<NUM> are handled in the simplified reaction model. Then, two Reaction Expressions (<NUM>) and (<NUM>) involving CeO<NUM> and Ce<NUM>O<NUM> defining the oxygen storage ratio are simultaneously solved. The oxygen storage ratio is defined by the ratio of the number of moles of CeO<NUM> and Ce<NUM>O<NUM> as shown in Reaction Expression (<NUM>). Each reaction expression is expressed in the form of an elementary reaction expression.

Ce<NUM>O<NUM> + O2 ⇒ CeO<NUM>.

The reaction rates R1 and R2 of CO and O<NUM> are calculated by the following Expressions (<NUM>) and (<NUM>), respectively. Expressions (<NUM>) and (<NUM>) are described in Arrhenius type, and the catalyst temperature, the molar concentration of CO and O<NUM>, the oxygen storage ratio Ψ, and the oxygen storage capacity Ψcap are considered. The value of the oxygen storage capacity Ψcap may be changed depending on the catalyst deterioration degree. For example, different values are set depending on whether the catalyst is new, not so new, or old. Here, A is a frequency factor, E is activation energy, R is a general gas constant, and Teat is a catalyst temperature. Since the expression is an exponential function, the calculation load can be reduced by replacing the expression with the table calculation regarding the catalyst temperature Teat. Data set in the table is set in advance on the basis of an experimental result.

As long as the performance such as the processing speed of the ECU <NUM> is improved, the oxygen storage ratio is not limited to Reaction Expressions (<NUM>) and (<NUM>) for CO and O<NUM>, and the oxygen storage ratio may be calculated in consideration of other reaction formulas. In this case, it is expected that the calculation result of the oxygen storage ratio approaches an actual value.

<FIG> is a diagram for explaining the calculation content of the oxygen storage ratio in the oxygen storage ratio calculation unit <NUM> shown in <FIG>. The oxygen storage ratio calculation unit <NUM> includes an exhaust gas concentration calculation unit <NUM>, an O<NUM> reaction rate/concentration calculation unit <NUM>, a CO reaction rate/concentration calculation unit <NUM>, and an oxygen storage ratio transition calculation unit <NUM>.

The exhaust gas concentration calculation unit <NUM> calculates the concentration of the exhaust gas chemical species (upstream exhaust gas concentration) based on the detection value of the air-fuel ratio sensor <NUM> (front air-fuel ratio sensor) provided on the upstream side of the catalyst <NUM>.

As the upstream exhaust gas concentration, for example, an upstream O<NUM> concentration and an upstream CO concentration are calculated. The catalyst upstream exhaust gas composition can be accurately calculated by performing calculation based on the measured behavior of the actual exhaust gas and the exhaust gas composition information obtained in advance at the development stage. The calculated upstream exhaust gas concentration is output to the catalyst temperature calculation unit <NUM>. Note that the exhaust gas concentration calculation unit <NUM> only needs to have the function as the ECU <NUM>, and is not limited to the inside of the oxygen storage ratio calculation unit <NUM>.

The O<NUM> reaction rate/concentration calculation unit <NUM> calculates the reaction rate and the concentration of O<NUM> based on the exhaust gas flow rate calculated by the catalyst boundary condition calculation unit <NUM>, the upstream exhaust gas concentration (upstream O<NUM> concentration, upstream CO concentration), the catalyst temperature calculated by the catalyst temperature calculation unit <NUM>, the past oxygen storage ratio value, and the oxygen storage capacity calculated by the oxygen storage ratio calculation unit <NUM>.

Similarly, the CO reaction rate/concentration calculation unit <NUM> calculates the reaction rate and the concentration of CO based on the exhaust gas flow rate, the upstream exhaust gas concentration (upstream O<NUM> concentration, upstream CO concentration), the catalyst temperature, the past oxygen storage ratio value, and the oxygen storage capacity. As described later, the O<NUM> reaction rate/concentration calculation unit <NUM> and the CO reaction rate/concentration calculation unit <NUM> are configured to be able to set the oxygen storage capacity according to the catalyst deterioration degree.

Further, the oxygen storage ratio transition calculation unit <NUM> calculates the temporal transition of the oxygen storage ratio based on the reaction rate and the concentration of O<NUM> calculated by the O<NUM> reaction rate/concentration calculation unit <NUM> and the reaction rate and the concentration of CO calculated by the CO reaction rate/concentration calculation unit <NUM>.

As described above, the oxygen storage ratio calculation unit (oxygen storage ratio calculation unit <NUM>) of this embodiment includes the oxygen concentration calculation unit (O<NUM> reaction rate/concentration calculation unit <NUM>) that calculates the oxygen concentration on the upstream side of the catalyst based on the detection value of the first exhaust gas sensor on the upstream side of the catalyst, and the carbon monoxide concentration calculation unit (CO reaction rate/concentration calculation unit <NUM>) that calculates the carbon monoxide concentration on the upstream side of the catalyst based on the detection value of the first exhaust gas sensor on the upstream side of the catalyst.

With such a configuration, the temporal transition of the oxygen storage ratio having the hysteresis can be accurately calculated in consideration of the air-fuel ratio change on the upstream side of the catalyst <NUM> and the catalyst deterioration state.

In this embodiment, the catalyst temperature is included in the input parameter of the catalytic reaction model (oxygen storage ratio calculation unit <NUM>). As described, as described for the reaction process of the three-way catalyst using Expressions (<NUM>) to (<NUM>), each reaction process strongly depends on the catalyst temperature. Therefore, the oxygen storage ratio can be accurately calculated by using the catalyst temperature as the input parameter of the catalytic reaction model.

<FIG> is a diagram for explaining a model configuration for improving calculation accuracy of the catalyst temperature and the oxygen storage ratio in the catalyst temperature calculation unit <NUM> and the oxygen storage ratio calculation unit <NUM> illustrated in <FIG>. A catalyst temperature/oxygen storage ratio calculation unit <NUM> illustrated in <FIG> includes exhaust gas temperature/catalyst temperature calculation units <NUM>, <NUM>, and <NUM>, O<NUM>/CO concentration/oxygen storage ratio calculation units <NUM>, <NUM>, and <NUM>, and an average oxygen storage ratio calculation unit <NUM>.

The exhaust gas temperature/catalyst temperature calculation units <NUM>, <NUM>, and <NUM> are models obtained by spatially dividing the catalyst <NUM>. That is, the exhaust gas temperature/catalyst temperature calculation units <NUM>, <NUM>, and <NUM> are division models when the catalyst <NUM> is divided into three from the catalyst inlet to the catalyst outlet. The internal functions of the exhaust gas temperature/catalyst temperature calculation units <NUM>, <NUM>, and <NUM> are the same, and configured by the catalyst temperature calculation unit <NUM> illustrated in <FIG>.

Similarly, the O<NUM>/CO concentration/oxygen storage ratio calculation units <NUM>, <NUM>, and <NUM> are also division models when the catalyst <NUM> is divided into three from the catalyst inlet to the catalyst outlet. The internal functions of the O<NUM>/CO concentration/oxygen storage ratio calculation units <NUM>, <NUM>, and <NUM> are the same, and are configured by the oxygen storage ratio calculation unit <NUM> illustrated in <FIG>.

The exhaust gas temperature/catalyst temperature calculation unit <NUM> at the catalyst inlet outputs the calculation result of the exhaust gas temperature to the exhaust gas temperature/catalyst temperature calculation unit <NUM> in the center of the catalyst and outputs the calculation result of the catalyst temperature to the O<NUM>/CO concentration/oxygen storage ratio calculation unit <NUM> at the catalyst inlet. The O<NUM>/CO concentration/oxygen storage ratio calculation unit <NUM> outputs the calculation result of the O<NUM>/CO concentration to the exhaust gas temperature/catalyst temperature calculation unit <NUM> and the O<NUM>/CO concentration/oxygen storage ratio calculation unit <NUM> in the center of the catalyst, and outputs the calculation result of the oxygen storage ratio to the average oxygen storage ratio calculation unit <NUM>.

In addition, the exhaust gas temperature/catalyst temperature calculation unit <NUM> in the center of the catalyst outputs the calculation result of the exhaust gas temperature to the exhaust gas temperature/catalyst temperature calculation unit <NUM> at the catalyst outlet and outputs the calculation result of the catalyst temperature to the O<NUM>/CO concentration/oxygen storage ratio calculation unit <NUM> in the center of the catalyst. The O<NUM>/CO concentration/oxygen storage ratio calculation unit <NUM> outputs the calculation result of the O<NUM>/CO concentration to the exhaust gas temperature/catalyst temperature calculation unit <NUM> and the O<NUM>/CO concentration/oxygen storage ratio calculation unit <NUM> at the catalyst outlet, and outputs the calculation result of the oxygen storage ratio to the average oxygen storage ratio calculation unit <NUM>.

Further, the exhaust gas temperature/catalyst temperature calculation unit <NUM> at the catalyst outlet calculates the exhaust gas temperature (downstream exhaust gas temperature), and outputs the calculation result of the catalyst temperature to the O<NUM>/CO concentration/oxygen storage ratio calculation unit <NUM> at the catalyst outlet. The O<NUM>/CO concentration/oxygen storage ratio calculation unit <NUM> outputs the calculation result of the O<NUM>/CO concentration (downstream O<NUM>/CO concentration) to the exhaust gas temperature/catalyst temperature calculation unit <NUM>, and outputs the calculation result of the oxygen storage ratio to the average oxygen storage ratio calculation unit <NUM>.

The average oxygen storage ratio calculation unit <NUM> calculates an average oxygen storage ratio of the entire catalyst <NUM> from the calculation results of the oxygen storage ratios obtained from the O<NUM>/CO concentration/oxygen storage ratio calculation units <NUM>, <NUM>, and <NUM>.

As described above, the inlet to the outlet of the catalyst <NUM> are spatially divided into a plurality of spaces, and the catalyst temperature and the oxygen storage ratio are calculated for each space, so that it is possible to consider the variation of each value for each space inside the catalyst <NUM>. With such a model configuration, it is possible to consider the time delay behavior and the hysteresis of the catalyst <NUM> and to improve the calculation accuracy. The division number can be arbitrarily determined in consideration of a balance between calculation load and accuracy.

<FIG> is a diagram for explaining hysteresis of oxygen sensor characteristics. In <FIG>, the horizontal axis represents the equivalent ratio, and the vertical axis represents the oxygen sensor output.

The static characteristics of the oxygen sensor are as described in <FIG>. On the other hand, a catalyst material is also used for the oxygen sensor, and has hysteresis due to detection delay. That is, when the gas rapidly changes from the lean state to the rich state, the equivalent ratio (broken line) corresponding to the change timing of the oxygen sensor output shifts to the rich side. When the gas rapidly changes from the rich state to the lean state, the equivalent ratio (one-dot chain line) corresponding to the change timing of the oxygen sensor output shifts to the lean side. Further, the behavior described above is affected by characteristic degradation of the material constituting the oxygen sensor.

Therefore, in the control model of this embodiment, it is desirable to consider not only the deterioration of the catalyst <NUM> but also the change in the dynamic characteristics of the entire post-processing system including the deterioration of the oxygen sensor <NUM> on the downstream side of the catalyst <NUM>.

Next, the statistical model <NUM> in the statistical model calculation unit <NUM> illustrated in <FIG> will be described with reference to <FIG>. In the engine system of this embodiment, the deterioration of the catalyst <NUM> and the oxygen sensor <NUM> is treated as a time-varying system, and a method of on-board approximation by sequential system identification of a linear transient model is adopted.

<FIG> illustrates a configuration of a nonlinear linear transient model adopted in the statistical model <NUM>. The statistical model <NUM>, which is a nonlinear linear transient model, includes three elements, and includes an input nonlinear element <NUM>, a linear transient model <NUM>, and an output nonlinear element <NUM>. An oxygen storage ratio is set as an input of the statistical model <NUM>, and a rear oxygen sensor detection value (estimated value) is set as an output of the statistical model <NUM>.

A function representing a characteristic of the input nonlinear element <NUM> is set as f(u), the oxygen storage ratio u(t) is set as an input, and v(t) is set as an output. When a term considering the hysteresis characteristic of the rear oxygen sensor is set to the function f (u) of the input nonlinear element <NUM>, a model equation of the output v is expressed by Expression (<NUM>) as an example. <NUM>] <MAT>.

As described above, in the input element (input nonlinear element <NUM>) of the statistical model <NUM> of this embodiment, the hysteresis term reflecting the characteristic in which the response speed of the output is different between when the input increases and when the input decreases is set. With this configuration, the hysteresis (see <FIG>) of the exhaust gas sensor (for example, the oxygen sensor <NUM>) on the downstream side of the catalyst <NUM> can be approximated. As described in <FIG>, the hysteresis of the exhaust gas sensor on the downstream side of the catalyst is affected by, for example, characteristic deterioration of the material constituting the oxygen sensor <NUM>.

In addition, the hysteresis term of this embodiment includes a term (β1·y·|u|+β2·u) defined by a first input parameter (u) that takes both a positive value and a negative value and a second input parameter (y·|u|) defined by the product of the output parameter (y) that takes a positive value and the absolute value (|u|) of the first input parameter. In this case, the hysteresis of the oxygen sensor <NUM> can be accurately reflected in the input element of the statistical model <NUM>.

In addition, a function representing a characteristic of the output nonlinear element <NUM> is g(x), an input from the linear transient model <NUM> is x(t), and an output is a rear oxygen sensor detection value y(t). When a sigmoid function for approximating the rear oxygen sensor static characteristic (see <FIG>) is set to the function g(x) of the output nonlinear element <NUM>, a model equation of the output y is expressed by Expression (<NUM>) as an example. <NUM>] <MAT>.

As described above, the output element (output nonlinear element <NUM>) of the statistical model <NUM> of this embodiment is set with a sigmoid function. With this configuration, the static characteristics of the exhaust gas sensor (for example, the oxygen sensor <NUM>) on the downstream side of the catalyst <NUM> can be approximated.

Note that the sigmoid function used in the output nonlinear element <NUM> can be replaced with table calculation in order to reduce the calculation load. In addition, the function g(x) of the output nonlinear element <NUM>
is not limited to the sigmoid function as long as the output is large when the equivalent ratio is in the rich state and the output is small when the equivalent ratio is in the lean state.

The model equation of the linear transient model <NUM> is expressed by Expression (<NUM>) as an example. The left side of Expression (<NUM>) represents the output side, and the right side represents the input side. [k] corresponds to a future, and [k-<NUM>] corresponds to the present (strictly, the past one step before). <NUM>] <MAT>.

Here, Expression (<NUM>) representing the linear transient model <NUM> can be transformed as Expression (<NUM>). θ on the right side in Expression (<NUM>) is a transposed matrix of a matrix of <NUM> row and <NUM> columns represented by Expression (<NUM>), and φ is a matrix of <NUM> row and <NUM> columns represented by Expression (<NUM>). a1, b1, and b2 in Expression (<NUM>) are arbitrary coefficients (examples of tuning parameters).

An auto-regressive with eXogenous (ARX) model generally used in control is applied to the linear transient model <NUM>. The order of the ARX model can be selected by a balance between accuracy and calculation load. A tuning parameter is set to the linear transient model <NUM>, and the tuning parameter is updated online. As a result, the latest system state can be always reflected in the statistical model <NUM>.

Since the statistical model <NUM> includes such components, it is possible to consider dynamic characteristics of the entire post-processing system and temporal changes thereof.

Note that the configuration of the statistical model used in the engine system of this embodiment is not limited thereto. For example, since the catalyst temperature and the exhaust gas flow rate are influencing factors of the catalytic reaction rate, it is possible to improve the accuracy of the statistical model by adding influencing factors such as the catalyst temperature and the exhaust gas flow rate to the input variables. In addition, the oxygen concentration on the downstream side of the catalyst can be set instead of the oxygen storage ratio. In the post-processing system C (<FIG>) including the air-fuel ratio sensor on the downstream side of the catalyst, the output nonlinear element <NUM> of the statistical model can be omitted, and the air-fuel ratio sensor detection value can be directly set.

<FIG> is a flowchart illustrating an example of a system identification algorithm in the system identification unit <NUM> illustrated in <FIG>. In this embodiment, a successive least squares approximation algorithm with a variable forgetting element is adopted for the successive identification of the tuning parameters of the statistical model <NUM>. The forgetting element is to exponentially forget the past history according to the age. By adopting the forgetting element, the influence of the latest state of the engine system can be appropriately considered for the tuning parameter.

First, the system identification unit <NUM> determines whether the statistical model <NUM> is in a state in which the tuning parameter can be identified (S2101). The control unit <NUM> of the ECU <NUM> determines whether the identification is possible based on the operating state of the internal combustion engine <NUM>. The system identification unit <NUM> refers to an identification permission flag that is set when identification is possible, and determines that the identification processing is permitted, for example, when the identification permission flag is "<NUM>". Here, whether it is the timing for identification may be used as the determination condition. If the system identification unit <NUM> determines that identification is impossible (NO in S2101), the processing of this flowchart ends.

Next, when determining that identification is possible (YES in S2101), the system identification unit <NUM> calculates Expressions (<NUM>) to (<NUM>) by the system identification algorithm (S2102 to S2106). <NUM>] <MAT>
[Math. <NUM>] <MAT>
[Math. <NUM>] <MAT>
[Math. <NUM>] <MAT>
[Math. <NUM>] <MAT>.

ε(k) in Expression (<NUM>) in Step S2102 represents a difference between the current measured value of the oxygen sensor <NUM> and the estimated value (prediction value) of the rear oxygen sensor. Equations (<NUM>) to (<NUM>) in Steps S2103 to S2106 are calculated based on the value of ε(k). θ^(k) (^ is on θ in the drawing) of Expression (<NUM>) in Step S2104 is a parameter vector of the statistical model <NUM>, and means an instantaneous value that changes sequentially. λ(k) in Expression (<NUM>) in Step S2105 represents a variable forgetting coefficient. Further, σ in Expression (<NUM>) is a parameter for adjusting the degree of variation of the variable forgetting coefficient λ(k). P(k) in Expression (<NUM>) in Step S2106 represents a covariance matrix defined by the regression vector.

These expressions, tuning parameters, and the like in the system identification are an example. Further, in Expressions (<NUM>) to (<NUM>), unlike Expression (<NUM>), a matrix of φ is transposed instead of θ, but any matrix may be used.

After the process of Step S2106, the process proceeds to the determination process of Step S2101. Then, the system identification unit <NUM> updates the tuning parameter of the statistical model <NUM> on the basis of the calculation result.

As described above, the statistical model learning unit (system identification unit <NUM>) of this embodiment is obtained by applying the sequential least squares algorithm. Since the sequential least squares algorithm is a method (sequential calculation equation) of correcting the immediately preceding estimated value each time new data is obtained, it is possible to perform online estimation and real-time estimation of the rear oxygen sensor. In particular, since the successive least squares approximation algorithm with a variable forgetting element sets an evaluation function that is more important for data closer to the identification time, applicability to a system in which parameters change as in this embodiment is high.

Note that, although the sequential least squares algorithm is adopted as the parameter identification algorithm of this embodiment, the present invention is not limited thereto. That is, even when other optimization methods such as a gradient method and a genetic algorithm are applied as the parameter identification algorithm, a similar or close effect is obtained.

<FIG> is a diagram for explaining the air-fuel ratio correction amount calculation in the air-fuel ratio correction amount calculation unit <NUM> shown in <FIG>. In each graph of <FIG>, the horizontal axis represents time, and the vertical axis represents an equivalent ratio (upper graph), a rear oxygen sensor output (middle graph), and a catalyst downstream NOx concentration (lower graph).

In each graph, a solid line indicates a behavior during past control, and a broken line indicates a behavior by future control. Here, a calculation method and an effect of the air-fuel ratio correction amount calculation unit <NUM> will be described using the air-fuel ratio correction after the fuel cut as an example.

The air-fuel ratio correction amount calculation unit <NUM> predicts the rear oxygen sensor output behavior in the n-step future based on the statistical model <NUM>. As the tuning parameter of the statistical model <NUM> used at this time, the latest dynamic characteristic of the post-processing system is considered by system identification. The air-fuel ratio correction amount calculation unit <NUM> controls the completion timing of the rich correction after the fuel cut (n-step future from the current time k) in consideration of the timing (time k + n) at which the prediction value of the rear oxygen sensor output crosses A threshold Th set in advance. The current output value of the statistical model <NUM> is expressed by Expression (<NUM>), and the future output value of the statistical model <NUM> is expressed by Expression (<NUM>). <NUM>] <MAT>
[Math. <NUM>] <MAT>.

The air-fuel ratio correction amount calculation unit <NUM> may calculate a prediction value of the rear oxygen sensor output by the statistical model <NUM> for each step and determine whether to perform the rich correction on the basis of the calculation result. Alternatively, the prediction value of the rear oxygen sensor output by the statistical model <NUM> may be calculated collectively for n steps to determine whether to perform the rich correction. Alternatively, the prediction value of the rear oxygen sensor output may be calculated collectively for n steps, and in a case where the prediction value of the rear oxygen sensor output exceeds the threshold, a step ((n-m) in the future) in which the prediction value of the rear oxygen sensor output does not exceed the threshold may be determined as the rich correction completion timing.

As described above, the air-fuel ratio correction amount calculation unit (air-fuel ratio correction amount calculation unit <NUM>) of this embodiment calculates the future catalyst downstream exhaust gas concentration information based on the learned statistical model (statistical model <NUM>), and calculates the current air-fuel ratio correction amount or the correction period (completion timing) of the air-fuel ratio based on the variation amount (difference, ratio, or the like) between the future catalyst downstream exhaust gas concentration and the target value of the catalyst downstream exhaust gas concentration.

According to this embodiment configured as described above, in the method of stopping the air-fuel ratio correction after the output of the exhaust gas sensor (for example, the rear oxygen sensor) on the downstream side of the catalyst has actually changed, the problem of excessive rich correction can be appropriately avoided, and NOx emission can be suppressed, and excessive rich correction can also be prevented. Therefore, the emission of CO and HC can also be suppressed.

Further, since the completion timing of the rich correction is optimized on-board according to the deterioration state of the catalyst <NUM> and the exhaust gas sensor on the downstream side of the catalyst, the adaptation work for each deterioration catalyst is unnecessary, and the adaptation development man-hours can be greatly reduced. Further, the catalyst purification efficiency can be maintained at a high level by appropriate air-fuel ratio control, which also contributes to a reduction in the amount of catalyst used, that is, a reduction in the cost of the post-processing system.

<FIG> is a flowchart illustrating an example in which the ECU <NUM> calculates an air-fuel ratio correction amount based on a flow rate sensor detection value and a catalyst upstream/downstream exhaust gas sensor detection value and executes fuel injection amount control.

First, the control unit <NUM> of the ECU <NUM> detects the intake air amount taken into the internal combustion engine <NUM> by the flow rate sensor <NUM> (S2301). Next, the catalyst boundary condition calculation unit <NUM> of the ECU <NUM> detects the exhaust gas composition (catalyst upstream air-fuel ratio, catalyst downstream oxygen presence/absence state) by the exhaust gas sensors (Air-fuel ratio sensor <NUM> and oxygen sensor <NUM> in <FIG>) assembled upstream and on the downstream side of the catalyst <NUM> (S2302).

Next, the catalyst temperature calculation unit <NUM> calculates the catalyst temperature based on the current operating state of the internal combustion engine <NUM> (S2303). Next, the oxygen storage ratio calculation unit <NUM> calculates the oxygen storage ratio on the basis of the information of the catalyst upstream exhaust gas composition, the catalyst temperature, the exhaust gas flow rate, and the oxygen storage capacity (S2304).

Next, the control unit <NUM> of the ECU <NUM> determines whether system identification can be performed on the basis of the operating state of the internal combustion engine <NUM> (S2305). The control unit <NUM> sets the identification permission flag to "<NUM> (ON)" when identification is possible, and sets the identification permission flag to "<NUM> (OFF)" when identification is impossible.

In a case where the control unit <NUM> determines that identification is impossible (NO in S2305), the processing proceeds to the determination processing in Step S2309. For example, an exhaust gas sensor operating state, a catalyst temperature range, an execution state of diagnosis processing, and the like are set as determination criteria of identification propriety. For example, the system identification processing is stopped before the start of the exhaust gas sensor immediately after the start of the engine. In addition, when the catalyst temperature is outside the predetermined range, the catalyst does not function normally, so that accurate system identification cannot be executed, and the identification processing is stopped. Further, the identification processing is stopped also at the time of diagnosis of the exhaust gas sensor related to the catalyst or the post-processing system.

As described above, the statistical model learning unit (system identification unit <NUM>) of this embodiment determines whether to learn the statistical model (statistical model <NUM>) on the basis of any one of the operating state of the first exhaust gas sensor (for example, the air-fuel ratio sensor <NUM>) and a second exhaust gas sensor (for example, the oxygen sensor <NUM>), the catalyst temperature state, the catalyst deterioration diagnosis state, and the diagnosis states of the first exhaust gas sensor and the second exhaust gas sensor, and stops learning of the statistical model when it is determined that learning is impossible. As a result, it is possible to efficiently perform learning while avoiding a period in which the learning effect does not increase.

In a case where the control unit <NUM> determines that identification is possible (Yes in S2305), the system identification unit <NUM> sequentially calculates the tuning parameter of the statistical model <NUM> on the basis of the successive least squares approximation algorithm with variable forgetting element to perform system identification (S2306).

In this embodiment, as an example, the system identification unit <NUM> (statistical model learning unit) is configured to learn the tuning parameter of the statistical model <NUM> during the fuel cut operation according to the determination result of the control unit <NUM>. At the time of fuel cut, the oxygen concentration in the exhaust gas increases, and the state of input and output fluctuates, so that highly accurate learning can be performed. Therefore, the rich correction (air-fuel ratio control) after the fuel cut can be more appropriately performed.

Next, the system identification unit <NUM> determines whether system identification is completed based on the update status of the tuning parameter (the temporal change amount of the tuning parameter is equal to or less than a predetermined value) (S2307). In a case where the identification is not completed, the tuning parameter changes as the input and the output change. Therefore, even if the input and the output fluctuate, if the fluctuation of the tuning parameter is small, it can be determined that the identification is completed. Next, in a case of determining that the system identification is not completed (NO in S2307), the system identification unit <NUM> returns to Step S2306 and continues the system identification processing.

Next, when it is determined that the system identification is completed (YES in S2307), the system identification unit <NUM> updates the tuning parameter of the statistical model <NUM> (S2308).

Next, in a case where it is determined in Step S2305 or after the process of S2308, the deterioration degree determination unit <NUM> determines the deterioration state of the catalyst <NUM> based on the identified tuning parameter (S2309). The tuning parameters include information on the dynamic characteristics of the catalyst <NUM> and the exhaust gas sensor, and it is possible to determine whether the deterioration state of the catalyst <NUM> is appropriate (whether it is within an assumption of the statistical model <NUM>) by capturing a change in the dynamic characteristics.

When it is determined that the deterioration of the catalyst <NUM> has progressed (NO in S2309), the deterioration degree determination unit <NUM> decreases the value of the oxygen storage capacity in the oxygen storage ratio calculation in Step S2304. Then, the processing of Steps S2304 to S2309 is executed again.

On the other hand, when it is determined that the catalyst <NUM> is not deteriorated (YES in S2309), the statistical model calculation unit <NUM> calculates (predicts) the future exhaust gas sensor output on the downstream side of the catalyst based on the statistical model <NUM> configured by the updated tuning parameter (S2310).

Next, the air-fuel ratio correction amount calculation unit <NUM> calculates the air-fuel ratio correction amount based on the future exhaust gas sensor output on the downstream side of the catalyst by the statistical model <NUM> (S2311).

Next, the target air-fuel ratio calculation unit <NUM> calculates the target air-fuel ratio by adding the air-fuel ratio correction amount to the air-fuel ratio control amount. Then, the control unit <NUM> of the ECU <NUM> calculates the fuel injection amount by the fuel injection valve <NUM> based on the intake air amount so that the air-fuel mixture in the internal combustion engine <NUM> becomes the target air-fuel ratio (S2312). After the process of Step S2312 ends, the process of this flowchart ends.

Note that, in a case where the statistical model <NUM> is not updated, the processing regarding system identification in Steps S2305 to S2308 can be omitted. However, as described in the description of <FIG>, by controlling the air-fuel ratio using the statistical model <NUM> updated to the latest state based on the measured value of the catalyst downstream exhaust gas sensor, the oxygen storage state in the catalyst <NUM> can be always maintained in a suitable state.

When the deterioration of the catalyst <NUM> is not considered, the deterioration determination process of the catalyst <NUM> in Step S2309 can be omitted.

As described above, the internal combustion engine control device (ECU <NUM>) of this embodiment is an internal combustion engine control device that controls an internal combustion engine including the catalyst (exhaust purification catalyst <NUM>) disposed in the exhaust pipe, the first exhaust gas sensor (for example, air-fuel ratio sensor <NUM>) disposed on the upstream side of the catalyst, and the second exhaust gas sensor (for example, oxygen sensor <NUM>) disposed on the downstream side of the catalyst. This internal combustion engine control device includes: an oxygen storage ratio calculation unit (oxygen storage ratio calculation unit <NUM>) that calculates an oxygen storage ratio of a catalyst based on a catalytic reaction model (for example, Reaction Expressions (<NUM>) to (<NUM>)) having at least a detection value of a first exhaust gas sensor as an input; a statistical model calculation unit (statistical model calculation unit <NUM>) that predicts a catalyst downstream exhaust gas concentration using a statistical model (statistical model <NUM>) having an oxygen storage ratio (u) as an input and a catalyst downstream exhaust gas concentration (y) as an output; and an air-fuel ratio correction amount calculation unit (air-fuel ratio correction amount calculation unit <NUM>) that calculates an air-fuel ratio correction amount of an air-fuel mixture of the internal combustion engine based on a future catalyst downstream exhaust gas concentration calculated by the statistical model calculation unit.

In this embodiment configured as described above, the oxygen storage ratio based on the catalytic reaction model (physical model) is input to the statistical model, and the catalyst downstream exhaust gas concentration is predicted. By using the catalytic reaction model which is a physical model, the catalyst downstream exhaust gas concentration can be predicted up to a slightly earlier time (for example, several <NUM> seconds) according to the physical law. On the other hand, by using the statistical model, fine matching can be performed for a certain limited time (shorter time). As described above, by combining the catalytic reaction model (physical model) and the statistical model, the catalyst purification efficiency can be kept high, and deterioration of emission performance can be prevented.

In addition, the internal combustion engine control device (ECU <NUM>) of this embodiment further includes a statistical model learning unit (system identification unit <NUM>) that learns (system identification) tuning parameters of the statistical model (statistical model <NUM>) based on the detection value of the second exhaust gas sensor (for example, the oxygen sensor <NUM>). With this configuration, the statistical model can always be updated to the latest state.

With such a configuration, the statistical model <NUM> (tuning parameter) having the oxygen storage ratio of the catalyst <NUM> calculated based on the catalytic reaction model as an input and the catalyst downstream exhaust gas concentration as an output is always updated to the latest state based on the measured value of the catalyst downstream exhaust gas sensor. Then, since the air-fuel ratio can be controlled based on the future catalyst downstream exhaust gas concentration predicted by the updated statistical model <NUM>, the oxygen storage state in the catalyst <NUM> can be always maintained in a suitable state. As a result, even when the catalyst <NUM> or the catalyst downstream exhaust gas sensor deteriorates, deterioration of emission performance can be appropriately prevented.

In addition, the internal combustion engine control device (ECU <NUM>) of this embodiment further includes a deterioration degree determination unit (deterioration degree determination unit <NUM>) that determines the deterioration degree (catalyst deterioration degree) of the catalyst based on the tuning parameter of the learned statistical model (statistical model <NUM>). For example, the deterioration degree determination unit may output the catalyst deterioration diagnosis result based on the deterioration degree and the deterioration state assumed in the statistical model.

The determination of the degree of catalyst deterioration in this manner contributes to setting of the oxygen storage capacity (ψcap) according to the catalyst deterioration degree in the oxygen storage ratio calculation unit (oxygen storage ratio calculation unit <NUM>). For example, by outputting the stepwise evaluation of the catalyst deterioration diagnosis result (for example, new, not so new, old) to the oxygen storage ratio calculation unit, the oxygen storage capacity (ψcap) can be set stepwise in the oxygen storage ratio calculation unit, and the oxygen storage ratio can be calculated based on the state of the catalyst.

In addition, the oxygen storage ratio calculation unit (oxygen storage ratio calculation unit <NUM>) of this embodiment has, as a parameter, an oxygen storage capacity corresponding to the amount of metal (platinum, rhodium, etc.) carried in the catalyst <NUM>, and sets the oxygen storage capacity according to the catalyst deterioration degree.

As described above, by setting the oxygen storage capacity according to the catalyst deterioration, the accuracy of the calculation result of the oxygen storage ratio is improved, and a highly accurate oxygen storage ratio can be input by the statistical model (statistical model <NUM>).

Further, the present invention is not limited to each of the above-described embodiments, and it goes without saying that various other application examples and modifications can be made within the scope of the appended claims.

For example, the above embodiments describe the engine system in detail and specifically for easy understanding the present invention, but not necessarily limited to those including all the components described. In addition, some of the configurations of each embodiment may be omitted, replaced with other components, and added to other configurations.

Claim 1:
An internal combustion engine control device that controls an internal combustion engine (<NUM>) using a calculated air-fuel ratio correction amount, the internal combustion engine (<NUM>) including a catalyst (<NUM>) disposed in an exhaust pipe, a first exhaust gas sensor disposed on an upstream side of the catalyst (<NUM>), and a second exhaust gas sensor disposed on a downstream side of the catalyst (<NUM>), the internal combustion engine control device (<NUM>) comprising:
An oxygen storage ratio calculation unit (<NUM>) that calculates an oxygen storage ratio of the catalyst (<NUM>) based on a catalytic reaction model having at least a detection value of the first exhaust gas sensor as an input;
a statistical model calculation unit (<NUM>) that predicts a catalyst downstream exhaust gas concentration by using a statistical model (<NUM>) having the oxygen storage ratio as an input and the catalyst downstream exhaust gas concentration as an output; and
an air-fuel ratio correction amount calculation unit (<NUM>) that calculates an air-fuel ratio correction amount of an air-fuel mixture of the internal combustion engine (<NUM>) based on a future catalyst downstream exhaust gas concentration calculated by the statistical model calculation unit (<NUM>),
characterized in that
the catalytic reaction model is defined by a reaction rate between at least oxygen on the upstream side of the catalyst (<NUM>), carbon monoxide on the upstream side of the catalyst (<NUM>), and a metal carried in the catalyst, wherein
the oxygen storage ratio calculation unit (<NUM>) includes:
an oxygen concentration calculation unit (<NUM>) that calculates an oxygen concentration on the upstream side of the catalyst (<NUM>() based on a detection value of a first exhaust gas sensor on the upstream side of the catalyst (<NUM>) ; and
a carbon monoxide concentration calculation unit (<NUM>) that calculates a carbon monoxide concentration on the upstream side of the catalyst (<NUM>) based on the detection value of the first exhaust gas sensor on the upstream side of the catalyst (<NUM>).