РАДИАЦИОННАЯ БИОЛОГИЯ. РАДИОЭКОЛОГИЯ, 2007, том 47, № 5, с. 517-522
ИЗУЧЕНИЕ ПОСЛЕДСТВИЙ АВАРИИ НА ЧЕРНОБЫЛЬСКОЙ АЭС
УДК 614.876:631.039.58
ЗАБОЛЕВАЕМОСТЬ РАКОМ ЩИТОВИДНОЙ ЖЕЛЕЗЫ СРЕДИ ЛИКВИДАТОРОВ ПОСЛЕДСТВИЙ АВАРИИ НА ЧЕРНОБЫЛЬСКОЙ АЭС: ПЕРИОД НАБЛЮДЕНИЯ 1986-2003 гг.
© 2007 г. В. К. Иванов*, С. Ю. Чекин, В. В. Кащеев, М. А. Максштов, К. А. Туманов
Медицинский радиологический научный центр РАМН, Обнинск
В последние годы в научной литературе было опубликовано значительное число работ, посвященных росту заболеваемости раком щитовидной железы (РЩЖ) среди детского населения территорий Белоруссии, России и Украины, загрязненных радионуклидами вследствие аварии на ЧАЭС. Вместе с тем в крайне ограниченном числе исследований рассматривается проблема заболеваемости РЩЖ среди ликвидаторов. В представленной работе изучается заболеваемость РЩЖ в когорте ликвидаторов (103427 человека), проживающих в 6 регионах России: Северо-Западном, Волго-Вятском, Центрально-Черноземном, Поволжском, Северо-Кавказском и Уральском. За период с 1986 по 2003 г. в этой когорте выявлено 87 случаев заболевания РЩЖ. Установлено статистически значимое превышение регистрируемой заболеваемости РЩЖ среди ликвидаторов над характерным для мужского населения России уровнем, SIR = 3.39 (95%-ный ДИ: 2.73; 4.16). Наиболее высокая заболеваемость РЩЖ (SIR = 6.49) зарегистрирована среди ликвидаторов, принимавщих участие в работах в апреле-июле 1986 г. Вместе с тем показано отсутствие статистически значимой зависимости показателя заболеваемости РЩЖ от дозы внешнего облучения.
Рак щитовидной железы, радиационный риск, ликвидаторы последствий аварии на Чернобыльской АЭС, Российский государственный медико-дозиметрический регистр, дозы внешнего облучения.
В результате аварии на Чернобыльской АЭС радиационному воздействию подверглись две основные группы: участники ликвидации последствий аварии, которые выполняли восстановительные работы в 30-километровой зоне Чернобыльской АЭС (ликвидаторы) и жители загрязненных радионуклидами территорий Белоруссии, России и Украины.
В последние годы в научной литературе опубликовано большое число работ, отражающих драматический рост заболеваемости раком щитовидной железы (РЩЖ) населения радиационно-загрязненных территорий [1-4]. Дается обоснование роли радиационного фактора в росте показателей заболеваемости РЩЖ за счет инкорпорированного облучения 131I, прежде всего для детского (на момент аварии) населения. Отмечается возможная роль и других факторов, обусловливающих рост заболеваемости: эффект скрининга (углубленная диспансеризация населения загрязненных территорий), йодный дефицит на этих территориях, генетическая предрасположенность, недостаточная йодная профилактика сразу же после аварии на ЧАЭС.
*Адресат для корреспонденции: 249036 Калужская обл., г. Обнинск, ул. Королева, 4, МРНЦ РАМН; тел.: (495) 95694-12, (48439) 9-33-90; факс: (495) 956-14-40; e-mail: nrer@obninsk.com.
Все указанные исследования относятся к жителям загрязненных территорий Белоруссии, России и Украины. Вместе с тем крайне ограниченное число работ посвящено проблеме изучения заболеваемости РЩЖ среди ликвидаторов [5, 6]. В наших предыдущих работах [7, 8] отмечалось повышение частоты заболеваемости РЩЖ среди ликвидаторов. В частности отмечалось, что увеличение частоты заболеваемости (в 2-3 раза) было выявлено в первые годы после аварии на ЧАЭС, что может быть объяснено повышенным уровнем диспансеризации ликвидаторов (эффект скрининга). В настоящей работе делается попытка оценки радиационных рисков в индукции заболеваемости РЩЖ среди ликвидаторов с использованием технологии внутреннего сравнения (internal comparison).
В настоящее время (на 01.01.2006 г.) в Российском государственном медико-дозиметрическом регистре (РГМДР) имеется информация на 638334 человек, подвергшихся радиационному воздействию в результате аварии на ЧАЭС и проживающих в России, в том числе на 187451 участника ликвидации ее последствий (ликвидатора).
Численность 50000 г 46258
40000
30000
20000
10000
0
1986
37111
20058
1987 Год въезда
1988-1990
Рис. 1. Распределение ликвидаторов по дате въезда ] 30-километровую зону.
Численность
40000г 35000 30000 25000 20000 15000 10000 5000 0
32228
33582
11625
7904
4279
п
9680
, 3^6, 963 ■
.и_I 11—11
<20 25-30 35-40 45-50 20-25 30-35 40-45 >50 Возраст на момент въезда
Рис. 2. Распределение ликвидаторов по возрасту на момент въезда в 30-километровую зону.
Численность
30000г 25000 20000 15000 10000 5000 0
27554
24915
- 14388
9403
- | | | 1
0-50 50-150 150-200 200+ Доза, мГр
Средняя доза, мГр 180 160 140 120 100 80 60 40 20 0
1986
168
- 93
- 33
- | | |
1987 Год въезда
1988-1990
Рис. 3. Распределение ликвидаторов по дозе внешнего облучения.
Рис. 4. Распределение средней дозы внешнего облучения по годам въезда в 30-километровую зону.
МАТЕРИАЛЫ И МЕТОДИКА
В настоящей работе рассматривается когорта ликвидаторов 103427 человек, проживающих в 6 регионах России: Северо-Западном, Волго-Вятском, Центрально-Черноземном, Поволжском, Северо-Кавказском и Уральском. Эти регионы отличаются повышенным качеством индивидуальной медицинской информации, поступающей ежегодно в РГМДР. Уровень охвата ликвидаторов ежегодной диспансеризацией на этих территориях составляет около 85%. Рассматривается заболеваемость РЩЖ указанной когорты ликвидаторов с 1986 по 2003 г.
На рис. 1 показано распределение ликвидаторов по дате въезда в 30-километровую зону: 1986, 1987, 1988-1990 гг. Как известно, восстановительные работы продолжались до 1991 г. В рассматриваемой когорте из 46258 ликвидаторов принимали участие в восстановительных работах в 1986 г., 22408 человека - в апреле-июле (когда было возможно инкорпорированное облучение щитовидной железы 1311) и 23850 человек - в августе-декабре.
На рис. 2 показано возрастное распределение изучаемой когорты ликвидаторов на момент их
работы в 30-километровой зоне Чернобыльской АЭС. Средний возраст ликвидаторов на этот момент времени составлял 34 года. Как видно из рис. 2, указанное возрастное распределение ликвидаторов дает принципиальную возможность для проведения долгосрочных эпидемиологических исследований.
На рис. 3, 4 показано распределение ликвидаторов по дозам внешнего облучения. Из рассматриваемой когорты ликвидаторов (103427 человека) 76260 имеют индивидуальные данные о полученных дозах внешнего облучения. Как видно из рис. 4, наибольшие дозы внешнего облучения получили ликвидаторы 1986 г. (средняя доза 168 мГр). Средняя доза ликвидаторов 1988-1990 гг. составила 33 мГр.
Важно отметить, что РГМДР располагает данными только о дозах внешнего облучения ликвидаторов. РГМДР не имеет информации о возможных дозах облучения щитовидной железы среди ликвидаторов за счет инкорпорированного облучения 1311. Поэтому оценка радиационных рисков проведена с использованием индивидуальных данных о дозах внешнего облучения.
Для нахождения величины стандартизованного отношения заболеваемости (SIR) в изучаемой когорте были использованы следующие модели:
X = exp (ft, )Xras, (1а)
X = exp 2 T12 )X rus, (16)
где Xrus - ежегодные (с 1986 по 2003 годы) повозрастные показатели заболеваемости РЩЖ взрослого мужского населения России в 5-летних возрастных интервалах; в0, 2 - неизвестные константы; Тх, 2 - категориальная переменная, введенная путем стратификации даты въезда на два уровня (первый: 1.08.1986 г. - 31.12.1990, второй: 26.04.1986 г. - 31.07.1986); SIR = exp(p0) в (1а) или SIRx, 2 = exp(P1, 2) в (16).
Эти две модели использовались для выявления гетерогенности SIR по дате въезда. В первом случае (1а) SIR не зависит от даты въезда в зону. Во втором случае, когда SIR зависит от даты въезда, оцениваются два значения SIR^ 2 по периоду въезда, SIRX с августа 1986 г. по декабрь 1990 г. и SIR2 с 26 апреля по конец июля 1986 г. В качестве статистического теста на гетерогенность использовался тест отношения правдопобий для моделей (1а) и (16).
Наиболее адекватная модель радиационного риска выбиралась посредством теста отношения правдоподобий при почленной редукции, начиная с членов взаимодействий факторов, в иерархическом семействе моделей следующего общего вида:
X = exp (aA + b,R, + c, (AR ))x
x( 1+ dD + e( DA) + f (DR,)),
где A - возраст на момент въезда; R, - административный регион (категориальная переменная, имевшая первоначально 6 уровней); D - доза внешнего облучения в Гр; (AR,), (DA), (DR,) - члены попарного взаимодействия между факторами; a, bj, c j, d, e, f - оцениваемые параметры, i = 1, ..., 6.
В процессе идентификации модели радиационного риска было выявлено, что члены взаимодействия являются избыточными, а первоначальные шесть административных регионов группируются в три региональные группы по уровню регистрируемой фоновой заболеваемости. Окончательно модель приобрела следующий вид:
X = exp(aA + bjRj )(1 + dD), (3)
где A - возраст на момент въезда; Rj - группы административных регионов по уровню фоновой заболеваемости (три уровня, j = 1, 2, 3); D - доза внешнего облучения в Гр; a, bj, d - оцениваемые параметры, среди которых d - величина избыточного относительного риска на единицу дозы (ERR/Гр).
Оценки всех параметров и тесты отношения правдоподобий вычислялись с использованием программного пакета EPICURE [9].
РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ
В модели (1а), без разделения когорты по периодам въезда, величина SIR оказалась статистически значимой и равной 3.39 с 95%-ным доверительным интервалом (ДИ) (2.73; 4.16). Это означает, что в среднем за все годы наблюдения регистрируемая заболеваемость ликвидаторов РЩЖ достоверно превышает ожидаемый спонтанный (национальный) уровень.
Кроме того, тест отношения правдоподобий для моделей (1а) и (16) показал, что величина SIR является гетерогенной по периодам начала работ ликвидаторов в зоне Чернобыльской АЭС: апрель-июль 1986 г. и август 1986 г. - конец 1990 г.
Наибольшая величина SIR выявлена для ликвидаторов, начавших участие в восстановительных работах в апреле-июле 1986 г. Среди этой группы ликвидаторов
Для дальнейшего прочтения статьи необходимо приобрести полный текст. Статьи высылаются в формате PDF на указанную при оплате почту. Время доставки составляет менее 10 минут. Стоимость одной статьи — 150 рублей.