CN110909958A - 一种计及光伏并网功率的短期负荷预测方法 - Google Patents
一种计及光伏并网功率的短期负荷预测方法 Download PDFInfo
- Publication number
- CN110909958A CN110909958A CN201911237147.1A CN201911237147A CN110909958A CN 110909958 A CN110909958 A CN 110909958A CN 201911237147 A CN201911237147 A CN 201911237147A CN 110909958 A CN110909958 A CN 110909958A
- Authority
- CN
- China
- Prior art keywords
- power
- data
- load
- grid
- short
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Abstract
本发明提供了一种计及光伏并网功率的短期负荷预测方法,包括如下步骤:调取电网的历史负荷数据,根据历史负荷数据变化情况及变化时刻的外部影响情况,确定电网电力负荷的关键影响因素;构建短期电力负荷预测模型;从电网的历史电力负荷数据库中选取训练数据;对选取的负荷序列进行初步数据清洗;对数据进行约束化处理,对电力负荷预测模型的参数进行设置,根据训练数据对模型进行训练;根据模型训练后得到的负荷预测方程,对电网的电力负荷进行短期预测;对光伏并网的功率进行短期预测;将预测的电网电力负荷与光伏并网功率结果进行叠加,得到计及光伏并网功率的电网短期等效预测负荷。
Description
技术领域
本发明属于电力负荷预测技术领域,特别涉及一种计及光伏并网功率的短期负荷预测方法。
背景技术
电力负荷预测在现代电力系统中具有非常重要的意义,是影响电力系统安全性、经济性的重要因素。电力负荷是指一个地区单位时间内消耗的电能。在一个区域内,只要是联网的用电设备都会消耗电能,消耗的总能量即为该地区的总负荷。现代经济社会的发展不仅对电力的需求量越来越高,同时也对电能的质量提出了更高的要求。影响电能总量需求主要是各类工业、农业和居民用电需求,这种需求不仅受时间的影响,还受天气、节假日因素的影响,具有非常复杂的变化规律。影响电能质量的因素是各发电厂的上网电能的波动,这种波动性虽然受发电设备的影响较大,但影响波动性的关键因素是用电户的用电习惯。
近些年来,光伏发电技术越来越受到到人们的关注,因为光伏发电随机性、不确定性的独有属性,光伏发电接入电网会引起安全和稳定问题,准确的光伏发电功率预测技术能够对电网提供更加准确的数据,有利于电网的进一步分配与调度。
发明内容
本发明提供一种计及光伏并网功率的短期负荷预测方法,对电网电力负荷和光伏并网功率分别进行短期预测,并将预测的电网电力负荷与光伏并网功率结果进行叠加,得到计及光伏并网功率的电网短期等效预测负荷。
本发明具体为一种计及光伏并网功率的短期负荷预测方法,所述短期负荷预测方法具体包括如下步骤:
步骤(1)、调取电网的历史负荷数据,根据历史负荷数据变化情况及变化时刻的外部影响情况,确定电网电力负荷的关键影响因素;
步骤(2)、根据步骤(1)中确定的电网电力负荷的关键影响因素构建短期电力负荷预测模型;
步骤(3)、从电网的历史电力负荷数据库中选取训练数据;
步骤(4)、对选取的负荷序列进行初步数据清洗,包括序列噪声抑制、缺失数据修补以及异常数据修正;
步骤(5)、对数据进行约束化处理,对电力负荷预测模型的参数进行设置,根据训练数据对模型进行训练;
步骤(6)、根据模型训练后得到的负荷预测方程,对电网的电力负荷进行短期预测;
步骤(7)、采用最小二乘支持向量机模型对光伏并网的功率进行短期预测;
步骤(8)、将预测的电网电力负荷与光伏并网功率结果进行叠加,得到计及光伏并网功率的电网短期等效预测负荷。
所述步骤(1)中的电网电力负荷的关键影响因素包括气温、湿度、天气类型、季节类型、日期信息、政策信息、上一时刻的负荷值、前一天同一时刻的负荷值。
所述步骤(2)中的短期电力负荷预测模型基于最小二乘支持向量机回归模型和径向基核函数确定:y=(q,φ(X))+b,目标函数为其中q为要寻求的最优权值,b为线性函数阈值,X=[x1,x2,…,x8]为一个8维向量,表示输入量;y为输出数据,单位为MW;ei为允许误差,F为惩罚因子;x1为预测时刻的预测气温;x2为预测时刻的预测湿度;x3为预测当天的天气类型,分别用数字量0-5表示晴天、多云、阴天、雨天、雪天、台风;x4为预测时刻的季节,分别用数字量0-3表示春天、夏天、秋天、冬天;x5为预测当天的日期信息,表示当天是否为周末或节假日,分别用数字量0、1表示工作日、休息日;x6为预测当天的政策信息,表示当天是否有重大事件,分别用数字量0、1表示没有重大事件和有重大事件;x7为预测时刻上一时刻的负荷值;x8为预测前一天同一时刻的负荷值。
所述步骤(4)中的序列噪声抑制采用模态分解方式,最先分解得到的模态具有最小的时间尺度,由于电力负荷序列呈现随机分布特性,将第一模态认定为噪声将其去除。
所述步骤(4)中的缺失数据修补从相似日方面和时间序列方面两个方面进行修补,对于相似日预测,将连续一段时间的数据按一天进行划分,则负荷序列从一个行向量转变为矩阵形式:
所述步骤(4)中的异常数据修正具体包括明显异常数据剔除和缺失数据修补,对于明显异常的数据将其剔除出序列,并采用缺失数据修补的方式重新给出新的数据来代替异常数据。
所述步骤(7)中采用最小二乘支持向量机模型对光伏并网的功率进行短期预测具体包括:
获取光伏并网功率历史数据,并对数据进行预处理,对缺测数据进行填补,对不合理数据进行修正;
对数据集进行天气类型的划分,寻找近期日类型相同的数据,将相同类型相同时刻的历史输出功率值、以及温度湿度气象信息,形成训练样本,作为输入数据;
对数据进行归一化处理;
选择径向基核函数,采用参数搜索和交叉验证,选择最佳的惩罚系数和核函数参数;
用得到的最佳的惩罚系数和核函数参数进行训练,得到最小二乘支持向量机功率预测模型;
选择晴天、多云、雨天和雾霾天四种类型的天气,分别对各整点的输出功率进行预测,验证并分析预测结果。
附图说明
图1为本发明一种计及光伏并网功率的短期负荷预测方法的流程图。
具体实施方式
下面结合附图对本发明一种计及光伏并网功率的短期负荷预测方法的具体实施方式做详细阐述。
如图1所示,本发明的预测方法包括以下步骤:步骤(1)、调取电网的历史负荷数据,根据历史负荷数据变化情况及变化时刻的外部影响情况,确定电网电力负荷的关键影响因素;步骤(2)、根据步骤(1)中确定的电网电力负荷的关键影响因素构建短期电力负荷预测模型;步骤(3)、从电网的历史电力负荷数据库中选取训练数据;步骤(4)、对选取的负荷序列进行初步数据清洗,包括序列噪声抑制、缺失数据修补以及异常数据修正;步骤(5)、对数据进行约束化处理,对电力负荷预测模型的参数进行设置,根据训练数据对模型进行训练;步骤(6)、根据模型训练后得到的负荷预测方程,对电网的电力负荷进行短期预测;步骤(7)、采用最小二乘支持向量机模型对光伏并网的功率进行短期预测;步骤(8)、将预测的电网电力负荷与光伏并网功率结果进行叠加,得到计及光伏并网功率的电网短期等效预测负荷。
电网电力负荷的关键影响因素包括气温、湿度、天气类型、季节类型、日期信息、政策信息、上一时刻的负荷值、前一天同一时刻的负荷值。用一个8维向量X=[x1,x2,…,x8]来表示输入量,x1为预测时刻的预测气温;x2为预测时刻的预测湿度;x3为预测当天的天气类型,分别用数字量0-5表示晴天、多云、阴天、雨天、雪天、台风;x4为预测时刻的季节,分别用数字量0-3表示春天、夏天、秋天、冬天;x5为预测当天的日期信息,表示当天是否为周末或节假日,分别用数字量0、1表示工作日、休息日;x6为预测当天的政策信息,表示当天是否有重大事件,分别用数字量0、1表示没有重大事件和有重大事件;x7为预测时刻上一时刻的负荷值;x8为预测前一天同一时刻的负荷值。
支持向量回归理论中,假设xi∈Rn为输入,yi∈R为相应的输出,回归问题就是寻找从输入到输出的映射f:Rn→R,使得f(x)=y。简单的线性回归问题为y=f(x)=q·x+b,其目的就是寻找最优的权值q,使得拟合曲线尽可能反应数据集的变化规律,对于优化问题,就是以精度ε可拟合所有训练样本。回归问题转化为优化问题为:约束条件为|<q,xi>+b-yi|≤ε,i=1,…,l。
本发明中的短期电力负荷预测模型基于最小二乘支持向量机回归模型和径向基核函数确定:y=(q,φ(X))+b,目标函数为其中q为要寻求的最优权值,b为线性函数阈值,X=[x1,x2,…,x8]为一个8维向量,表示输入量;y为输出数据,单位为MW;ei为允许误差,F为惩罚因子。约束条件为|<q,xi>+b-yi|≤ε,i=1,…,l及其中,ε为精度。
序列噪声抑制采用模态分解方式,最先分解得到的模态具有最小的时间尺度,由于电力负荷序列呈现随机分布特性,将第一模态认定为噪声将其去除。缺失数据修补从相似日方面和时间序列方面两个方面进行修补,对于相似日预测,将连续一段时间的数据按一天进行划分,则负荷序列从一个行向量转变为矩阵形式:每一列表示相似日序列;假设和分别表示缺失数据xt在垂直和水平两个方向的修正结果,则最终修正结果为对于非连续缺失数据的修补,仅从相似日方面进行修补。异常数据修正具体包括明显异常数据剔除和缺失数据修补,对于明显异常的数据将其剔除出序列,并采用缺失数据修补的方式重新给出新的数据来代替异常数据。
本发明同样采用上述最小二乘支持向量机模型对光伏并网的功率进行短期预测,具体包括:获取光伏并网功率历史数据,并对数据进行预处理,对缺测数据进行填补,对不合理数据进行修正;对数据集进行天气类型的划分,寻找近期日类型相同的数据,将相同类型相同时刻的历史输出功率值、以及温度湿度气象信息,形成训练样本,作为输入数据;对数据进行归一化处理;选择径向基核函数,采用参数搜索和交叉验证,选择最佳的惩罚系数和核函数参数;用得到的最佳的惩罚系数和核函数参数进行训练,得到最小二乘支持向量机功率预测模型;选择晴天、多云、雨天和雾霾天四种类型的天气,分别对各整点的输出功率进行预测,验证并分析预测结果。基于24小时天气预报,可以获得预测日的天气状况,包括最高温度、最低温度和平均温度,以及湿度等气象信息。判断天气类型,选择类型相同的前五天相同时刻的输出功率值作为预测模型输入,对晴天、多云、雨天和雾霾天分别经行预测,对于晴天的预测,我们选取同样是晴天并且和预测日时间最相近的五天的数据作为预测模型的输入变量,对于多云的预测,我们选取同样是多云并且和预测日时间最相近的五天的数据作为预测模型的输入变量,对于雨天的预测,我们选取同样是雨天并且和预测日时间最相近的五天的数据作为预测模型的输入变量,对于雾霾天的预测,我们选取同样是雾霾天并且和预测日时间最相近的五天的数据作为预测模型的输入变量。
在分别得到电网的短期电力负荷预测数据与光伏并网的短期功率预测数据后,将预测的电网电力负荷与光伏并网功率结果进行叠加,得到计及光伏并网功率的电网短期等效预测负荷。此等效预测负荷与不考虑光伏并网功率而直接预测的电力负荷相比,预测数据更为准确实用。
最后应该说明的是,结合上述实施例仅说明本发明的技术方案而非对其限制。所属领域的普通技术人员应当理解到,本领域技术人员可以对本发明的具体实施方式进行修改或者等同替换,但这些修改或变更均在申请待批的权利要求保护范围之中。
Claims (8)
1.一种计及光伏并网功率的短期负荷预测方法,其特征在于,所述短期负荷预测方法具体包括如下步骤:
步骤(1)、调取电网的历史负荷数据,根据历史负荷数据变化情况及变化时刻的外部影响情况,确定电网电力负荷的关键影响因素;
步骤(2)、根据步骤(1)中确定的电网电力负荷的关键影响因素构建短期电力负荷预测模型;
步骤(3)、从电网的历史电力负荷数据库中选取训练数据;
步骤(4)、对选取的负荷序列进行初步数据清洗,包括序列噪声抑制、缺失数据修补以及异常数据修正;
步骤(5)、对数据进行约束化处理,对电力负荷预测模型的参数进行设置,根据训练数据对模型进行训练;
步骤(6)、根据模型训练后得到的负荷预测方程,对电网的电力负荷进行短期预测;
步骤(7)、采用最小二乘支持向量机模型对光伏并网的功率进行短期预测;
步骤(8)、将预测的电网电力负荷与光伏并网功率结果进行叠加,得到计及光伏并网功率的电网短期等效预测负荷。
2.根据权利要求1所述的一种计及光伏并网功率的短期负荷预测方法,其特征在于,所述步骤(1)中的电网电力负荷的关键影响因素包括气温、湿度、天气类型、季节类型、日期信息、政策信息、上一时刻的负荷值、前一天同一时刻的负荷值。
3.根据权利要求2所述的一种计及光伏并网功率的短期负荷预测方法,其特征在于,所述步骤(2)中的短期电力负荷预测模型基于最小二乘支持向量机回归模型和径向基核函数确定:y=(q,φ(X))+b,目标函数为其中q为要寻求的最优权值,b为线性函数阈值,X=[x1,x2,…,x8]为一个8维向量,表示输入量;y为输出数据,单位为MW;ei为允许误差,F为惩罚因子;x1为预测时刻的预测气温;x2为预测时刻的预测湿度;x3为预测当天的天气类型,分别用数字量0-5表示晴天、多云、阴天、雨天、雪天、台风;x4为预测时刻的季节,分别用数字量0-3表示春天、夏天、秋天、冬天;x5为预测当天的日期信息,表示当天是否为周末或节假日,分别用数字量0、1表示工作日、休息日;x6为预测当天的政策信息,表示当天是否有重大事件,分别用数字量0、1表示没有重大事件和有重大事件;x7为预测时刻上一时刻的负荷值;x8为预测前一天同一时刻的负荷值。
4.根据权利要求3所述的一种计及光伏并网功率的短期负荷预测方法,其特征在于,所述步骤(4)中的序列噪声抑制采用模态分解方式,最先分解得到的模态具有最小的时间尺度,由于电力负荷序列呈现随机分布特性,将第一模态认定为噪声将其去除。
6.根据权利要求5所述的一种计及光伏并网功率的短期负荷预测方法,其特征在于,所述步骤(4)中的异常数据修正具体包括明显异常数据剔除和缺失数据修补,对于明显异常的数据将其剔除出序列,并采用缺失数据修补的方式重新给出新的数据来代替异常数据。
8.根据权利要求7所述的一种计及光伏并网功率的短期负荷预测方法,其特征在于,所述步骤(7)中采用最小二乘支持向量机模型对光伏并网的功率进行短期预测具体包括:
获取光伏并网功率历史数据,并对数据进行预处理,对缺测数据进行填补,对不合理数据进行修正;
对数据集进行天气类型的划分,寻找近期日类型相同的数据,将相同类型相同时刻的历史输出功率值、以及温度湿度气象信息,形成训练样本,作为输入数据;
对数据进行归一化处理;
选择径向基核函数,采用参数搜索和交叉验证,选择最佳的惩罚系数和核函数参数;
用得到的最佳的惩罚系数和核函数参数进行训练,得到最小二乘支持向量机功率预测模型;
选择晴天、多云、雨天和雾霾天四种类型的天气,分别对各整点的输出功率进行预测,验证并分析预测结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911237147.1A CN110909958A (zh) | 2019-12-05 | 2019-12-05 | 一种计及光伏并网功率的短期负荷预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911237147.1A CN110909958A (zh) | 2019-12-05 | 2019-12-05 | 一种计及光伏并网功率的短期负荷预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110909958A true CN110909958A (zh) | 2020-03-24 |
Family
ID=69822596
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911237147.1A Withdrawn CN110909958A (zh) | 2019-12-05 | 2019-12-05 | 一种计及光伏并网功率的短期负荷预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110909958A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111931994A (zh) * | 2020-07-20 | 2020-11-13 | 深圳供电局有限公司 | 一种短期负荷及光伏功率预测方法及其系统、设备、介质 |
CN112132351A (zh) * | 2020-09-28 | 2020-12-25 | 华北电力大学 | 电动汽车充电负荷聚类预测方法及装置 |
CN112801374A (zh) * | 2021-01-29 | 2021-05-14 | 广东晨兴智能科技有限公司 | 模型训练方法、用电负荷预测方法、装置及设备 |
CN116167508A (zh) * | 2023-01-16 | 2023-05-26 | 南通电力设计院有限公司 | 一种基于气象因素分解的短期光伏出力快速预测方法及系统 |
CN117477581A (zh) * | 2023-12-26 | 2024-01-30 | 佛山市达衍数据科技有限公司 | 一种电力系统负荷均衡控制方法及电力系统 |
-
2019
- 2019-12-05 CN CN201911237147.1A patent/CN110909958A/zh not_active Withdrawn
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111931994A (zh) * | 2020-07-20 | 2020-11-13 | 深圳供电局有限公司 | 一种短期负荷及光伏功率预测方法及其系统、设备、介质 |
CN112132351A (zh) * | 2020-09-28 | 2020-12-25 | 华北电力大学 | 电动汽车充电负荷聚类预测方法及装置 |
CN112801374A (zh) * | 2021-01-29 | 2021-05-14 | 广东晨兴智能科技有限公司 | 模型训练方法、用电负荷预测方法、装置及设备 |
CN116167508A (zh) * | 2023-01-16 | 2023-05-26 | 南通电力设计院有限公司 | 一种基于气象因素分解的短期光伏出力快速预测方法及系统 |
CN116167508B (zh) * | 2023-01-16 | 2024-03-15 | 南通电力设计院有限公司 | 一种基于气象因素分解的短期光伏出力快速预测方法及系统 |
CN117477581A (zh) * | 2023-12-26 | 2024-01-30 | 佛山市达衍数据科技有限公司 | 一种电力系统负荷均衡控制方法及电力系统 |
CN117477581B (zh) * | 2023-12-26 | 2024-03-26 | 佛山市达衍数据科技有限公司 | 一种电力系统负荷均衡控制方法及电力系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110729764B (zh) | 一种含光伏发电系统的优化调度方法 | |
CN110880789B (zh) | 一种风电、光伏联合发电系统经济调度方法 | |
CN110909958A (zh) | 一种计及光伏并网功率的短期负荷预测方法 | |
CN112257941B (zh) | 基于改进型Bi-LSTM的光伏电站短期功率预测方法 | |
CN110866658A (zh) | 一种城市电网中长期负荷预测方法 | |
CN110097220B (zh) | 一种风力发电月度电量预测方法 | |
CN110991747A (zh) | 一种计及风电场功率的短期负荷预测方法 | |
CN111628499A (zh) | 一种计及多风险因素的配电网新能源消纳能力评估的方法 | |
CN110601250B (zh) | 一种基于新能源电站接入点静态电压稳定性评估方法 | |
CN111525552B (zh) | 一种基于特征信息的三阶段短期风电场群功率预测方法 | |
CN114004139A (zh) | 一种光伏发电功率预测方法 | |
CN114372360A (zh) | 用于电力负荷预测的方法、终端及存储介质 | |
CN114021420A (zh) | 一种分布式光伏超短期发电功率预测方法及系统 | |
CN114792156A (zh) | 基于曲线特征指标聚类的光伏输出功率预测方法和系统 | |
CN116227637A (zh) | 一种面向有源配电网的精细化负荷预测方法和系统 | |
CN115222106A (zh) | 一种自适应模型的用户日前负荷预测方法 | |
CN115965132A (zh) | 基于ga-bp神经网络的分布式光伏数字孪生系统功率预测方法 | |
CN115481918A (zh) | 一种基于源网荷储的单元状态主动感知及预测分析系统 | |
CN115545333A (zh) | 一种多负荷日类型配电网负荷曲线预测方法 | |
CN110991748A (zh) | 一种城市电网短期负荷预测方法 | |
CN109858668B (zh) | 一种雷电气候下电力负荷区域协调预测方法 | |
CN112307672A (zh) | 基于布谷鸟算法优化的bp神经网络短期风功率预测方法 | |
CN117273195A (zh) | 一种蒸汽供热系统需求负荷预测方法 | |
CN116485443A (zh) | 一种电力现货市场价格信息预测优化方法及系统 | |
CN115860797A (zh) | 一种适应新电价改革形势的电量需求预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20200324 |
|
WW01 | Invention patent application withdrawn after publication |