CC BY
29
_ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА_
2019 Геология Том 18, № 4
ГЕОЛОГИЯ, ПОИСКИ И РАЗВЕДКА НЕФТЯНЫХ И ГАЗОВЫХ МЕСТОРОЖДЕНИЙ
УДК 551.734.5 (550.4)
Комплексный геолого-геохимический анализ отложений саргаевского горизонта Удмуртской Республики
И.С. Хоптаа' b, Д.Д. Кожановс
аПермский государственный национальный исследовательский университет 614990, Пермь, ул. Букирева, 15. E-mail: regional.psu@yandex.ru ЬАО «КамНИИКИГС»
614016, Пермь, ул. Краснофлотская, д.15. E-mail: KhoptaIS@yandex.ru с Московский государственный университет им. М.В. Ломоносова 119234, Москва, ГСП-1, Ленинские горы. E-mail: dimitriykozz@gmail.com (Статья поступила в редакцию 6 сентября2019)
Представлен анализ битумоидов саргаевского горизонта среднефранского подъяруса верхнего девона. Происхождение пород связано с накоплением осадков в относительно глубоководном морском бассейне, в условиях восстановительной и резко восстановительной среды. С целью генетического анализа проведена геохимическая интерпретация. Многочисленные следы перемещения свободных битумоидов свидетельствуют о миграционно-аккумуляционных трансформациях.
Ключевые слова: саргаевский горизонт, доманикиты, органическое вещество, битумоиды, Удмуртская Республика.
DOI: 10.17072/psu.geol.18.4.394 Введение
Отложения саргаевского горизонта Волго-Уральской нефтегазоносного басейна (ВУНБ) в связи с их высокой битуминозно-стью давно привлекают внимание исследователей. Начиная с ранних этапов изучения и по настоящее время проявляется особый интерес к высокоуглеродистым отложениям этого горизонта. Несмотря на то, что геологоразведочные работы на нефть в Удмуртии начались ещё в 1945 г., изученность её намного ниже, чем соседних с нею территорий Татарстана и Башкирии.
В настоящее время Удмуртия обеспечивает 9% объёма нефтедобычи в ПФО. Добывая 2% суммарной российской нефти, республика занимает 13-е место среди регионов страны. Степень освоения запасов - 47%.
На фоне того, что региональный фонд нефтегазоносности практически исчерпал свои возможности, следует прежде всего обратить внимание на необходимость выявления нетрадиционных залежей нефти в верх© Хопта И.С., Кожанов Д.Д., 2019
недевонско-турнейских отложениях Удмуртской Республики (УР). Для решения вопроса о новом источнике углеводородов (УВ) в продуктивных комплексах осадочного чехла следует использовать результаты новых геохимических исследований ОВ с использованием современных методов интерпретации.
Основные черты геологического строения Удмуртской Республики. Особенности осадконакопления отложений саргаевского горизонта ф^г)
Главной структурно-тектонической особенностью территории УР является наличие двух структурных мегакомплексов: консолидированного кристаллического фундамента архейско-раннепротерозойского возраста и верхнепротерозойско-кайнозойского осадочного чехла. В тектоническом строении выделены два крупных элемента: северная вершина Татарского свода (СТс) и Верхнекамская впадина (ВКВ). По строению отложений верхнедевонско-турнейского комплекса
СТс включает большую часть Можгинской палеовпадины Камско-Кинельской системы (ККСВ) и ее западный борт. В северную часть также попадают впадины Камско-Волжской системы (КВС). СТс в верхнеде-вонско-турнейском комплексе представлен Сарапульской впадиной, а на севере также впадинами КВС (Геология ..., 1976; Кожанов, Хопта, 2017).
В палеозое, по мнению Р.О. Хачатряна (Хачатрян и др., 1978), ВУНБ отличался сравнительно спокойным тектоническим режимом с преобладанием унаследованных структуроформирующих движений. В развитии палеозойских структур ведущая роль принадлежала механизму глыбовых опусканий с последовательным расширением территорий впадин и прогибов и сокращением сводов. Иными словами, палеозойский этап развития ВУНБ характеризовался постепенным погружением и расширением очагов генерации УВ и локализацией палеозон их концентрации.
Доманиковая фация УР приурочена к верхнедевонско-турнейским отложениям, которые представлены тремя типами разрезов: 1) впадинным; 2) сводовым; 3) переходным или бортовым.
Многие исследователи сходятся во мнении, что в саргаевское, доманиковое, мен-дымское и верхнефранское время в зоне ККВ скорость тектонических опусканий на территории развития собственно доманиковых фаций заметно превышала темп осадкона-копления, т.е. эти фации отражают условия некомпенсированных прогибаний (Можгин-ская и Сарапульская палеовпадины) (Кожанов, Хопта, 2017). Застойный режим морского бассейна, обогащение органическими остатками, восстановительная и резко восстановительная обстановка (сидеритно-сульфидные фации), частые смены трансгрессивного и регрессивного режимов седиментации привели к образованию кремнисто-битуминозных темно-серых известняков и известково-глинистых сланцев (нефтегазо-материнских толщ, перекрытых породами-покрышками). В процессе эволюции нефте-материнские породы ВУНБ не испытывали воздействия глубокого катагенеза (Баженова и др., 2012).
Активная тектоническая дифференциация доманикового бассейна в связи со специфическим сочетанием разнофациальных условий седиментации сопровождалась его па-леогеоморфологическим расчленением. При этом расположение структурно-фациальных зон позднедевонского бассейна в виде ККСВ рифового обрамления оказалось закономерно унаследованным относительно пале-оструктурных и палеогеоморфологических элементов доманикового бассейна. Иначе говоря, предпосылки резкого структурно-фациального расчленения позднедевонского бассейна были заложены в саргаевское время, явившееся началом развития качественно новых тектоно-седиментационных режимов в ВУНБ (Шеходанов, Армишев, 1970).
Таким образом, формирование нефтема-теринского потенциала пород саргаевского горизонта УР происходило в благоприятных условиях (трансгрессивный режим, накопление тонкопелитовых глинистых осадков (Павленко, 2015), оптимальная скорость осадконакопления, восстановительные условия, формирование ОВ за счет низших организмов планктона). Дифференцированный рельеф дна, существовавший в верхнедевонское время, способствовал созданию на западе активного седиментационного режима на месте унаследованных депрессий.
Результаты исследований
Геохимическая характеристика разреза саргаевского горизонта УР составлена по результатам интерпретации лабораторных исследований 242 образцов из интервала глубин 1656-2222 м (Нефти ..., 1977; Нефти ..., 1989).
Исследование битумоидов, полученных из карбонатно-глинистых отложений саргая верхнего девона на территории УР, позволяет выделить несколько групп гомологических рядов нормальных алканов. Каждая из них характеризует разность протекающих процессов внутри одного стратиграфического горизонта, покрывающего всю площадь УР, а также различие в условиях образования исходного ОВ пород.
Первая группа (рис. 1) по характеру конфигурации гомологического ряда (площади Дебесская, Сушинская, Балезинская, Есеней-
ская) определяется в области среднемолеку-лярных н-алканов (пики содержания в точках н-С19; н-С22; н-С23; н-С24), данные показатели характерны для смешанного, преимущественно сапропелевого исходного ОВ. Это подтверждается и наличием преобладания фитана над пристаном в общей структуре (Pr/Ph=0,47-0,59). Интересно также то, что отношения i-C19/n-C17; i-C20/n-C18 разнятся для битумоидов разных площадей. Так, у битумоидов Дебесской площади значения данного коэффициента равны 2,2 и 1,5, что характеризует наличие зоны генерации превращенных (преобразованных) нефтей. Остальные образцы менее подвержены процессам преобразования, коэффициент метаморфизма н-алканов в них ближе к 0,5, составляя в среднем 0,53 и 0,64 для каждой пары отношения. В составе н-парафинов мета-но-нафтеновой фракции помимо катагенно-генетического можно проследить влияние и седиментационного фактора. В первой группе отчетливо выявляется восстановительная и резко восстановительная обстановка мелководно-морского осадконакопления с застойным режимом, что достаточно точно подтверждает полученные ранее сведения о палеогеографии бассейна седиментации сар-гаевского времени франского века верхнего девона.
Вторая группа (рис.1) по результатам интерпретации образцов Лозолюкской и Наря-довской площадей также характеризуется общей приуроченностью ряда к среднемоле-кулярным гомологам (En-C17.23/En-C24.30 =2,5). Пики содержания н-алканов относятся к значениям содержания н-С19, н-С21, в среднем равные 14 и 12 % на сумму н-алканов. Отношение пристан/фитан повсеместно для образцов менее 1 (Pr/Ph=0,5-0,61). При наличии более высокой степени метаморфизма (i-C19/n-C17= 2,2 - 2,23; i-C2o/n-C18=1,28 - 1,3) можно с уверенностью заключить, что исходным являлось ОВ сапропелевой природы, накопление которого происходило в обстановках относительно глубоководного морского бассейна, в условиях восстановительной и резко восстановительной среды. Миграционная составляющая формирования битумоидов данных площадей, возможно, эпигенетична, т.к. центр гомологического ряда смещен по направлению
к средне- и даже низкомолекулярным гомологам в сравнении с членами ряда первой группы. Битумоид указанных площадей можно определить как миграционно-остаточный.
Третья конфигурационная группа ряда н-алканов битумоидов, полученных с Нарядов-ской и Аксеновцевской площадей (рис. 1), определяется повышенным содержанием более высокомолекулярных н-алканов в сравнении с предыдущими площадями (отношение En-C17.23/En-C24.30 близко к 1). Данный факт вместе с низким значением коэффициента метаморфизма свидетельствует о наличии здесь зон генерации превращенных и слабопревращенных нефтей (i-C19/n-C17= 1,17; i-C20/n-C18=1,49). Значение отношения н-алканов с нечетным числом атомов к н-алканам с четным примерно равно 1 (НЧ/Ч=0,88-1,03), пристан/фитан=0,39-0,35. В связи с этим битумоид имеет сапропелевую природу исходного вещества.
Характер распределения гомологического ряда н-алканов площадей Решетниковская и Восточно-Красногорская (рис. 1) определил наличие 4-й конфигурационной группы. В сравнении с предыдущими данная ассоциация образцов наиболее сильно обогащена высокомолекулярными н-алканами, пики содержания достигают значений 12,6-33,3 % на сумму н-парафинов в интервале от н-С26 до н-С29, НЧ/Ч=0,96. Эти особенности позволяют отнести исходное ОВ к отложениям прибрежно-морской относительно мелководной зоны осадочного бассейна, что отражает его смешанный гумусово-сапропелевый характер (Белоконь, Фрик, 1993). Однако отношение Pr/Ph=0,48 при условии того, что органическое вещество находилось в условиях генерации слабопревращенных нефтей (Милешина и др., 1983; Чахмахчев, 1983) (i-C19/n-C17=0,57; i-C20/n-C18=0,76), дает повод для сомнений в полученном выводе. По-видимому, пики гомологического ряда н-алканов в данном случае являются неинформативными в связи со своей слабой преобра-зованностью и отсутствием деструкции высокомолекулярных н-алканов до легких, и использование отношения Pr/Ph в реконструкции палеогеографических условий седиментации является наиболее целесообразным (Айзенштат, 2006).
Число атомов углерода в молекулах Число атомов углерода в молекулах
Рис. 1. Распределение п-алканов в битумоидах пород саргаевского горизонта по площадям
Так, по наблюдаемым значениям этого отношения необходимо сделать вывод о том, что в образовании битумоида существенную роль сыграло сапропелевое ОВ в условиях восстановительной и резко восстановительной среды относительно глубоководного морского бассейна.
Таким образом, в результате изучения характера распределения значений содержания н-алканов с разным числом атомов, а также рассмотрения основных характеристик, связанных с этими веществами, можно заключить, что исходное ОВ было образовано в условиях относительно глубоководного и мелководного застойного моря, в восстановительной и резко восстановительной среде. Прослеживается наличие как сингенетично-го, так и эпигенетичного битумоида, зон генерации слабопревращенных и превращенных нефтей. Различия в конфигурации гомологических рядов, по-видимому, обусловлены не спецификой исходного ОВ, а характером его преобразования в процессе диагенеза в совокупности с изменением тектонического плана ВУНБ в верхнефранское (саргаев-ское) время.
Был произведен анализ распространения основных показателей битумоида - Сорг, БХЛ, Бпэ, Бспб, Vop, Nip, Рхб - на территории УР как основных параметров при установлении нефтеносности отложений.
Органический углерод (Сорг) концентрируется как в пределах СТс, так и территории ВКВ. В основном данный показатель принимает максимальные значения в прибортовых зонах ККСВ, причем в областях, приуроченных к Сарапульской впадине, значения выше (12%), чем в областях, приуроченных к Мо-жгинской палеовпадине (8-9%) (рис. 2,А).
Характер распределения хлороформенно-го битумоида (Бхл) проявляется по-иному из-за того, что значительная часть УВ перемещенная. Бхл, представляющий собой приближенные по составу к нефтям фракции ОВ, на территории Удмуртской Республики концентрируется по большей части на территории ВКВ, чуть меньше - на северной вершине СТс (рис. 2,Б). Данная особенность, скорее всего, определяется повсеместным распространением восстановительной обстановки преобразования УВ и преобладанием ее над окисленной. Последняя преобладает в основном на границе ВКВ и СТс, а также близ границы с Бирской седловиной.
Рис. 2. Карта распространения: А - Сорг; Б - БХ
Подобным образом распределяется содержание гуминовых кислот (рис. 3,Б). При интерпретации следует учитывать переход гуминовых и фульвокислот в гумин в процессе диагенеза, который является предшественником керогена. То есть имеет место катагенетическая расчлененность пород саргаевского горизонта УР, при которой часть битумоида, достигшая начала главной фазы нефтеобразования, будет иметь низкие значения содержания гуминовых кислот. В местах, где битумоидный коэффициент понижен, а гуминовые кислоты имеют высокие значения, нефтегенераци-онный потенциал развит слабо вследствие недостаточности воздействия термобарического фактора. При этом зачастую возникает ситуация, когда оба значения достаточно сильно повышены, обычно это происходит с одновременным повышением Сорг. В таком случае можно сделать вывод о мощном нефтематеринском потенциале исследуемых пород, который в данное время активно реализуется, но частично (Айзенштат и др., 1989). То есть данная толща является максимально активной в плане генерации сланцевой нефти. К сожалению, по данному факту проведено мало исследований и
достоверно определить его на данном этапе можно лишь при наличии дополнительных факторов.
Большие значения рХБ саргаевского горизонта, характеризующие долю битумоида в ОВ и составляющие в среднем 40%, свидетельствуют о том, что в породах присутствует в основном паравтохтонный битумоид, согласно закономерности Успенского -Вассоевича. На территории УР наиболее крупные значения рХБ прослеживаются в восточной, юго-восточной и центральной части, преимущественно на территории Сара-пульской впадины. Увеличение рХБ в породах обусловлено нарастающим уровнем ка-тагенетического преобразования отложений и, следовательно, восстановленности биту-моидов. Широкое распространение парав-тахтонного битумоида связано с тем, что данные толщи представляют собой в основном закрытую систему, из которой эмиграция ОВ затруднена (рис. 3,А).
Одним из основных принципов геохимического анализа является использование биомаркеров (Vop, Nip) (Белоконь, Фрик, 1993). Данные химические структуры весьма широко распространены на территории
Удмуртской Республики в отложениях, приуроченных к саргаевскому горизонту. Площади их расположения не схожи друг с другом. Так, никелевые порфирины в основном концентрируются в пределах северной вершины СТс, тогда как ВКВ обогащена вана-диловыми комплексами. Здесь важную роль играют условия осадконакопления изучаемых отложений, а также специфика поведения этих компонентов в процессе миграции. По-видимому, обогащение никелевыми пор-фиринами происходило в условиях при-брежно-морского осадконакопления, развитого в саргаевское время в северной и центральной части Удмуртской Республики, а ванадиловыми - в относительно глубоководных морских условиях, распространенных в восточной и северо-восточной ее части. Общая обогащенность порфиринами, возможно, вызвана концентрированием их за счет потери легких фракций в процессе эмиграции.
Выводы
Таким образом, обоснованы перспективы нетегазоносности отложений саргаевского горизонта УР, проведена детальная характеристика индивидуального углеводородного
состава образцов битумоида отложений сар-гаевского горизонта.
Выяснено, что доманиковые отложения саргаевского горизонта в пределах УР являются нефтематеринскими (от средне - до высокопродуктивных) и нефтепроизводящими с высоким генерационным потенциалом. Многочисленные следы перемещения свободных битумоидов свидетельствуют о ми-грационно-аккумуляционных трансформациях, предопределяющих развитие залежей нефти в близлежащих проницаемых пластах (преимущественно в прибортовых зонах впадин), а также в отложениях самого горизонта. Благоприятные условия для нефтеоб-разования создавались в областях развития зарождавшейся ККСВ (по большей части на юге и в центральной части УР), где повсеместно зафиксированы высокие концентрации битумоидов.
Библиографический список
Баженова О.К., Бурлин Ю.К., Соколов Б.А., Хаин В.Е. Геология и геохимия нефти и газа: учебник для вузов. М.: Изд-во МГУ, 2012. С.78-179.
Белоконь Т.В., Фрик М.Г. Применение биомаркеров в нефтегазовой геологии. М.: Геоин-формарк, 1993. 47с.
Геология и нефтеносность Удмуртской АССР. Ижевск: Изд-во «Удмуртия», 1976. 128 с.
Кожанов Д.Д., Хопта И.С. Анализ распределения, условий накопления и генерации УВ в верхнедевонских отложениях Удмуртской Республики // Геология, геоэкология и ресурсный потенциал Урала и сопредельных территорий. Уфа: Альфа-Реклама, 2017. С. 292-297.
Милешина А.Г., Калинко М.К., Сафонова Г.Л. Изменение нефтей при фильтрации через породы. М.: Недра, 1983. 175 с.
Нефти, газы и битумоиды Пермского Прикамья и сопредельных районов: каталог физических свойств / под ред. С.А. Винниковского, А.З. Кобловой. Пермь, 1977. 568 с.
Нефти, газы и ОВ пород севера Урало-Поволжья: каталог физико-химических свойств / под ред. А.З. Кобловой, Ю.А. Дулепова/ Кам-НИИКИГС. Пермь, 1989. Вып. III. 685 с.
Павленко Т.В. Влияние повышенной битуми-нозности на петрофизические свойства и ее связь
с определенными литологическими типами (по данным исследования доманиковых отложений) // Вестник Пермского университета. Геология. 2015. №3, С. 75-83. DOI: 10.17072/psu.geol.28.75
Хачатрян Р.О., Алиев М.М., Батанова Г.П. Девонские отложения Волго-Уральской нефтегазоносной провинции. М.: Недра, 1978. 216 с.
Чахмахчев В.А. Геохимия процесса миграции углеводородных систем. М.: Недра, 1983. 231 с.
Шеходанов В.А., Армишев В.М. Геологическое строение и нефтеносность Камско-Кинельских прогибов в Удмуртской АССР // Геология и нефтеносность Камско-Кинельских прогибов / Казан. гос. ун-т. Казань, 1970. С. 49-57.
Aizenshtat Z. The Analysis of Alkanes and Cy-cloalkanes// The Chemistry of Alkanes and Cycloal-kanes. 2006. P. 289-349. D0I:10.1002/0470034378. ch7
Aizenshtat Z., Feinstein S., Miloslavski I., Ya-kubson Z. and Yakubson C.I. Oil-oil correlation and potential source rocks for oils in Paleozoic reservoir rocks in the Tataria and Perm basins, Russia// Organic Geochemistry. 1998. Vol. 29 (1-3 -3 pt 1). P. 701-712. D0I:10.1016/S0146-6380(98)00065-5
Comprehensive Geological and Geochemical Analysis of Sediments of the Sargaev Horizon of the Udmurt Republic
I.S. Khoptaa' b, D.D. Kozhanovc
a Perm State University
15 Bukireva Str., Perm 614990, Russia. E-mail: regional.psu@yandex.ru bJSC Kama Research Institute for Complex Deep and Supereep Borehole Investigation 15 Krasnoflotskaya Str., Perm 614016, Russia. E-mail:KhoptaIS@yandex.ru cLomonosov Moscow State University
GSP-1, 1 Leninskiye Gory, Moscow 119991, Russia. E-mail:dimitriykozz@gmail.com
The analysis of the bitumoids of the Sargaevskiy horizon of the Middle Frasnian sub-stage of the Upper Devonian is presented. The origin of the rocks is associated with the accumulation of sediments in a relatively deep sea basin, in a reducing and sharply reducing environment. For the purpose of genetic analysis, a geochemical interpretation has been carried out. Numerous traces of the movement of free bitumoids indicate migration-accumulation transformations.
Keywords: Sargaevskiy horizon; Domanicites; organic matter; bitumoids; Udmurt Republic.
Reference
Bazhenova O.K., Burlin Yu.K., Sokolov B.A., Khain V.E. 2012. Geologiya i geokhimiya nefti i gaza [Geology and geochemistry of oil and gas]. MGU, Moskva, pp.78-79. (in Russian)
Belokon T.V., Frik M.G. 1993. Primenenie bi-omarkerov v neftegazovoy geologii [The use of bi-omarkers in oil and gas geology]. Geoinformark, Moskva, p. 47. (in Russian)
Geologiya i neftenosnost Udmurtskoy ASSR [Geology and oil potential of the Udmurt ASSR]. Udmurtiya, Izhevsk, 1976, p.128 (in Russian).
Kozhanov D.D., Khopta I.S. 2017. Analiz raspre-deleniya, usloviy nakopleniya i generatsii UV v verkhnedevonskikh otlozheniyakh Udmurtskoy Respubliki [Analysis of the distribution, conditions of accumulation and generation of hydrocarbons in the Upper Devonian sediments of the Udmurt Republic]. In: Geologiya, geoekologiya i resursnyy potentsial Urala i sopredelnykh territoriy, AlfaReklama, Ufa, pp. 292-297 (in Russian).
Mileshina A.G., Kalinko M.K., Safonova G.L. 1983. Izmenenie neftey pri filtratsii cherez porody [Oil change during filtration through rocks]. Nedra, Moskva, p.175. (in Russian)
Nefti, gazy i bitumoidy Permskogo Prikamya i sopredelnykh rayonov. Katalog fizicheskikh svoystv. [Oil, gases and bitumoids of Perm Pri-kamye and adjacent areas. Physical Property Catalog]. S.A. Vinnikovskiy, A.Z. Koblova (Eds.), Perm, 1977, p. 568. (in Russian)
Nefti, gazy i OV porod severa Uralo-Povolzhya. Katalog fiziko-khimicheskikh svoystv. [Oil, gases and OM of rocks of the north of the Ural-Volga region. Catalog of physicochemical properties]. A.Z. Koblova, Yu.A. Dulepov, 1989, KamNIIKIGS, Perm, p. 685. (in Russian)
Pavlenko T.V. 2015. Vliyanie povyshennoy bi-tuminoznosti na petrofizicheskie svoystva i ee svyaz s opredelennymi litologicheskimi tipami (po dannym issledovaniya domanikovykh otlozheniy) [Influence of increased bituminosity on the petrophysical prop-
erties and its relationship with specific lithological types according to results of study of Domanik formations] . Vestnik Permskogo Universiteta. Geologiya. 3(28):75-83. (in Russian). doi: 10.17072/ psu.geol.28.75
Khachatryan R.O., Aliev M.M., Batanova G.P. 1978. Devonskie otlozheniya Volgo-Uralskoy nefte-gazonosnoy provintsii [Devonian sediments of the Volga-Ural oil and gas province]. Nedra, Moskva, p. 216. (in Russian)
Chakhmakhchev V.A. 1983. Geokhimiya protsessa migratsii uglevodorodnykh sistem [Geochemistry of the process of migration of hydrocarbon systems]. Nedra, Moskva, p. 231. (in Russian)
Shekhodanov V.A., Armishev V.M. 1970. Geolog-icheskoe stroenie i neftenosnost Kamsko-Kinelskikh progibov v Udmurtskoy ASSR [Geological structure and oil potential of the Kama-Kinelskiy troughs in the Udmurt Republic]. In: Geologiya i neftenosnost Kamsko-Kinelskikh progibov, Kazan University, Kazan, pp. 49-57. (in Russian)
Aizenshtat, Z. 2006. The Analysis of Alkanes and Cycloalkanes. The Chemistry of Alkanes and Cyclo-alkanes. pp. 289-349 doi: 10.1002/0470034378.ch7
Izenshtat Z., Feinstein S., Miloslavski I., Ya-kubson Z., Yakubson C.I., 1998. Oil-oil correlation and potential source rocks for oils in Paleozoic reservoir rocks in the Tataria and Perm basins, Russia. Organic Geochemistry, 29(1-3):701-712. doi: 10.1016/S0146-6380(98)00065-5
