实施方案
[0025] 下面结合附图进一步说明本发明。
[0026] 如图1所示，一种天然气掺水燃烧发动机，包括天然气发动机本体和掺水系统；掺水系统由水罐22、水泵21、高压油泵8和注水器10组成；注水器10的结构可以与现有技术的喷油器的结构完全相同；液态水存储于水罐22中，水罐22的出水口接水泵21的入口，水泵21的出口接高压油泵8的入口，高压油泵8的出口接注水器10的入口，注水器10的喷射口插入发动机气缸的燃烧室；液态水经水泵21和高压油泵8加压后，由注水器10喷射进入发动机气缸的燃烧室；水泵21和高压油泵8的启停以及注水器10内电磁阀的开闭均由电子控制单元7控制；电磁阀的开闭决定注水器10的喷射时刻与喷射脉宽。
[0027] 天然气发动机本体包括空气滤清器1、涡轮增压器、进气压力温度传感器3、节气门4、燃气喷射器5、凸轮轴位置传感器6、电子控制单元(ECU) 7、火花塞9、排气温度传感器11、氧传感器12、柴油氧化催化器(DOC) 13、天然气罐14、减压阀15、燃气驱动器16、燃气压力温度传感器17、冷却水温传感器18、爆震传感器19和转速传感器20；进气压力温度传感器3可以采用博世的型号为0261230245的进气压力温度传感器。
[0028] 空气滤清器1的出气口接进气管管段一的进气口，进气管管段一的出气口接涡轮增压器的压气端2-1进气口；新鲜空气经空气滤清器1过滤后，进入涡轮增压器的压气端2-1进气口，由涡轮增压器进行增压；涡轮增压器的压气端2-1出气口接进气管管段二的进气口，进气管管段二的出气口接发动机气缸的进气口；进气管管段二上设有节气门4和两个进气压力温度传感器3，两个进气压力温度传感器3设置在节气门4两端；进气压力温度传感器检测进气管管段二内压力和温度，节气门4控制进气管的气流量。
[0029] 发动机气缸的排气口接排气管管段一的进气口，发动机气缸燃烧后的废气经由排气管管段一排出，排气管管段一的出气口接涡轮增压器的涡轮端2-2进气口，涡轮增压器的涡轮端2-2出气口接排气管管段二的进气口，排气管管段二的出气口接柴油氧化催化器13的进气口，柴油氧化催化器13 的出气口开放设置；排气管管段二上设有排气温度传感器11和两个氧传感器12；其中，涡轮增压器的涡轮端为增压提供能量，排气温度传感器11测量排气温度，氧传感器12检测混合气空燃比，柴油氧化催化器13处理尾气中的碳氢污染物。
[0030] 天然气存储于天然气罐14，天然气罐14的出气口经减压阀15接燃气驱动器16，燃气驱动器16的出口接燃气喷射器5；燃气喷射器5的喷嘴插入进气管管段二内腔；燃气驱动器16的出口处设有燃气压力温度传感器17；减压阀15将天然气的压力降低到喷射压力，随后通过燃气驱动器16与燃气喷射器5，将天然气喷射进入进气管管段二；燃气压力温度传感器17检测天然气的温度与压力。
[0031] 此外，发动机气缸的水套上设有冷却水温度传感器18，监测冷却水温度；发动机气缸的缸盖上设有爆震传感器19，监测发动机是否发生爆震；凸轮轴两端均设置凸轮轴位置传感器6，为发动机点火、注水器喷射等提供相位信号；电子控制单元7通过节气门4和设置在曲轴上的转速传感器20，确定发动机工作工况，通过查找MAP图的方式，确定注水器10的喷射时刻与喷射脉宽；同时，电子控制单元7利用PID算法，对注水器10的喷射脉宽进行修正，以达到当前工况下的预期排温。
[0032] 如图2所示，该天然气掺水燃烧发动机的掺水控制方法，具体如下：
[0033] 电子控制单元内的PID控制器根据发动机转速(转速传感器测得)和节气门开度(节气门开度传感器测得)，从MAP图中读取当前工况 n(PID控制器的第n个采样周期)的预期排温Th(n)以及注水器10的喷射时刻和喷射脉宽初值，当n＝1或n＝2时，电子控制单元直接控制注水器内的电磁阀按喷射脉宽初值在喷射时刻打开，实现发动机气缸内直接喷水。当n ≥3时，按以下步骤对喷射脉宽进行修正：
[0034] 第一步：PID控制器根据排气温度传感器11传来的信号获取当前工况n 的实际排温Tf(n)。
[0035] 第二步：PID控制器计算当前工况n的预期排温与实际排温之间的温差 Td(n)＝Th(n)-Tf(n)。
[0036] 第三步：将Td(n)作为输入参数，通过PID控制器的PID控制算法计算注水器10的喷射脉宽修正值Δt(n)，对注水器的喷射脉宽(喷射的时间) 进行修正，然后控制注水器内的电磁阀按喷射脉宽初值与喷射脉宽修正值Δt(n)之和的喷射脉宽在喷射时刻打开，实现发动机气缸内直接喷水。当前工况的喷射脉宽初值更新为修正值后的喷射脉宽初值并写入MAP图中，PID 控制既能保证控制精度，又能提高系统的鲁棒性，PID控制算法如下：
[0037] Δt(n)＝a0Td(n)+a1Td(n-1)+a2Td(n-2)
[0038] 其中，工况n的前一工况的温差Td(n-1)＝Th(n-1)-Tf(n-1)，工况n的前两个工况的温差Td(n-2)＝Th(n-2)-Tf(n-2)，Th(n-1)和Th(n-2) 分别为工况n的前一工况和前两个工况的预期排温，Tf(n-1)和Tf(n-2)分别为工况n的前一工况和前两个工况的实际排温，各工况的温差计算好后便存于PID控制器中，以便计算下一工况时调用；Kp表示比例系数；ti表示积分时间，单位为秒；td表
示微分时间，单位为秒；t表示PID控制器的采样周期，单位为秒。
