Sunday, July 27, 2014

1 Критические уровни в процессах развития биологических систем

АКАДЕМИЯ НАУК СССР дальневосточный научный центр
институт биологии моря
A.B.ЖИРМУНСКИЙ, В.И.КУЗЬМИН
КРИТИЧЕСКИЕ
УРОВНИ
В ПРОЦЕССАХ
РАЗВИТИЯ
БИОЛОГИЧЕСКИХ
СИСТЕМ
ИЗДАТЕЛЬСТВО «НАУКА» МОСКВА 1982












Жирмунский а. В., Кузьмин В. и. Критические уровни в процессах развития биологических систем. м.: Наука, {982,180 с.
Монография, написанная членом-корреспондентом АН СССР А. В. Жирмунским и доктором технических наук В. И. Кузьминым, посвящена рассмотрению различных биологических систем с точки зрения соотношения эволюционных и революционных (критических) периодов в их развитии. На основе анализа аллометрических зависимостей развития систем авторами выявлено соотношение между последовательными критическими значениями аргумента, при котором система скачком переходит в новое устойчивое состояние. Проверка этого соотношения, проведенная в ряде работ с участием специалистов по эмбриологии, популяционной биологии и геологии, показала его действенность. Сделано заключение, что иерархически связанные процессы развития систем имеют среди критических уровней такие, соотношение последовательных значений которых равно ее = 15, 15....
Рассчитана на биологов широкого профиля, экологов, эволюционистов, эмбриологов, палеонтологов, геологов, кибернетиков, биофизиков, специалистов, интересующихся применением теории систем в биологии.
Табл. 24, ил. 77, библ. 359 назв.
-Ответственный редактор доктор биологических наук А. В. ЯБЛОКОВ
О   Издательство «Наука», 1982

ПРЕДИСЛОВИЕ
Предлагаемая книга биолога А. В. Жирмунского и специалиста по моделированию систем В. И. Кузьмина посвящена рассмотрению процессов развития биологических систем на разных уровнях организации живой материи. Рассматриваемая проблема является центральной как в биологии, так и вообще в естествознании. Развитие природных систем, с одной стороны, характеризуется большим разнообразием специфических свойств и механизмов их реализации, ас другой — некими общими тенденциями, которые уже давно обнаружены в философии и естествознании. Данная книга посвящена количественному анализу общих закономерностей процессов развития.
G давних времен в естествознании имеются противоборствующие направления, одно из которых отдает предпочтение эволюционным тенденциям в процессах развития, а другое обращает внимание на систематические революционные преобразования («катастрофы»). В последнее время намечается объединение этих точек зрения, при котором эволюция и революция принимаются как последовательные фазы развития, связанные с переходом количественных изменений в качественные.
До недавних пор исследования процессов развития носили описательный характер либо основывались на эволюционных моделях. Вместе с тем были отдельные попытки, не получившие завершения, использовать непрерывные модели на отдельных этапах развития систем. Границы применимости модели определялись по анализу фактического морфофункционального материала. При этом оставались открытыми вопросы о том, насколько специфичны моменты наступления новой фазы развития системы и существуют ли некоторые общие закономерности смены фаз в развитии. Исследованию именно этой проблемы посвящена предлагаемая читателю книга.
Особенностью подхода авторов к проблеме является учет релаксационных (запаздывающих) реакций на формирование текущих характеристик процесса развития. В книге рассматриваются факторы и закономерности формирования этапов развития, далее

вводится математическая модель, результаты которой проверяются на большом материале.
Основной результат, который представляют на обсуждение авторы, состоит в определении некоторых постоянных соотношений между характеристиками развивающихся систем в момент смены характера развития. Если авторы правы, то вывод имеет чрезвычайно большое теоретическое и практическое значение. Например,! знание количественных значений критических численностей могло бы стать теоретической базой для проведения ряда мероприятий, связанных с охраной природы и рациональной эксплуатацией биологических ресурсов. Выведенные закономерности7 могут представить существенный интерес также для геологии, в частности при построении единой геохронологической шкалы.
Наряду с дальнейшей проверкой полученных результатов в % биологии существенный интерес, по нашему мнению, представляет анализ аналогичных явлений, выходящих за пределы биологии и геологии.
Академик М. С. Гиляров Академик Б. С. Соколов

ВВЕДЕНИЕ
Процессы развития организмов, их сообществ и всей биосферы Земли постоянно привлекают внимание биологов. Каковы закономерности развития этих и других биологических систем во времени и в пространстве, что общего между функционированием и развитием биологических систем разных уровней организации и интеграции и чем они принципиально различаются, как соотносятся медленные (эволюционные) и быстрые (революционные, скачкообразные) этапы развития? Эти проблемы представляют большой теоретический интерес и решение их может иметь весьма важное значение для практики.
Анализируя рост морских моллюсков и других животных (по литературным данным), авторы обратили внимание на изменение режимов развития животных и начали изучать их соотношения. При этом была обнаружена закономерная смена периодов алло-метрического роста более короткими по времени перестройками. Использование представлений о запаздывающих процессах позволило нам вывести математические зависимости и предложить модель смены эволюционных аллометрических этапов развития критическими периодами.
При моделировании роста используют два основных подхода. Первый и наиболее часто применяемый состоит в попытках ввести единую кривую роста, описывающую весь период развития [Вин-берг и др., 1968; Bertalanffy, 1968; и др.].
Второй подход связан с выделением в развитии участков, на которых оно происходит в рамках постоянных (или близких к постоянным) экзогенных и эндогенных условий [Brody, 1927; Шмальгаузен, 1935]. Остроумов [1912а,б, 1918] отмечал, что рост неоднороден на протяжении жизни организма и должен быть расчленен на отдельные этапы, на каждом из которых определяющую роль играют различные факторы. Детально возможности этого направления были исследованы Шмальгаузеном (1935) в рамках модели параболического роста. В этом случае процесс роста описывается степенной зависимостью
у=аТ\
где у — размер системы, Т — ее возраст, а и Ъ — константы. Величина константы Ъ показывает, на сколько процентов изменится функция у при изменении аргумента на 1 %.
При изменении лимитирующих рост факторов происходит, как показал Шмальгаузен (1935), скачкообразное изменение величины константы Ь. Примером такой смены величины константы Ъ
5

2,0 -|
1,0 -Щ
Ä
•|   l  1   i  i   i   i  i   i   i  i  i   111
ill   i   if
0
5 9 12
17
21 Годы
Рис. 1. Изменение интенсивности роста мальчиков с возрастом по абсциссе — возраст после рождения; по ординате — константа роста — показатель степени аллометрической зависимости б, характеризующий процент прироста массы на 1% прироста возраста. Кг — К* — критические возраста, в которых изменяется константа роста [по: шмальгаузен, 1935]
является рассмотренный Шмальгаузеном (1935) процесс роста массы тела человека (рис. 1). Видно, что величины относительных приростов массы на единицу относительного прироста возраста (Ь) остаются постоянными на отдельных участках, а в момент рождения и в возрасте 1, 12 и 19 лет скачком меняются.
Именно в этих возрастах в организме человека происходят изменения, приводящие к коренным перестройкам обменных процессов. Это: переход от эмбрионального периода, когда питание происходит за счет веществ, поступающих от материнского организма, к постэмбриональному лактационному периоду, затем переход к самостоятельному питанию, стимуляция роста, связанная с половым созревацием, и прекращение роста.
Критические состояния в развитии системы (перестройки) можно продемонстрировать на кривой зависимости объемов яйцеклеток морских ежей Strongylocentrotus nudus от возраста клеток (рис. 2). Видны точки излома кривой роста (точки 1—4); будем
too
h
Рис. 2. Изменение объема яйцеклеток морского ежа Strongylocentrotus nudus с возрастом
10
по абсциссе — возраст (месяцы годичного цикла); по ординате — объемы яйцеклеток (шкала логарифмическая), 1—4 — критические точки; tx — t4 — критические возраста; Х\ — х4 — критические уровни. развитие , между точками 2 и 3 — критический период (фаза пег рестройки) [по: гнездилова и др., 1976]
tf ti tj ц
i  i м i i i i  i i i i
IX  XI I III  V  У11
Месяцы
6

называть их критическими точками, Сортветствующие им возраста (h—h) назовем критическими возрастами, или критическими аргументами, а соответствующие значения функций (#4—хк) — критическими уровнями. Переходную стадию между двумя последовательными периодами роста размеров яйцеклеток (системы), на которой эти размеры убывают, назовем критическим периодом (фазой), или фазой перестройки. Как видно из рис. 2, две последовательные фазы роста яйцеклеток 1—2 и 3—4 прерываются фазой перестройки 2—3.
Таким образом, возможность выделения по морфофункциональ-ным признакам фаз развития, на которых поддерживается определенный механизм регуляций, определяющий характер роста, очевидна. Открытым является вопрос о возможности определения моментов изменения тенденций роста непосредственно по результатам наблюдений над количественными характеристиками растущей системы задолго до наступления критического периода. Решение именно этой проблемы и составило предмет исследований, результаты которых излагаются в книге. Одной из задач проведенных нами. исследований было моделирование моментов наступления критических периодов в развитии биологических систем, исходя из количественной информации о предыстории их развития.
Далее были произведены проверки предложенной модели на биологических системах, различающихся по уровню интеграции и типу развития. Это потребовало использования и осмысливания материала из разных разделов биологии и даже геологии, в которых авторы, естественно, недостаточно компетентны. Поэтому отдельные работы, опубликованные в ряде статей, выполнены в сотрудничестве со специалистами — эмбриологом к.б.н. С. Г. Васец-ким, популяционным биологом д.б.н. А. В. Яблоковым, палеонтологами и геологами акад. Б. С. Соколовым и чл.-корр. АН СССР В. Д. Наливкиным. Авторы выражают им глубокую благодарность за активное и увлеченное участие в рассмотрении имеющихся данных с точки зрения нетрадиционных представлений и за согласие на включение материалов совместных работ в настоящую книгу.
Материалы, вошедшие в книгу, обсуждались на ученых советах и научных семинарах Института биологии моря и Тихоокеанского института географии ДВНЦ АН СССР (Владивосток), Института биологии развития АН СССР (Москва), Всесоюзного нефтяного научно-исследовательского геолого-разведочного института (Ленинград), а также с академиками А. П. Александровым, Г. И. Марчуком, А. А. Баевым, М. С. Гиляровым, А. А. Дородницыным, Н. Н. Красовским, В. А. Коптюгом, Е. М. Крепсом, Ю. А. Овчинниковым, М. А. Садовским, В. Е. Соколовым, членами-корреспондентами АН СССР Г. Г. Винбергом, И. И. Гительзоном и Д. Г. Кнорре, д.ф.-м.н. Ю. И. Журавлевым и с Г. Т. Вороновым. Авторы приносят перечисленным ученым и всем, принимавшим участие в обсуждении, искреннюю признательность за внимание, поддержку,  ценные  критические  замечания  и рекомендации.
7

1. ПРОБЛЕМА КРИТИЧЕСКИХ УРОВНЕЙ В РАЗВИТИИ СИСТЕМ
Под развитием понимается направленное, закономерное изменение материи и сознания; в результате развития возникает новое качественное состояние объекта — его состава или структуры. «Различают две формы развития, между которыми существует диалектическая связь: эволюционную, связанную с постепенными количественными изменениями объекта, и революционную, характеризующуюся качественными изменениями в структуре объекта» [Сов. энциклопед. словарь. М.: Сов. энциклопедия, 1980, с. 1109].
К такому пониманию сущности и состава процессов развития человеческая мысль пришла в результате длительной борьбы идей. Вместе с тем элементы таких воззрений известны с древних времен. Так, в древнем Вавилоне наряду с представлением о сотворении мира божеством и его неизменности было распространено учение о циклическом развитии природы, связанном с изменением положения небесных тел, о «великом годе», по прошествии кото- j рого на Земле должны повториться те же события, которые происходили в начале этого «великого года». i
Представления о стадиях развития вселенной и Земли, связанных с «великим годом», развиваются далее древнегреческими; философами Гераклитом и Эмпедоклом. Энгельс отмечает, что \ «первоначальный, наивный, но по сути дела правильный взгляд на мир был присущ древнегреческой философии и впервые ясно выражен Гераклитом: все существует и в то же время не существует, так как все течет, все постоянно изменяется, все находится в постоянном процессе возникновения и исчезновения» [Маркс, Энгельс. Поли. собр. соч. 2-е изд., т. 20, с. 20]. Энгельс отмечает далее, что диалектика «была исследована более или менее точным .< образом лишь двумя мыслителями: Аристотелем и Гегелем» (там же, с. 367), и, несмотря на идеалистические посылки, «именно Гегель первый дал всеобъемлющее и сознательное изображение ее всеобщих форм движения» (там же, с. 371). Маркс и Энгельс соединили диалектику «как самое всестороннее, богатое содержанием и глубокое учение о развитии» (Ленин. Поли. собр. соч., т. 26, с. 53) с материалистическим воззрением на природу и общество, обосновав новое мировоззрение — диалектический материализм.
Хотя к диалектическому пониманию природы можно прийти на основе накапливающихся фактов естествознания, но,—как указывал Энгельс —«его можно легче достигнуть, если к диалектическому характеру этих [естественнонаучных] фактов подойти с пониманием  законов   диалектического   мышления»   (Маркс   и
8

Энгельс. Поли. собр. соч. 2-е изд., т. 20, с. 14). История развития идей в естествознании, и в частности в геологии, палеонтологии и биологии, полностью подтверждает этот тезис Именно поэтому авторы сочли необходимым для обсуждения проблемы развития природных систем наряду с данными математики, биологии, палеонтологии и геологии рассмотреть также кратко историю борьбы воззрений в этой области в перечисленных выше конкретных науках в свете представлений диалектико-материалистической философии.
В разделе 1.1 обсуждаются вопросы о биологических системах, иерархических уровнях их интеграции и общей теории систем Берталанфи. В разделе 1.2 излагается кратко история взглядов на соотношение эволюции и скачков в развитии. Раздел 1.3 посвящен месту, которое занимает проблема критических уровней в изучении развития сложных систем.
1.1. Биологические системы и иерархические уровни их интеграции
В последние годы в философии и конкретных науках, в частности биологии, все большее значение уделяется представлениям о системах. Под системой понимается «множество элементов, находяг щихся в^отношениях и связях друг с другом, которое образует определенную целостность, единство» [Садовский, 1976, с. 463]. Система имеет определенную структуру, т. е. взаимодействие между элементами влияет на элементы, а элементы, в свою очередь, оказывают влияние на систему. Система взаимодействует с другими системами, в том числе с системами более высокого порядка, являющимися для нее внешней средой. Таким образом, системы характеризуются иерархичностью, многоуровневостью: они состоят из элементов, которые, как правило, сами являются системами, и вместе с тем являются элементами систем более высокого порядка.
Еще в прошлом веке в связи с созданием клеточной теории было осознано наличие различных ступеней, или уровней организованности живой материи [Schleiden, 1846; Вирхов, 1865; Бер-нар, 1878; и др.]. В первой половине XX века эта идея получила дальнейшее развитие как среди философов [Brown, 1926; Sellars, 1926], так и среди биологов. Обсуждались разные уровни организации биологических систем, их функционирования и интеграции, а также подходы к исследованиям на различных уровнях.
В лекциях, читанных в 1912—1913 гг. в Военно-Медицинской академии, И. П. Павлов говорил о трех этажах, на которых должны проводиться исследования в физиологии, образно уподобляя последнюю зданию: в верхнем этаже осуществляется «анализ взаимоотношений и высочайшего уравновешивания организма и окружающего мира» [И. П. Павлов, 1952а, с. 58], средний этаж — это системная и органная физиология, наиболее развитая в то
9

время, а нижний этаж.-*- физиология клетки, которая составляв «дно, фундамент жизни» [И. П. Павлов, 1952а, с. 65]. В друго выступлении Павлов пишет о немецком физиологе Гейденгайне -~; представителе «той физиологии, которая должна сменить нашу современную органную физиологию и. которую можно считать предвестником последней ступени в науке о жизни — физиологии живой молекулы» [Й. П. Павлов, 1949, с. 162].
Таким образом, Павлов выделяет четыре уровня, подлежащих физиологическому исследованию,— физиология целостного организма, органная физиология, физиология клетки и физиология живой молекулы.
В середине XX века рядом биологов, социологов и философов была разработана теория интегративных уровней организации. Основные положения этой теории сформулированы в статье биолога Новикова [Novikoff, 1945]: «Концепция интегративных уровней организации — это общее описание, эволюции материи через последовательные и все более высокие порядки сложности и интеграции. Она рассматривает развитие материи, начиная от космологических изменений, приводящих к образованию Земли, до социальных изменений в обществе, как непрерывное, потому что оно никогда не кончается, и как прерывистое, потому что оно проходит через ряд разных уровней организации — физический, Химический, биологический и социологический. ...Эта концепция... не сводит явления высшего уровня к явлениям низшего, как в механизме, ни описывает явления более высокого уровня в туманных нематериалистических терминах, заменяющих понимание,! как это делает витализм» [с. 209]. Новиков различает следующие уровни организации живой материи: клетки, ткани, органы, системы органов, организмы, популяции.
Весьма большую роль в развитии представлений о наличии в живой природе пространственной иерархии взаимосвязанных систем сыграл Берталанфи [Bertalanffy, 1952; и др.]. Предложенную им классификацию иерархических порядков «живого уровня» [Bertalanffy, 1949] Сетров [1971, с. 41] коротко выражает следующим рядом: поле, плазмаэлементарные частицыатомы -^молекулы-^агрегаты молекул, мицеллы-»-коллоиды-*- клетки -^ткани-^органы->организмы->колонии, сообщества живых организмов.
Значительный интерес к проблеме уровней проявляют биологи разных специальностей: зоологи Гиляров [1954], Беклемишев [1964], Маркевич [1968]; ботаник Лавренко [1965]; физиолог Баллок [Bullock, 1955]; медики Давыдовский [1958], Шумаков и: др. [1971]; цитофизиологи Александров [1975], Ушаков [1963]; биохимики Сцент-Дьердьи [1964], Поглазов [1970], Энгельгардт [1970]; цитологи де Робертис и др. [1973]; морской биолог Зенкевич [1967]; эволюционисты Шмальгаузен [1961], Тимофеев-Ресовский и др. [1977]; экологи Дайс [Dice, 1952], Наумов [1964] и Одум [1975] и другие. Они выделяют различные ряды структурных*    функциональных,    организационных,    интегративных
*0

уровней «живого», «жизни» или живой материи * особенно деталь* но обсуждая области своих интересов.
Не меньшее внимание этой области уделяют философы. Выпускаются специальные монографии [например, Кремянский, 1969; Сетров, 1971], сборники [Проблема целостности в современной биологии, 1968; Проблема уровней и систем в научном познании, 1970; Развитие концепции структурных уровней в биологии, 1972; и др.], создаются сложные схемы, авторы которых пытаются объединить все известные уровни в единую систему [Веденев и др., 1972, с. 67; Сетров, 1972, с. 315], развертываются дискуссии о числе уровней «живого», начальных и конечных уровнях, количествах и типах иерархий систем.
f}zgg^^ RnTTP0f- п ттжте1? грата-
це жизни. Можно ли говорить о «физиологии живой молекулы», как это" делал И. П. Павлов [1949, с. 162], выделять молекулярную [Энгельгардт, 1970], физико-химическую [Овчинников, 1980] и даже субмолекулярную [Сцент-Дьердьи, 1964] биологию, определять жизнь, как «самореализацию потенциальных возможностей электронных состояний атомов» [Бернал, 1969, с. 17]. Интересную пресс-конференцию по этому вопросу провела редакция научно-популярного журнала «Наука и жизнь» (1962, № 4, с. 7—13) — она, показал а значительное разнообразие точек зрения известных советских биологов по вопросу о том, являются ли живыми вирусы и белковые молекулы.
Не менее острые разногласия проявились в области экологии.
Возникшее на стыке цитологии и экологии направление исследований, получившее название гщтоэ^^лоти [Александров, 1952, 1975; Ушаков, 1963; Библь, 1965ГяШрм^ский, 1966; и др.], изучающее клеточные адаптации и клеточные механизмы надклеточ-ных адаптации, по мнению некоторых ученых, не может быть отнесено к экологии [Зенкевич, 1967; Гиляров, 1973], так как она, с их точки зрения, должна исследовать лишь организменные и надорганизменные приспособления.
По мнению Макфедьена [1965], главный объект экодагии--популящи, другие исследователи ставятТво главу ^гла^^одш^ы, биогеоценозы [Сукачев, 1957] или дя$осистемы [Tansley,^l935], Наконец, в"связи с привлечением внимания if взаимоотношениям человека и биосферы, многочисленные авторы, в подавляющем числе небиологи, стали говорить и писать о «эдобальной экологии» [см., например, Будыко, 1977]. "
Другая группа спорных вопросов возникла в связи с осмысливанием мно^с^венности путе^до к^
,е^§§?15ё1^Кляются лйТколонии одноклеточных многоклеточными организмами? Что такое инфузория — клетка или организм? Следует ли относить популяции и виды к разным или единому «по-пуляционно-видовому уровню»?
Ушаков [1963], обсуждая роль различных уровней в процессах адаптации и эволюции, предложил различать следующие И координационных уровней: молекулярный, органоидный, клеточный,

органно-тканевой, системы органов (корреляты), организменный, семейно-стадный (репродукционный), популяционный, видовой и межвидовой (ценотический). Последний уровень Ушакова обозначен Лавренко [1965] как биоценотический. Им также выделен биостроматический уровень, соответствующий «слою» всего живого вещества на нашей планете.
Завадский [1966] подверг схемы Ушакова [1963] и Наумова [1964] критике, находя, что их подход к классификации уровней формален, предложенный ими ряд «линейной последовательности» — «сборный» и «искусственный» и что в их иерархиях «сосуществуют принципиально разнородные системы — первичные и универсальные (организм, вид, биоценоз, биосфера) со вторичными и частными (орган, ткань, семья, стадо и т. п.)» (с. 41). Он предложил свою схему: «Основные формы организации живого, главные ступени их эволюции и структурные компоненты каждой ступени» (с. 42). Однако эта схема не только не пригодна для использования ввиду чрезвычайной сложности определений орга-низменной и «популяционно-видовой» ступеней «первичных систем», но и потому, что реально существующие, известные и доступные исследованию уровни организации живой материи в этой системе исчезают.
Основываясь на представлениях В. И. Вернадского [Вернадский, 1978] об одновременном появлении сложного комплекса форм (организмов, которые одновременно были видами, составляли биоценозы и в совокупности биосферу), Завадский выкидывает из своей иерархии «основных, первичных форм организации живого» клетку — этот основной кирпич, из которого построена живая природа1.
Обвиняя биологов в «искусственной» линейной последовательности, он сам «впадает в противоречия», не приняв во внимание соотношения различных систем иерархии, разработанные Вудже-ром [Woodger, 1937] и Берталанфи [Bertalanffy, 1960].
По Берталанфи можно выделить несколько иерархических рядов систем, а именно: иерархия деления, например, при размножении тканевых клеток или простейших; пространственная иерархия многоклеточного организма, состоящего из органов, тканей, клеток, ядер; генетическая иерархия (например, филогенетическая) — соотношение организмов, относящихся к виду, с организмами различных родов, семейств и т. д.; физиологическая иерархия процессов, примером которой может служить осуществление функции нервной системы на разных уровнях — интимные процессы, протекающие в нервных клетках и волокнах, которые лежат в основе реализации рефлекторных актов; последние в свою
1 Кстати сказать, по представлениям Опарина [1957; и др.] и большинства современных биологов, жизнь на Земле возникла из гетерогенных коллоидных частиц, коацерватов, которые не имели такой сложной структуры, как современные одноклеточные организмы или тканевые клетки, но уж если их можно считать организмами, то, конечно, организмами одноклеточными (или, как часто говорят, доклеточными формами).
12

очередь регулируются деятельностью центральной нервной системы.
Как отмечает Сетров [1966, 1972], важнейшими для биологии необходимо признать пространственную и генетическую иерархии. Прежде всего из-за их столкновения возникают противоречия при классификации биологических систем.
Как указывалось выше, в схеме Завадского отсутствует уровень клеток, которые он обозначает «структурным компонентом» организма. Действительно, клетки — структурный компонент многоклеточного организма, но и в этой функции они являются полноправными членами пространственной иерархии — ряда, в который входят и клетка, и организм. Вместе с тем клетки простейших, микроорганизмов, половые клетки организмов с внешним оплодотворением (многие морские беспозвоночные), зиготы несут двойственную функцию — клетки и организма [Полянский, 1959, 1971]. В частности, они непосредственно сталкиваются с действием факторов внешней среды и потому, как отмечает Ушаков [1963], подвержены действию отбора.
Сетров [1972] отмечает, что подобная же ошибка допускается Завадским [1961, 1968] при объединении популяции и вида в единый «популяционно-видовой уровень». По всей видимости, такое объединение производится рядом исследователей потому, что изредка вид бывает представлен единственной географической популяцией. Кроме того, в некоторых случаях особи популяции, относящиеся к одному виду, длительно изолированы и, казалось бы, не имеют никаких связей друг с другом. Однако, отмечает Сетров, при этом не учитывается, что связь может быть не только пространственной, но и временной. Популяции, длительное время разделенные барьерами, сохраняют свои наследственные признаки, характерные для вида в целом, которые они получили от общих предков.
Убедительным подтверждением этой точки зрения являются данные Ушакова [1959], который установил, что теплоустойчивость мышц (генотипический видовой признак) морских ежей Strongylocentrotus intermedins из изолированных друг от друга в течение нескольких тысяч лет популяций из прибрежных вод Курильских островов, Сахалина и Приморья является одинаковой.
Другим примером генетической связи между поколениями могут служить опыты Жирмунского и др. [1967], показавших методами изотопной термометрии, что верхние температуры обитание ископаемых и современных моллюсков одного и того же вида остаются теми же самыми.
В работе Калабушкина [1976] установлена высокая стабильность популяционной системы брюхоногих моллюсков Littorina squalida лагуны Буссе [о-в Сахалин], выражающаяся в сохранении средней частоты генов окраски раковин на протяжении 2000 поколений.
Сетров [1971] приходит к выводу, что отрицание правомерности построения «линейных» рядов иерархии не может быть
13

признано правильным, так как эти иерархии реально существ вуют; Схема организационных уровней/предложенная Сетровым [1972, с. 315], представляется наиболее удачной из известных нам «всеобъемлющих» схем. Она объединяет живую материю с неорганической природой, учитывает существование различных видов* иерархий и двойственную природу клеток и популяций, являющихся одновременно системами разных иерархий. Естественно, что она может быть подвергнута критике разными специалистами» в том что касается названия и числа конкретных уровней, но это связано прежде всего с недостаточной разработанностью этих разделов биологии.
Следует отметить, что некоторые авторы [Завадский, 1966J предлагают сравнительно небольшое число «основных» уровней. У Одума (1975) это 6 «главных уровней организации жизни»: ген, клетка, орган, организм, популяция, сообщество2. В «Большой Советской Энциклопедий» [Астауров и др., 1970] даются другие 6 уровней: молекулярный, клеточный, организменный, популяци-онно-видовой, биоценотический, биосферный.
В 4 разделе рассматриваются размерные диапазоны различных структур организма, которые должны быть использованы при построении новой схемы структурно-функциональных уровней организации живой материи, исходящей как из функционального значения различных структур, так и из их линейных размеров.
Развитие представлений об организованности материи, ее структурированности, выражающейся в существовании иерархической системы взаимосвязанных систем разных уровней интеграции, привело к поискам общих и частных закономерностей, которые характеризуют эти системы.
Этой проблеме посвящены работы Богданова (1912, 1925—29), в которых разработана применимая к любым явлениям «всеобщая организационная наука» — тектология, основным принципом которой является математическое моделирование, ^ктрдощя, рассматривает явления с точки зрения непрерывно протекающих^ процессов организации и дезорганизации и выявления регулдрую-йдаг^йГэд^ Для живой протоплазмы Богданов
Г1925*—11529Х считает характерным подвижное равновесие, которое регулируется специфическим «бирегулятором».
Аналогичные представления позже развивали Бауэр [1935], которой разрабатывал теорию (^устойчивого неравновесия», и Насонов и Александров [1940], выдвин^шйе концепцию «подаиж-mrroj^ живой протоплазмы, в основе которого лежат
денатурационные   и   ренативационные   изменения   клеточных
Бауэр [1935] полагал задачей теоретической биологии «объединить и выразить в форме одного или нескольких законов то, что
2 Заметим, что первый член этого ряда относится в соответствии с классификацией Берталанфи к генетической иерархии.
14

свойственно всем системам, которые мы обозначаем исак живые, и что характерно только для них» (с. 22).
Сходную задачу ставит перед биологией автор общей теории систем ^-Берталанфи [1952]. Он определяет ее как «логико-материалистическую область, предметом исследования которой является формулирование и выведение тех принципов, которые имеют силу для систем вообще» (с. 199). Берталанфи [1960, 1968] полагает, что имеются общие принципы, поддерживающие существование любых систем, независимо от природы составляющих элементов и отношений между ними. Эти принципы могут быть определены математическим языком. Нахождение этих принципов и математических закономерностей и составляет главную задачу общей теории систем. Они позволят поставить на точную научную основу такие основные представления биологии, как целостность, динамика, организация.
1.2. О соотношении эволюции и скачков в развитии
Еще в древнем мире наряду с воззрениями о неизменности мира были распространены представления о его развитии, постепенном или неравномерном. Возрождение материалистической философии в XVII—XVIII веках привело к осознанию необходимости развития науки на основе опытных знаний [Ф. Бэкон] и способствовало развитию естественных наук. Однако господство теологических догм над естествоиспытателями накладывало отпечаток на развиваемые ими представления. В конце XVIII — начале XIX века последовательно формируются три главных направления в объяснении основных законов развития Земли и живой природы — катастрофам, униформизм, эволюционизм [Уэвелл, 1869; Т. Н. Huxley, 1883; Равикович, 1969; Завадский, Колчинский, 1977; и др.].
Хотя идеи катастрофизма высказывались в сочинениях отдельных натуралистов XVIII века, наиболее яркое его обоснование было дано Ж. Кювье в 1812 г. в работе «Рассуждения о переворотах на поверхности земного шара» [Кювье, 1937]. Кювье выдвинул гипотезу о периодически повторяющихся катастрофах не глобального характера, после которых от уцелевших организмов происходили новые сообщества. Через два десятилетия Лайель [1830—1833] опубликовал книгу «Основные начала геологии, илд новейшие изменения Земли и ее обитателей», которая явилась основой униформизма [Лайель, 1866]. Лайель приводил доказательства постепенного преобразования земной поверхности и условий жизни на ней. Взгляды Лайеля оказали огромное влияние на развитие естествознания. Как отмечал Энгельс [Маркс, Энгельс. Поли. собр. соч. 2-е изд., т. 20, с. 352], теория Лайеля «приводила непосредственно к учению о постепенном преобразовании организмов... к учению об изменчивости видов».
Представление об эволюции живых существ от простых форм к все более сложным под влиянием изменения внешних условий было подготовлено трудами ряда натуралистов, среди которых

выдающееся место занимал Ламарк. Однако ввиду противоречивости в объяснении причин эволюции и недостаточной обоснованности его теория не оказала революционизирующего влияния на развитие современной ему науки. И лишь опубликование «Происхождения видов» Дарвина [1859] привело к победе эволюционной теории в биологии.
Однако в XX веке отвергнутые было теории обретают на новой основе новых сторонников. В то время как А. Н. Северцов, Шмальгаузен, Добжанский, Майр, Тимофеев-Ресовский, Симпсон, Дж. Гексли и другие представители современной синтетической теории эволюции полагают, что «все эволюционные изменения в прошлом вызывались причинами, действующими в настоящее время и доступными поэтому экспериментальной проверке» [Dob-zhansky, 1937, с. 4] и потому относятся по классификации Завадского и Колчинского [1977] к «неоуниформистам», появляется много сторонников «автогенетического», «экогенетического», «циклического» и «ациклического» «неокатастрофизма» и их разновидностей, а также «вариаформизма», согласно которому утверждается «возможность различных форм изменений как общих, так и частных факторов и причин эволюции» [Завадский, Колчинский, 1977, с. 49].
Не имея возможности рассматривать и обсуждать все эти направления, отметим, что особенно много «неокатастрофистов» встречается среди геологов и палеонтологов. Так, согласно взглядам палеонтолога Соболева (1924), история органического мира слагается из малых и больших биогенетических циклов, включающих периоды эволюционного развития, прерываемые революционными эпохами. Он находит, что коренные преобразования фауны происходили на границах протерозоя и палеозоя, девона и карбона, перми и триаса, мезозоя и кайнозоя, совпадали по времени с началом циклов орогеыетических процессов и полагает, что они были детердшниррваны космическими циклами.
^"иГиндевольф [Schindewolf, 1954, 1963] также обращает внимание на резкие изменения в составе фауны на границах эр и периодов геологической истории и выделяет циклы из нескольких фаз, в которых скорость вымирания и появления новых форм различаются. По мнению Шиндевольфа, роль отбора — в элиминации типов, существование которых в данных условиях невозможно. Ускорение темпов эволюции в определенные фазы Шиндевольф "предположительно связывал с изменением космической и солнечной радиации. Аналогичные представления развивали Красовский и Шкловский [1957].
Личков [1945] отмечал, что «каждый раз революция в органическом мире всеземного масштаба происходила в эпоху образования складок земной коры или орогенеза» (с. 169).
Сушкин [1922], Яковлев [1922], Красилов [1973] и некоторые другие ученые развивали представление о решающем влиянии изменений климата на эволюцию жизни. При том смену таксономических групп первые два автора объясняли не вытеснением од
16

ной группы другой, а занятием места, освобожденного вследствие вымирания от неблагоприятных условий.
Ряд авторов связывал критические эпохи в развитии жизни с циклическим усилением вулканизма [например, А. П. Павлов, 1924], изменением интенсивности солнечного освещения [Голен-кин, 1959], орогенезом [Криштофович, 1957], усилением миграции радиоактивных элементов к поверхности Земли [Иванова, 1955].
Этот перечень воззрений можно значительно увеличить, однако и сказанного достаточно, чтобы заключить, что jBongoc о соотношении эволюционных и революционных фаз в гзазвитйГ^ Еашщ; и оргаА Мы закон-
чим его выдержками из недавно переведенной книги Гранта «Эволюция организмов» [1980]: «Давно известно, что геологическая история жизни отмечена эпизодами массового вымирания, разделенными длительными периодами постепенного эволюционного изменения. Палеонтологическая летопись заставляет предполагать нечто среднее между катастрофизмом Кювье и строгим униформизмом Лайеля — состояние эпизодической эволюции, В истории периодически возникают эпизоды массового вымирания одних групп, за которыми следует развитие и адаптивная радиация новых групп [Newell, 1967]... Существенная черта эпизодов вымирания — высокая скорость вымирания, наблюдаемая в совершенно не связанных между собой группах примерно в одно и то же время. Так, у аммоноидов и рептилий отмечается быстрое вымирание в поздней перми, позднем триасе и позднем мелу... Большинство этих эпизодов вымирания можно коррелировать с периодами поднятия материков и установления континентального климата. Такая корреляция позволяет считать, что возмущения, происходившие в среде, запускали сложные цепные реакции, приводившие в конечном счете к вымираниям в крупных масштабах [Newell, 1967; Axelrod, 1974]. ...За эпизодами массового вымирания следует образование и развитие новых групп... Необходимое условие для возникновения новой группы — наличие соответствующих экологических возможностей, которые и создаются в результате массового вымирания» 3 (с. 328—329).
Близкие взгляды развивает старейший советский геолог Д. В. Наливкин. Он отмечает, что «Земля является частью Вселенной, и поэтому крупные события, происходившие в Солнечной системе и в Галактике, влияли на ее развитие и строение. В истории Земли существовали периоды спокойного эволюционного развития и активных революционных перестроек... Эволюционные периоды по длительности превышают революционные... В революционные периоды континенты обычно поднимались, море отступало, т. е. происходили максимальные регрессии» [1980, с 6-7].
8 Последнее заключение, как видим, совпадает с ранее высказанными представлениями Сушкина [1922] и Яковлева [1922].
17

Борьба идей в отношении развития Земли й органического мира оказала доминирующее влияние на взгляды естествоиспытателей. Воспринятая Ламарком, Лайелем и другими учеными точка зрения Лейбница, что «природа не делает скачков»/была принята рядом ведущих математиков, в результате чего возникли соображения типа приведенных Лоренцом [1897] в защиту необходимости использования только непрерывных моделей: «...если даже принять возможность разрывного изменения в природе, то все же можно при исследовании подобного явления заменить его процессом изменения непрерывного и потом определить, как отразится на полученных результатах увеличение скорости изменения». Именно такие соображения привели к созданию основного современного математического аппарата, например уравнений математической физики. Такие подходы к моделированию, описывая процесс развития на эволюционном участке, не позволяют рассматривать траектории развития, включающие переходы через точки качественных изменений.
Решающее влияние на осознание принципиальной роли скачков в развитии биологических систем связано с введением де Фризом понятия о мутациях, названных Шредингером [1972] «квантовыми скачками» в генной молекуле, которые являются скачкообразными изменениями, происходящими без промежуточных состояний.
Позже идея квантованности физических процессов привела к созданию квантовой физики. С. И. Вавилов [1943] писал: «Новая физика в некоторых пунктах отказалась от идеи непрерывности, идея атомизации, скачков, прерывностей глубоко проникла в современную науку. Атомируется масса, электрический заряд, энергия, действие...» (с. 173).
Шредингер [1972] считал, что время ожидания мутации с порогом энергии W определяется зависимостью t=% exp W/kT, где т — период колебаний, к — постоянная Больцмана, Т — температура.
В связи с тем, что pK=const при 77=const, где X — длина свободного пробега, р — давление, a p=RT/V (V — объем, R — газовая постоянная) и Я=£/У, где U — скорость частицы, т* — время движения между двумя соударениями, получим
RT/V-Ut*=C или
%«=CVlRXJT-"=CT-\
При фиксированном пороге энергии t==r exp ат* время между двумя последовательными мутациями будет экспоненциальной функцией времени цикла функционирования частицы. Отсюда т*=1/а In tlx.
По экспериментальным данным, приведенным Поповой [1971] для щуки, ^(^Зе1,4^. Величина т* есть время запаздывания, в свя
18

зи с чем последовательность мутаций связана с характеристиками «памяти» системы об ее предыстории. Здесь речь идет о принципиальном влиянии «памяти» о предыдущих во времени состоя^ яиях на текущие характеристики системы.
Аналогичным образом многократно отмечается влияние пространственной «памяти» — действия соседних элементов на развитие данного [Гурвич, 1977; Дьюкар, 1978; Светлов, 1978; и др.] •
Таким образом, уже из термодинамической модели мутаций Шредингера [1972] следует необходимость учета характеристик памяти систем и их связей с соседними элементами для описания процессов типа мутаций — скачков в развитии. В то же время основные модели, используемые в настоящее время для описания развития биологических систем, как правило, включают только зависимость скорости роста системы от ее размера и возраста н данный момент и не учитывают влияния предыстории ГМина, Клевезаль, 1976].
Практика моделирования процессов развития показывает, что модели определенного вида оказываются справедливыми во вполне определенном диапазоне изменения параметров объекта и времени его развития. Связано это, как правило, с переменным характером внутренних и внешних условий развития объекта.
Установлено, например, что при постоянных внутренних и внешних условиях развития расширенное воспроизводство происходит по закону слржных процентов или соответствует экспонен-. циальному закону с постоянным значением относительных приростов. С изменением этих условий, начиная с некоторого момента времени, потребуются адаптационные изменения в структуре системы и характере ее регуляций, направленные на приспособление к новым условиям развития [Эшби, 1962].
Таким образом, здесь находит отражение диалектический закон перехода количественных изменений в качественные, который должен найти выражение в структуре используемых для моделирования экономических процессов описаний. В связи с этим необходима разработка моделей:
— описывающих изменение характеристик объекта при определенном принципе его действия (эволюционные модели);
— предназначенных для определения сроков использования определенного принципа действия, и в том числе моментов перехода к действию нового принципа.
19

1.3. Критические уровни в развитии сложных систем
Моделирование критических уровней развития биологических систем, под которыми понимаются моменты коренных качественных изменений, представляет существенный теоретический и практический интерес. Критические уровни определяют диапазоны, в которых система имеет различную морфофункциональную структуру, принципы регулирования, а переход через этот уровень связан с соответствующими структурными и регуляторными изменениями.
Между последовательными критическими уровнями система сохраняет свои качественные свойства, в связи с чем здесь для нее должна быть характерна низкая чувствительность к внешним и внутренним изменениям условий развития (высокая устойчивость). В таких диапазонах можно ожидать пониженной реакции на внешние воздействия, в связи с чем и модели для описания таких стадий развития могут быть предельно простыми (с малым числом параметров) без потери точности описания.
Диапазоны повышенной чувствительности к внешним воздействиям зафиксированы в экспериментальной эмбриологии [Светлив, 1960] и получили название критических. В эти периоды отмечается повышенная отбраковка дефектных эмбрионов при действии самых различных раздражителей. Такова, например, стадия гастг^ляции, отбраковка на которой сопоставима с отбраковкой прнГрождении.
В промежутке между гаструляцией и рождением имеется еще ряд критических периодов, однако их значимость существенно ниже. Отсюда возникают проблемы классификации типов критических уровней, установления рангов их значимости и определения характера изменений, свойственных каждому рангу.
Наличие критических уровней на популяционном уровне организации биологических систем наблюдается в виде стабильных амплитуд колебаний численности популяций («волны жизни» Четверикова), при этом верхняя и нижняя границы поддерживаются на приблизительно постоянном уровне в течение длительных промежутков времени (многих десятков и даже сотен лет). Выявление теоретических пределов развития биологических систем на различных иерархдческих уровнях, установление взаимосвязей между уровнями и общих количественных закономерностей их формирования представляют большое значение для теории развития биологических систем.
Практические результаты, связанные с определением критических уровней, необходимы для управления развитием биологических систем. Промысловая эксплуатация популяции до уровня ниже определенной критической численности приводит к ее распаду. Определение допустимых пределов эксплуатации — задача исследования. Одно и то же воздействие на организм может оказаться повреждающим в одной фазе его развития, лечебным в другой фазе, либо вообще не оказать заметного влияния — как клас-
20

инфицировать фазы развития по их чувствительности к управляющим воздействиям? Подобного рода вопросы, специализированные для систем каждого иерархического уровня, можно формулировать до достаточно детальных и конкретных постановок задач. Однако общим для них является необходимость определить критические уровни в развитии биологической системы на разных уровнях иерархии и детализировать особенности ее функционирования в рамках каждого критического диапазона.
Одной из задач данной книги является исследование вопроса об общих подходах к моделированию критических уровней развития биологических систем, о взаимосвязи характеристик критических уровней и формировании критических диапазонов развития для последовательных уровней иерархии.
Другая задача связана с проверкой действенности разработанных нами моделей на примерах оценки размерных диапазонов уровней биологической организации, критических уровней популя-ционной иерархии, онтогенетического и филогенетического развития, а также на примере развития земной коры и жизни.
В заключительном разделе книги рассматривается вопрос о возможности распространения найденных закономерностей на другие системы и формулируется закон критических уровней ал-лометрического развития систем.

2. ФОРМИРОВАНИЕ КРИТИЧЕСКИХ УРОВНЕЙ
В настоящем разделе излагаются основные механизмы формирования критических уровней (2.1), которые в дальнейшем положены в основу аксиоматики разрабатываемых моделей. Проведенная здесь классификация основных типов развития (2.2) связана с выделенными механизмами и предварительно характеризует принципиальную возможность моделирования критических уровней раз-
вития.
ГШ*.
В заключительной части раздела (2.3) рассматриваются влияния внешних и внутренних факторов на положение критических рубежей.
2.1. Формирование критических уровней развития биологических систем
Оценка функциональных возможностей биологических систем связана с необходимостью детального исследования механизмов их адаптации к изменяющимся условиям среды. Шмальгаузен [19463 писал, что «не существует никакой организации, даже самой глу боной, которая не была бы приспособительной уже в своей основе. Вне приспособления никакой организации не существует...» (с. 146).
Взаимодействие среды и биологической системы может определяться случайным характером ритмики параметров внешней среды, запаздывающим (релаксационным) характером реакции^оога-низма^на внешнее воздействие [Серов, 1974] и огг^ережающими ^хантамами отражения действительности в процессе адаптации I^hoxhhTI962; Наумов, 1963; Георгиевский и др., 1975].
Релаксационный характер биологических процессов, а также опережающий характер отражения действительности при адаптации моделируются дифференциально-разностными уравнениями [Вольтерра, 1976; Кузьмин, Ленская, 1974]. В пространстве параметров, характеризующих развитие системы, описываемой дифференциально-разностными уравнениями, существуют области, где решения являются: экспоненциальными с затухающими колебаниями; экспоненциальными с колебаниями, имеющими экспоненциально возрастающие амплитуды; колебательными.
колебания являются существенным свойством релаксационных систем, однако при этом важным фактором оказывается характер этих колебаний — затухающий или расходящийся — и характер изменения общей тенденции процесса, соответствующий росту или убыванию исследуемого фактора. Цянь Сюэ-сэнь [1956] и Эшби [1962] формулируют задачу обеспечения адаптации систе-
22

мы к новым условиям как определение области ее устойчивого функционирования и изменения ее параметров таким образом, чтобы цри выходе на границу этой области обеспечить возвращение системы внутрь области.
Колебания с экспоненциально возрастающими амплитудами соответствуют неустойчивому равновесию объектов и практически в нормальных процессах развития биологических объектов не реализуются, в связи с чем области в пространстве па--раметров, где могут иметь место такие случаи, можно считать запретными.
Переход объекта при выходе на границу области устойчивости внутрь этой области считают происходящим скачкообразно [Насонов, 1959; Эшби, 1962]. Именно таким моментам соответствуют механизмы адаптации к новым условиям развития объекта. Изменение параметров развивающейся системы происходит так, чтобы было обеспечено ее устойчивое существование в новых условиях, т. е. успешно решена проблема адаптации. Обеспечение устойчивого развития объекта в будущем может быть произведено только путем использования механизмов опережающего отражения действительности [Анохин, 1962]. Как показал Н. Н. Красовский [1968], решение задачи об определении текущего состояния объекта с запаздыванием по результатам измерения его характеристик в прошлом требует решения уравнений опережающего типа, т. е. связано с необходимостью опережающего отражения действительности.
Для механизмов управления, соответствующих опережающему отражению действительности, будет иметь место своя область устойчивости, границы которой в общем случае могут не совпадать с границами областей устойчивости развития биологического объекта в соответствии с запаздывающим характером его реакций.
Если границы критических уровней для опережающих механизмов управления шире, чем для процесса развития запаздывающего типа, то процесс может быть
Рис. 3. Схема формирования последовательности адаптационных перестроек По абсциссе — времена обменных и регуляторных процессов (запаздывание) — т; по ординате — темп экспоненциального роста — u. Заштрихованы запрещенные для развития зоны. а — область роста, б — область уменьшения размеров [Из: Кузьмин, Жирмунский, 1980]
Et
:/ -+
1 I  I  \ I Г I 1 . 5 7 9H1315171S Сутки
Рис. 4. Последовательные стадии экспоненциального роста эмбриона курицы (1-Ш)
По абсциссе — возраст; по ординате — масса (шкала логарифмическая) [По: Шмальгау-зен, 1935]
23

стабилизирован любым изменением параметров, обеспечивающим переход внутрь области устойчивости. Если же эти границы находятся внутри области для процесса развития запаздывающего типа, то для стабилизации процесса необходимо изменить его параметры так, чтобы обеспечить попадание внутрь области устойчивости для обоих процессов (рис. 3). В ряде практически реальных случаев имеет место именно такой вариант [Кузьмин, Ленская, 1974].
Таким образом, механизмы запаздывающего типа, характеризующие релаксационные механизмы реакций биологических объектов на внешние воздействия и внутренние процессы развития и опережающего отражения действительности, формируют области
устойчивого развития.
Шмальгаузен [1935] показал, что процесс роста животных начинается экспоненциальной фазой. Последовательность таких экспоненциальных фаз — это последовательность переходов между (^едаШи^моментами адаптационных перестроек. В кривой роста (риег 4) такое развитие представляет собой последовательность экспоненциальных фаз. При этом, если следовать аллометрическо-му закону роста, который соответствует линейному изменению времени запаздывания с возрастом, получаем ритмический процесс с модулем, соответствующим отношению времен двух последовательных адаптации и равным величине скачка в темпах развития между границами критических областей для запаздывающего характера развития и опережающего механизма адаптации. Такая последовательность может рассматриваться как ритмическая. Для процессов устойчивого роста была получена оценка величины скачка в темпах, равная 1,43 [Кузьмин, Ленская, 1974]..
Из представленных соображений следует, что изменение времен запаздывания реакций биологических объектов в определенном диапазоне значений может не приводить к изменению протекающих в них процессов. Адаптационные перестройки будут происходить только при достижении запаздываниями некоторых
критических уровней.
Интерес к критическим точкам в развитии биологических объектов нашел отражение в исследованиях Пэрна [1925], выделявшего «узловые точки», характерной для которых он считал особую чувствительность организма в этот период, повышенную ранимость, зависящую от «частных перестроек». К аналогичным выводам пришли Свобода [Swoboda, 1904], Флис [Fliess, 1925] у Томмен [Thommen, 1976], которые разработали представления об околомесячных биоритмах человека.
Величины запаздываний характеризуют времена обменных и регуляторных процессов, происходящих в биологических системах,, и зависят от эндогенных факторов: размер животного [Винберг и др., 1968], его видовая принадлежность, возраст [Стреллерг 1964], и экзогенных условий (освещенность, температура, соленость воды для морских организмов, количество пищи и т. п.). Таким образом, эндогенные условия могут формировать ста-
24

200 fi
Рис. 5. Развитие с промежуточной дифференцировкой
Изменение массы желудка мальчика с возрастом [По: Клейн, Кришун, 1967]
10
■ ■ I I I | ' ' | t ' ' ' ' » ' » ' ч /   J  5   7  9  11 13 15 17 19 Годы
Рис. 6. Задержка и последующее ускорение роста опухоли
I — первая положительная фаза опухоли; II — отрицательная фаза — объем опухоли
уменьшается; III — вторая положительная фаза [По: Эмануэль, Евсеенко, 1970]
Рис. 7. Циклический процесс. Изменение массы семенников (с придатками)
у черного хоря Mustela putorius в течение года
[По: Данилова, Русакова, 1969]
ционарную составляющую динамики изменения времен обменных процессов, тогда как экзогенные4 условия задают колебания этой величины.
В развитии биологических объектов наблюдаются случаи, когда уменьшение темпов роста системы (рис. 3, квадрант А) оказывается недостаточным для ее стабилизации, и необходима коренная перестройка ее морфофункциональной структуры. Такая перестройка происходит с уменьшением размеров (или веса), т. е. объект использует для ее реализации имеющиеся ресурсы. Перестройка позволяет достичь более высоких темпов роста, чем те, что были зафиксированы перед перестройкой. Такие случаи достаточно характерны для биологических объектов [например, Эмануэль, Евсеенко, 1970].
При исследовании процессов развития клеток Вермель и Пор-тугалов [1935] отметили некоторое уменьшение их размеров перед периодом интенсивного роста. Они также заметили, что существует пропорциональность между модами (максимумами) кривой изменчивости ядер клеток разных возрастов, стадий диф-ференцировок и видов тканей. Эта пропорциональность весьма близка к геометрической прогрессии со знаменателем 2, т. е. 1:2:4... Однако Вермель [1935] указывает на существование Рядов, где знаменатель прогрессии близок к 1,5. Стабильный ха-
25

рактер получаемых соотношений привлекает внимание возможной общностью механизмов, лежащих в основе таких процессов.
ч Исследование процессов роста биологических систем, имеющих режимы спада параметров в определенные моменты времениг представляет интерес для решения ряда задач управления развитием. Так, приведенная на рис. 5 кривая изменения веса желудка с возрастом у мальчиков имеет точки максимума, приходящиеся на 4 и 11 меряцев, 6, 12 и 19 лет, что достаточно хорошо соответствует переходным возрастам [Автандилов, 1973]. В эти периоды,, которые соответствуют перестройкам в организме, воздействия различного рода будут оказываться наиболее повреждающими,, так как здесь система дезорганизована для осуществления предстоящей перестройки, что фиксируется на кривой роста в виде участков спада.
Такой механизм обсуждается в ряде работ. Шапошникова [1973], например, считает, что «сопротивляемость организма к различным повреждающим воздействиям внешней среды изменяется и значительно уменьшается в отрицательной фазе» (с 57).
Шмальгаузен [1935] отмечал, что «периоды депрессии, во-первых, явно совпадают с периодами максимального накопления некоторых продуктов обмена в эмбриональном организме..., а во-вторых, сопровождаются более интенсивными формообразовательными процессами. При развитии некоторых органов (легкие,, конечности) удается особенно явно обнаружить правильное чередование периодов усиленного роста и преобладающего формообразования» (с 41).
Один из возможных механизмов, приводящих к необходимости дифференцировки, состоит в том, что возможности управления системой в рамках данной морфологической структуры ограничены.
Очевидно, в процессе развития наступает такой момент, когда возможности управления биологической системой исчерпываются, в связи с чем становится необходимой дифференцировка как способ сохранения жизнеспособности при последующем развитии системы.
В ряде случаев обращается внимание на наличие на кривой развития участков спада характеристик, но природа этого явления детально не анализируется. Так, Эмануэль и Евсеенко [1970] г рассматривая кинетическую кривую роста медленно развивающегося плоскоклеточного рака правой нижней доли легкого (рис. 6), констатируют, что, достигнув значительной величины, опухоль самопроизвольно уменьшалась в размерах, а затем после некоторого периода задержки рост возобновлялся с несколько большей скоростью.
Принципиальным в данном случае является неповторимый характер перестроек, происходящих в каждом периоде спада характеристик объекта для случая, когда размер системы в конце фазы развития с отрицательным темпом превосходит ее размер в начале развития в положительной фазе. 26 .

Когда в конце отрицательной фазы биологический объект достигает того же размера, какой он имел в начале положительной фазы (рис. 7), развитие становится циклическим, т. е. последовательные фазы повторяют друг друга. Здесь, в отличие от ранее рассмотренных случаев, имеет место повторяющаяся последовательность адаптационных перестроек.
Следовательно, в процессах развития наблюдаются адаптационные перестройки, ко^рыеНпрсш^^ повторяющиеся промежутки времени. Природа их принципиально Jjgg-v лична1 в связи с чем возможности^шГ использования для целей регуляции состояния также. принципиально различаются.
^ представляют ритмы функциональной актив-
ности, которые формируются за счет взаимодействия разных систем организма. Исследования в этом направлении достаточно фрагментарны и относятся к попыткам установления взаимосвязей между отдельными органами и системами организма. Нами [Жирмунский, Кузьмин, 1979] была предпринята попытка использовать в этих целях древнюю схему взаимодействия между функциональными системами организма, которая позволяет упорядочить разрозненные классические [Боткин, 1950; Захарьин, 1910] и современные [Корнева, 1965; Климов, 1976; Маркова, 1976] данные об этом явлении.
Серьезной проблемой в исследовании критических уровней развития является согласование их параметров на разных ступенях иерархии. К сожалению, в настоящее время публикаций об исследованиях биологических объектов одновременно на ряде основных уровней иерархии очень мало [Терскова, 1975; Гнездило-ва и др., 1976]. Тем не менее на имеющемся достаточно ограниченном материале удается наблюдать удивительное согласование критических характеристик процессов на различных уровнях.
Синхронизация динамических систем, состоящих из большого числа элементов со слабыми связями между ними, была математически обоснована Андроновым и др. [1967] и Блехманом [1971]. Количественные исследования параметров синхронизации процессов на разных уровнях представляют проблему для детального теоретического и экспериментального изучения на различных биологических объектах. Эмме [1967] обсуждает причины патологических состояний в длительных циклах, когда в норме клетки центра или ткани функционируют несинхронизованно, а шок или травма могут синхронизировать фазы всех клеток и придать им новую длительность циклов.
Дискретный характер адаптации^ в к^итетесщ^^чках нашел отражение в концепции кризисных явлевдн^ ГСерелом, кризис составлял важную часть* учения Тйппократа [Гиппократ, 1848]: каждая быстрая перемена, производящая в болезни излечение или смерть, значительное улучшение или усиление, называлась переломом. Считалось, что эти переломы имеют место преимущественно в определенные дни после начала заболевания. Принципиальную неоднородность процессов в организме отмечал Пэрна [1925],
27

который считал, что «волну характеризует не возрастание и убывание, а смена фаз» (с. 17).
Отмечают большую роль эндокринных факторов в происхождении периодических психозов, так как начало болезни в этих случаях обычно совпадает с каким-нибудь периодом гормональной перестройки: половое созревание, первая неделя после родов,, выкидыш, климактерический период [Сухарева, 1974].
Сам характер ритмики, имеющей место в организме, обычна существенно отличается по природе от гармонических колебаний, которые используются при построении моделей функционирования органов и систем в организме [Шумаков и др., 1971]. В частности, в циркадном ритме, как и в длиннопериодических ритмах,, активная фаза отделена от пассивной качественным переходом, скачком, переключением механизмов управления на принципиально другой режим работы, в связи с чем описание таких ритмов, как гармонических колебании, за внешним сходством характеристик скрывает принципиальное различие самих механизмов функционирования организма на различных стадиях развития.
Таким образом, в основу построения модели, позволяющей определять критические уровни развития систем, необходимо положить учет влияния релаксационных характеристик на развитие системы, а также их зависимость от внутренних и внешних факторов [Кузьмин, Жирмунский, 1980в].
Задачей исследования является определение степени влияния исходной информации на положение критических уровней развития.
2.2. Классификация основных типов развития
Релаксационные свойства оказывают существенное влияние на характеристики развивающейся системы и формирование ее критических уровней. Они выражаются в явлениях «памяти» о предыстории развития, если речь идет об анализе динамики системы, и «нелокальности», которая характеризует влияние соседних элементов на состояние данного (при этом анализируются процессы формообразования или развития системы в пространстве). Таким образом, при описании развития биологических систем необходимо использовать такие модели, которые учитывают влияние пространственных («нелокальность») или временных («память») запаздываний.
Показано   [Кузьмин, Ленская,  1974], что в качестве такой модели можно использовать уравнение вида
x=kx(t—т), (1)
где х — размер системы, т — запаздывание, к — константа скорости роста. Здесь скорость роста размера системы х зависит от ее размера в момент времени (t—т), при этом величина т характеризует время (или пространство), смещение состояния системы на которое определяющим образом влияет на развитие в данный
28

момент (либо в данной точке). Так, в случае, когда х — численность популяции, скорость ее роста пропорциональна не общей численности x(t), а количеству половозрелых особей, т. е. родившихся т единиц времени тому назад: x(t—т), где т — срок созревания.
В разделе 3 будет показано, что уравнение (1) обладает достаточно общими свойствами, которые позволяют рассматривать его как обобщение целого ряда моделей, используемых для описания роста.
Величина запаздывания является переменной в ходе процесса развития и определяется совокупностью внешних и внутренних условий. В большом количестве случаев переменный характер запаздывания непосредственно определяется изменением размеров системы. В частности, рост системы при фиксированной ее организации приводит к увеличению времени как на управление ею, так и на перенос внутри системы потоков вещества (подвод пищи, отвод продуктов жизнедеятельности и т. д.). При этом времена на прохождение управляющих команд ху растут с увеличением системы [Флейшман, 1971]. В линейном приближении будем считать, что ху=а±х.
Для переноса потока вещества потребуется время тем меньшее, чем больше потребляющих единиц, в связи с чем в первом приближении время, необходимое для доставки и усвоения потока вещества тм, можно принять обратно пропорциональным размеру системы, т. е. хм=а2/х. Полное время цикла системы будет представлено суммой времен x=xyJrxM—diX+a2/x. Параметры аА а2 являются функциями способа организации системы.
Когда в рамках одной организации размер системы х небольшой, преобладающей является вторая составляющая времени цикла тм, тогда как в случае больших размеров системы будет преобладать величина ху. В связи с этим рассмотрим два крайних случая.
Первый из них характеризуется ранними стадиями развития системы, т. е. здесь Тм>ту.
Уравнение процесса роста для него запишем в виде x=k2x(t—a2/x), (2)
тогда как во втором случае сложность системы такова, что ху>хМг тогда
x=kix(t—aix). (3)
Для уравнения (2) в соответствии с теоремой сравнения [Эльсгольц, Норкин, 1971] для к2<0 при
к2а2/х0<1/е
существует и конечен lim ^ф (t) e^MfaWxo)*
t -*оо
пРи Х+е~Ща2/Хо=0. Когда к2а2/х0^1/е, получим X+e~Ve=0, т. е. А=

=-ё, x*(*)^Ce~vV При k2ajz0*0 Я=-1 и x,(t) -*Се~*«'. Таким , образом, показатель экспоненты в области отрицательных темпов | «будет меняться не более чем в е раз. L
Произведем оценку значения я0, при котором предел имеет j место Хо>к2а2е< 1
Значит, существует минимальный размер системы, обеспечивающий   экспоненциальный   характер   решения   уравнения   (2)
при ц<0. ,
Использование результатов теоремы сравнения для уравнения (3) позволяет сказать, что для kt<0 при kiatx0<l/e существует и конечен
lim Xq> (t) exp к\К (kiaix0) tf
при Я+ехр (— Хк^Хо)^. Когда й1а1х0=1/е, получим результат, аналогичный предыдущему, т. е. А,+ехр (—К/е) =0 и Я=—е.
При kidtXo^O Я=—1, что также соответствует результату, полученному для предыдущего случая.
Произведем оценку значения x0l при которой предел имеет a:o<l/ea1Ai. Таким образом, существует максимальный размер системы, обеспечивающий экспоненциальный характер решения уравнения (3) при »<0. Значит, при больших размерах системы, в соответствии с уравнением (3), имеется максимальное значение, обеспечивающее экспоненциальный характер развития процесса в области отрицательных темпов, а при меньших размерах системы, в соответствии с уравнением (2), имеется минимальное значение размера системы, обеспечивающее экспоненциальный характер развития в области отрицательных темпов.
Здесь возможны следующие варианты взаимного расположения областей, из которых достижима область экспоненциальных режимов развития при отрицательных темпах:
1. ea2ft2<l/ea1ifc1 или еУА4А:2а1а2<1,
2. ea2^2>l/eö1Ä;1 или eVA1A2a1a2>l
3. ea2k2=l/ealki или ^к^к2а^2=\.
Первый из рассматриваемых вариантов соответствует наличию! области (рис. 6), где может происходить рост системы в области! положительных темпов (до границы 1/eaifci: моментti на рис. 8),| а затем будет происходить падение размера до границы еа2к2 (мо-| мент t2). Дальнейшее падение величины х в рамках модели экспо-1 пенциального роста становится невозможным. 5
Таким образом, при неизменных значениях параметров а{ и a2i\ предполагая, что могут реализовываться только экспоненциальные; режимы развития, получим класс циклических систем (рис. 9).;
Однако в связи с тем, что параметры ai и а2 являются функциями   организационной структуры системы, в отрицательной
30

А
1
вам
еагнг
еагпг
t
Рис. 8. Траектория развития системы внутри области 1
По абсциссе — возраст системы (t), по ординате — ее размер (х). Границы выделяют область развития системы. Объяснения см. в тексте
Рис. 9. Траектория циклического развития системы Обозначения см. на рис. 8
фазе развития, т. е. при уменьшении размеров системы, может произойти изменение организационной структуры, что приведет к изменению параметров at и а2, а значит, и к сдвигу границ роста размеров в очередной фазе развития.
Таким образом, предельные размеры системы после перестройки должны определяться величиной скачка в параметрах а4 и а%. В целом кривая роста с учетом возможной смены параметров at и а2 за счет перестройки в отрицательной фазе может иметь вид, представленный на рис. 10. ;
В соответствии с полученной моделью размер системы в процессе развития лимитируется сверху механизмами управления, а снизу — характеристиками материальных потоков в системе Отсюда возникает принципиальная возможность перехода в определенной фазе развития к циклическому процессу. *
Полный набор возможных переходов представлен на рис Щ
I — начальное развитие системы — до перестройки; II — развитие Системы со сменой параметров ах и а2 и переходом из области 1 в область 2\ III — циклическое развитие на новом уровне. Остальные обозначения см. на рис. 8
Рис. 10. Траектория развития системы с учетом смены параметров ai и а2
t
31

2.3. Влияние внешних и внутренних факторов на формирование критических уровней биологических систем
Исследованиям влияния внешних (экзогенных) и внутренних (эндогенных) факторов на развитие биологических систем посвящено большое количество работ [например, Шмальгаузен, 1946; Прос-сер, Браун, 1967; Тимофеев-Ресовский и др., 1977]. Так как задачей настоящей книги является анализ количественных закономерностей развития биологических систем, мы рассмотрим конкретные результаты измерений, связанных с влиянием внешних и внутренних факторов на формирование критических уровней развития биологических систем.
Влияние внешних факторов на развитие и функционирование биологических систем хорошо изучено на примере температуры. Известны видовые адаптации к температурным условиям обитания. Животные, населяющие более теплые географические или вертикальные зоны, переносят действие более высоких температур, нежели животные, обитающие в более холодных зонах [Kinne, 1954; Fry, 1957; Ушаков, 1960; Жирмунский, Шляхтер, 1963; Проссер, Браун, 1967; и др.].
Винберг [Винберг и др., 1968] приводит зависимость продолжительности жизни различных видов холоднокровных животных (D) от температуры (Т)
В = Ле-°'083Г. (4)
Стреллер [1964] для ряда беспозвоночных указывает зависимость скорости старения от температуры, имеющую коэффициент температурной чувствительности того же порядка, что и в соотношении (4). Аналогичная зависимость с таким же коэффициентом температурной чувствительности получается и для периодов дозревания рыб. Нами обработаны материалы Кошелева [1968] по зависимости времени наступления половой зрелости волжских лещей (т) от температуры. Результаты обработки, приведенные на рис. 11, показывают, что зависимость имеет вид
т = A1e-Q>o™T. (5)"
Таким образом, влияние температуры на продолжительность* жизни, характеристики старения и сроки созревания подчиняются зависимости одного вида.
Условия существования оказывают существенное влияние так-; же на темп роста животных. Селин [1980] исследовал рост миди" Грея (Crenomytilus grayanus) в обрастаниях на якорной цепи и на естественном субстрате на разных глубинах. Как видно на рис. 12, скорость роста оказалась выше для мидий из верхних торизонтов (6 и 12 м), ниже с 24 м из обрастаний цепи и наибо лее низкой с 20 м из естественных поселений на валунах. По-видимому, на скорость роста оказывают влияние как температур
i
32

Рис 11 Зависимость возраста наступления половой зрелости лещей от температуры окружающей среды По абсциссе — температура воды в июле; по ординате — возраст наступления половой зрелости (шкала логарифмическая) [По: Кошелев, 1968]
Годы} 7 5\
15
20
25   30°С
Рис. 12. Кривые линейного роста мидий Грея (Crenomytilus gray anus) из обрастаний якорной цепи на глубинах 6 м (1), 12 м (2), 24 м (3) и из естественного поселения на валунах (4)
По абсциссе — возраст моллюсков; по ординате — длина раковины. Вертикальные линии — 99,9%-ные доверительные интервалы [По: Селин, 1980]
Рис. 13. Зависимость температуры, приводящей к остановке мерцательного движения клеток жабр моллюсков из разных географических зон за 1 мин, от средней для географической зоны температуры воды на поверхности в августе По абсциссе — средняя для зоны температура воды; по ординате — температура, приводящая к остановке мерцательного движения за 1 мин. [По: Жирмунский, 1971]
50 40
30
ri I < ' > ' '—I— 4    S    12г16 Годы
14
Рис. 14. Численность населения Земного шара в зонах с разными температурными условиями
По абёциссе — средняя годовая температура; по ординате — численность населения [По: Бар-тон, Эдхолм, 1957]
18
ШН ЧЖШ
SOOh
200
100
24
28°С
10
20    Ш
2   А. В. Жирмунский, В. И. Кузьмин

воды, так и другие условия, на что указывает существенно более низкий темп роста мидий на валунах.
Аналогичные данные приводит Бейн [Bayne et al., 1976] для двух групп личинок Mytilus edulis, выращенных при разных температурах, и для пяти популяций этих моллюсков из разных мест обитания, изученных разными исследователями.
В большом числе исследований доказана связь между теплоустойчивостью клеток и температурой существования видов. На рис. 13 приводятся данные, характеризующие соответствие между средними температурами, приводящими за 1 мин к остановке мерцательного движения клеток жабр моллюсков из верхней сублиторали разных географических зон северо-западной части Тихого океана, и средними августовскими температурами воды на поверхности в южных частях ареалов этих моллюсков [Жирмунский, 1971]. Эти данные показывают связь видовых свойств клеток животных с условиями обитания видов.
Как и другие виды, человек живет в определенных температурных пределах. Бартон и Эдхолм [1957] приводят данные о численности населения Земли в зонах с разными температурными условиями. При температурах ниже 5° и выше 27° С численность населения резко сокращается (рис. 14).
Таким образом, температура как фактор среды оказывает существенное влияние на вырабатывающуюся в процессе эволюции наследственно закрепляемую приспособленность биологических систем к среде, а также на их способность изменять свойства при текущих сдвигах в среде.
Развитие биологических систем зависит также от внутренних факторов — организации* функционирования, интенсивности обменных и регуляторных процессов. (Хорошо изученной характеристикой, определяющей развитие организмов, является интенсивность обмена веществ. Как показал Винберг [Винберг и др., 1968], интенсивность обменных процессов может непосредственно использоваться как характеристика, определяющая рост организмов. Гипотеза о том, что величина относительного прироста пропорционально интенсивности обмен, приводит к хорошему согласованию получаемых кривых роста с реализуемыми размерами.
В литературе имеется большое количество исследований, посвященных описанию результатов экспериментов по измерению интенсивности обмена и влиянию на него различных факторов [см.: Винберг и др., 1968; Проссер, 1978]. Интенсивность обмена а при этом описывается аллометрической зависимостью от веса тела X
с] = аХ'\ (6)
Будем считать время обменных процессов т величиной, обратно пропорциональной интенсивности обмена т~1/д, откуда
34

с использованием зависимости (6) получим
Проссер [1977] приводит данные о убывании потребления кислорода с ростом температуры, что соответствует в принципе кривым Аррениуса и соотношению типа (5). Отсюда, имея в виду совместное действие температуры среды и размера биологической системы из (5) и (7), получим соотношение вида
В связи с тем, что экзогенные факторы, влияющие на развитие биологических систем, часто имеют циклический характер (колебания температуры, солнечной активности и т. д.), рассмотрим простейший случай совместного влияния эндогенного фактора (размера) и синусоидально меняющегося экзогенного фактора (например, температуры) на характер роста. Целью такого рассмотрения является определение основных особенностей кривых роста в этих условиях.
В разделе 2.2 было рассмотрено влияние размеров системы на формирование критических уровней развития. Рассмотрим процесс, в котором размер и возраст связаны аллометрической зависимостью (Шмальгаузен, 1935)
Тогда зависимость времени обменных процессов от размера и температуры (8) с учетом (9) и (10) перепишется в виде
Будем считать время обменных процессов определяющим запаздывание в уравнении процесса развития (1). Тогда динамика процесса будет описываться уравнением
T===aie~bTXb.
(8)
т = A exp (bo sin wt) tb.
(И)
процесс развития является устойчивым при выполнении условия к A exp (— bo sin wt) tb < 1/e,
откуда
(13)
Это неравенство представлено на рис 15, откуда видно, что при
2* 35

Рис. 15. Сезонное изменение температуры (а) и зависимость запаздывания от времени (б)
По абсциссе — время, t; по ординате: а — температура среды, Г; б — запаздывание, т; изменение запаздывания под действием эндогенных факторов (1), под действием эндогенных и экзогенных факторов (2). Рубежи А и Б соответствуют двум разным значениям правой части неравенства (13), при выходе на рубеж система попадает в критическое состояние. tu t2 — времена выхода в критическое состояние, когда действие экзогенных факторов, отсутствует. t\ t'2 — времена выхода в критическое состояние, когда на систему действуют одновременно эндогенные и экзогенные факторы
воздействии факторов внешней среды, которые уменьшают время запаздывания (повышают интенсивность обмена), длительность периода до перестройки растет по сравнению со средним уровнем. Повышение времени запаздывания (снижение интенсивности обмена) приводит к более раннему наступлению критического периода. Сезонная структура активности зимоспящих животных подтверждает эти результаты.
Слоним [1971] отмечает, что понижение температуры внешней среды — основной фактор, вызывающий спячку млекопитающих на фоне остальных необходимых условий. Если температура в норе становится ниже определенной для вида величины, наступает пробуждение животных. Такого рода вариант представлен в верхней части кривой рис. 15, где пониженная температура вызывает более раннее наступление критической точки — перехода системы в другое функциональное состояние.
Таким образом, если периодическим фактором является температура, то общий максимум длительности развития до критической точки будет иметь место, если животные находятся в области сезонного повышения температур.
В рассмотренном случае эндогенные условия определяют общую тенденцию развития процесса и диапазон, в котором будет происходить перестройка, тогда как экзогенные условия задают отклонения от основной тенденции. В зависимости от конкретных внешних условий определяется момент наступления критической точки ъ диапазоне, заданном эндогенными условиями развития. Очевидно, в принципе именно таким является соотношение влияний эндогенных и экзогенных факторов, влияющих на развитие биологических систем. Это обстоятельство определяет необходимость детальных исследований влияния основных эндогенных и экзогенных факторов на время запаздывания.
36

2 Критические уровни в процессах развития биологических систем


3. МОДЕЛЬ КРИТИЧЕСКИХ УРОВНЕЙ РАЗВИТИЯ СИСТЕМ
Исторически моделирование развития биологических систем шло в двух основных направлениях: первое было связано с использованием кусочно-непрерывных функций, описывающих процессы на отдельных участках без определения характера связи между отдельными участками; второе искало способы описания процесса в целом.
Сильной стороной кусочно-непрерывных моделей является возможность выделения диапазонов развития, где качественные свойства развивающейся системы не меняются. Переход к новым принципам и механизмам регулирования, качественное изменение мор-фофункциональных структур приводит к изменению характера роста, что и пытались отразить представители этого направления, Это была попытка выделить эволюционные фазы развития биологических систем, а скачки, качественные изменения свойств системы, предлагалось учитывать по изменению параметров траектории роста. При этом в рамках эволюционных участков развитие описывалось функциями с минимальным числом параметров, например показательными или степенными.
Основным недостатком этого направления являлось отсутствие возможности определения моментов, в которые происходит переход к новой зависимости. В связи с этим в течение достаточно длительного времени оно было вне поля основных интересов исследователей роста.
Успехи, достигнутые при описании всего процесса роста прерывными моделями, на некоторое время исключили из рассмот* рения недостатки этого подхода. В частности, основные трудности, которые предопределили потерю интереса к первому направлению, пытались преодолеть путем введения многопараметрических нелинейных моделей. Но в ряде случаев в них теряется смысл моделирования как способа оперировать на количественном уровне исходными измеримыми характеристиками системы, которые можно содержательно интерпретировать. Модели часто стали превращаться в многопараметрические аппроксимации с неопределенными областями и пределами применимости и количеством параметров, соизмеримым с количеством измерений.
Постановка вопроса о необходимости определения пределов использования модели определенного вида возникла достаточно давно. Так, еще в середине XVIII века Эйлер [1936] при разработ-т ке анализа бесконечно малых рассмотрел возможность использования показательных функций для моделирования роста численности популяций. Он рассмотрел модель роста численности населения «после потопа» и показал, что при постоянстве относи
37

тельных приростов в течение исторически непродолжительного промежутка времени получается такая расчетная численность населения, для прокормления которого не хватило бы ресурсов всей Земли. Таким образом, был поставлен вопрос о том, что модель роста в соответствии с законом сложных процентов при постоянном проценте имеет ограниченную область применения.
Эта же модель — рост численности организмов в геометрической прогрессии — была положена в основу теории Дарвина. Произведя соответствующие расчеты, Дарвин пришел к выводу, что численность даже медленно развивающихся организмов ограничивается гибелью потомства в результате борьбы за существование. По его расчетам, для слонов, с учетом того, что каждая слониха дает за свою жизнь 6 детенышей, через 750 лет потомство одной пары составило бы 10 млн. особей, тогда как на самом деле численность слонов невелика. Модель, использованная Дарвином, была геометрической прогрессией, которая меняла модуль при достижении некоторых критических уровней, определяемых давлением среды. Именно этот факт в последние годы был положен в основу глобального моделирования [Форрестер, 1978].
Возникают следующие вопросы, ответы на которые определяют направления развития методов моделирования.
— Каким образом минимизировать количество используемых при моделировании параметров без потери точности, свойственной экспериментальным данным?
— Как определить пределы использования моделей данного вида?
Первый из указанных вопросов связан с поиском унифицированных моделей роста, что нашло отражение в преимущественном использовании в научных исследованиях достаточно ограниченного набора функций — показательной, степенной, Гаусса. Необходимо также обсудить вопрос о том, насколько содержательны такие модели (не является ли их использование историческим пережитком) и каким мог бы быть набор унифицированных моделей.
Ответ на второй вопрос предопределяет качественную содержательность моделирования, так как гарантировать точность результатов моделирования можно лишь, зная диапазоны характеристик развивающейся системы, внутри которых она не меняет своих основных свойств.
Основное направление исследования, результаты которого приводятся в настоящем разделе, связано с анализом общих свойств моделей развития, определением обобщенного варианта модели, частным случаем которого являются многие модели, используемые для описания роста. Для этой модели путем исследования ее качественных свойств определяются простейшие модели, описывающие развитие.
Биологической основой целесообразности и необходимости тако-то исследования является принципиально малая чувствительность реальных траекторий роста на эволюционных участках к изменению характеристик внешней среды. Отражением этого является тот
38

факт, что оптимальные условия существования достаточно широки. В результате и модели для описания процессов развития также должны быть мало чувствительны к изменению их вида на эволюционных участках, тогда как при выходе на границы, где свойства развивающейся системы качественно меняются, чувствительность к исходной информации и виду модели должна резко повышаться.
Далее на основе исследования свойств решений типовых моделей определяются пределы их использования, а следовательно, и критические значения характеристик развивающейся системы, при достижении которых она меняет качество.
В последующих разделах на основе анализа количественной информации, связанной с развитием биологических систем, оценивается содержательность и конструктивность такого подхода.
3.1. Уравнение роста
Анализу основных тенденций роста биологических систем посвящено большое количество исследований [ J. S. Huxley, 1932; Шмаль-гаузен, 1935; Brody, 1945; Medawar, 1945; Bertalanffy, 1952; Вин-
берг и др., 1968].
Итог этих исследований подведен в работе Медавара [Medawar, 1945], в которой указывается, что о соотношении между размерами и возрастом организма можно сделать лишь одно общее фундаментальное утверждение, выражающееся уравнением:
dlnyldx=kf(z), (14)
где kf(x) — такое положительное число, что с увеличением х f(x) уменьшается, a df(x)/dx стремится к нулю. Кроме этого, отмечается особый случай, когда /(a:)=const, что соответствует экспоненциальному росту.
Зависимость [14] является достаточно общей, чтобы характеризовать практически любые непрерывные кривые роста. Свидетельством этого является тот факт, что при дробно-линейной функции
/(а)=г      ' + (15)
соотношение (14) соответствует уравнению Пирсона [Крамер, 1978], которое было использовано для нахождения общей зависимости, описывающей функции плотности распределения, наиболее, часто используемые для обработки экспериментальных данных. При этом Пирсон получил в качестве решений уравнения (14) ^ в зависимости от параметров функции (15) нормальное распределение, экспоненциальное, распределение Пуассона и т. д.
Шмальгаузен [1935] положил уравнение типа (14) в основу исследований, которые привели его к модели параболического роста.
Правая часть уравнения (14) в этом случае представляет собой Удельную скорость роста, т. е. прирост на единицу размера системы

за единицу временя. Шмальгаузен анализировал тенденции изменения удельных скоростей роста в функции возраста развивающейся системы, на основе чего по экспериментальным данным установил падение этой величины, обратно пропорциональное возрасту, т. е.
kf(x)=k/x. (16)
Такой метод построения модели роста представляется перспективным, так как позволяет на основе анализа динамики удельных скоростей роста определять вид функции f(x), а затем непосредственным интегрированием уравнения (14) получать кривую роста.
Однако использованные Пирсоном дробно-линейные функции (15) содержат не только участки, где функция f(x) убывает со временем или постоянна, как это указывал Медавар [1945], но и участки роста, на которых относительные приросты возрастают. Это соответствует возрастающим темпам экспоненциального роста. В процессах развития биологических систем экспериментально зафиксированы стадии роста, на которых величины относительных приростов могут возрастать. Эти стадии достаточно кратковремен-ны по сравнению с периодом развития системы. Наличие таких стадий нашло отражение в моделях динамики численности интенсивно растущих популяций [Бейли, 1970; Уатт, 1971].
В связи с изложенным представляется целесообразным рассмотреть общие свойства процессов, описываемых уравнением Пирсона, так как из него могут быть получены основные модели роста.
Покажем, что уравнение Пирсона является частным случаем уравнения с отклоняющимся аргументом вида
dy/dx = Xe-tey(z — x)1 * (17)
где у(х—х) — количественная характеристика состояния развивающейся системы (например, ее пространственные, временные или энергетические показатели), х — аргумент процесса, х — фазовое смещение переменной х (величина релаксации), К — параметр процесса.
В частности, когда аргумент х характеризует пространственную величину, то релаксация имеет смысл нелокальности, т. е. влияния соседних элементов на характеристики рассматриваемого.
Разложение правой части этого уравнения в ряд Тэйлора приводит после преобразований к уравнению
dy (х)   _    Хе** y(x)dx
а разложение экспонент в правой части с использованием двух членов разложения дает соотношение
dy (х)   _        + y(x)dx ~ 1 +   + А,2*2 •
40

При S~Ax получим *    •"'*      ;
dy(x) .. _ + _ "Ж + а;
или
(х)__а-\-Ьх
у (х) dx        (I2X2 -\- aix + а0 '
где
что соответствует уравнению Пирсона, т, е. функция, характеризующая относительные приросты / (ж), определяется дробно-линейной зависимостью типа (15). Таким образом, исследование свойств уравнения (17) позволит определить общие особенности решений достаточно широкого класса моделей роста.
Можно показать, что уравнение (17) действительно содержит информацию относительно большого количества законов распределения и кривых роста, лежащих в основе обработки экспериментальных данных.
Рассмотрим динамическую модель процесса с предысторией
достаточно общего вида
(amDm + a^D™-1 +'... + а0). у (х) +
+ (bpDp + b^D*'1 + ... + 60) у (х - Z) = 0,
где
a((i=0, l,2,...,m);    вД/=0,1,2, ...,р) -постоянные коэффициенты, D=d/dx — оператор дифференцирования. Будем искать решение этого уравнения в виде
y=ez*. (18)
Тогда получим характеристическое уравнение R (z) = amzm + am^1zm-1+ ... + Ж> +
+ {bvzv + b^zv-i + ... + b0) е-г* = 0. При нахождении решений этого уравнения для z>0 принимают
R(z) ~amzm+bpzpe-zlc=0.
Заметим, что основные законы распределения, используемые для обработки экспериментальных данных, соответствуют этому соотношению при конкретных значениях входящих в него параметров (табл. 1).
Дифференциальное уравнение, соответствующее данному характеристическому, имеет вид
у* (s)=_^Lp(P> m
41

Таблица 1. Параметры уравнения (17) и соответствующие им функция плотности вероятности

Параметры уравнения
Функции"плотности вероятности
Наименование закона распределения 

р—m
X  
1

0
1/252
l/3/2S.e-^2aJ
нормальное
2

0


экспонента
со
т\
к
1

Пуассона
4
1/С*Ц>-\)\
/с —1
. 1/е

гамма
5
к/с
fc —Г


Вейбулла
6

1
Х/2б*

Релея
7
1/2^2 Т(к/2)
(/с-2)/2
1/2
HM/2e-v2,

 


/2*/2Г (/с/2)

Примечание: о, к, С — параметры соответствующих распределений.
откуда с учетом (18)
zmezx=--v-zvezxe~zx.
(19)
Поделив обе части этого уравнения на zm 4, получаем
zez
: _ _     V   Xv~m+1ez(x-x)
и, снова учитывая (18), найдем
dy (х)
=--R-ZV-m+ly (х — х).
(20)
dx
Из (19)
Ьр =      егх
a 7V-m '
m z
и уравнение (20) перепишется в виде, соответствующем уравнению (17).
Таким образом, уравнение (17) является достаточно общим. Обозначим в уравнении (17)
Хех* = к,
тогда это уравнение перепишется в виде
dy/dx=ky (х—х). (21)
Назовем (21) уравнением воспроизводства. Более общий случай 42

уравнения (21) можно представить, если считать параметры к и х функциямивремени: dy/dx=k (х) у [х—х (х) ].
Уравнение (17) получило широкое распространение для описания процессов экспоненциального развития. Здесь к — константа скорости роста.
В уравнениях (17) и (21) наличие запаздывания приводит к необходимости задавать для интегрирования начальную функцию на интервале времени от начального момента до
у (х) =ф0 (х),    х0—х<х<х0. 1
Уравнение (21) представляет интерес для моделирования развития биологических систем с различных точек зрения. Прежде всего в связи с тем, что рост в постэмбриональный период начинается экспоненциальной фазой [Шмальгаузен, 1935; Bertalanffy, 1952; Брей, Уайт, 1959], что соответствует постоянству относительных приростов. Основное ограничение, приводящее к изменению темпов роста, связано с изменением внутренних и внешних условий роста.
Во всех основных моделях роста, используемых в настоящее время, описание развития биологических систем производится как функция только характеристик роста на данном уровне и возраста [Зотин, 1974; Мина, Клевезаль, 1976]. Такие модели называют марковскими или моделями без предыстории, т. е. характеристики системы в данный момент определяют только по состоянию системы на тот же момент вне зависимости от характера предыдущего развития. В то же время известно, что предшествующее развитие оказывает существенное влияние на состояние биологических систем и что процессы развития являются релаксационными.
Необходимость учета релаксационных свойств объектов, в том числе и биологических, давно обсуждается в литературе. Сущность влияния предыстории на развитие динамических процессов отмечал еще Пирсон [1911]:  «...мгновенная передача влияния тела Р, внезапно пришедшего в движение, на удаленное от него тело Q представляется мало правдоподобной; требуется некоторое время для того, чтобы изменение положения Р могло опираться на Q» (с. 637). Вольтерра [1976] при разработке математической теории борьбы за существование учитывал влияние последействия, или памяти на текущее состояние  системы. Гудвин [1966]  при разработке динамической теории внутриклеточных процессов указывал, что «переход от дифференциальных уравнений и интегральных инвариантов к функциональным уравнениям и инвариантным мерам сделал бы теорию значительно более мощной, позволив рассматривать гетерогенные системы вместо гомогенных, а также системы с запаздыванием и гистерезисными эффектами» (с. 32). В настоящее время направление, связанное с учетом влияния пространственных и временных предыстории на текущее состояние систем, получило настолько существенное развитие, что началась разработка механики сплошных сред, в основу которой

положена аксиома о существенном влиянии предыстории &а развитие процессов [Тру с дел л, 1975].
Уравнение (21) примечательно также тем, что его частйым случаем является экспоненциальный рост. Функцию запаздывающего аргумента у{х—х) для получения приближенного решения этого уравнения можно разложить в ряд Тэйлора с использованием двух членов разложения. Большего количества членов в этом случае брать нельзя, чтобы исключить появление колебаний, не свойственных основному исходному процессу [Эльсгольц, Нор-кин1971]:
y(x—x)=y(x)—xdy(x)/dx. (22)
Подставляя (22) в (21), получим
dy (x)/dx .
l + ie
(23)
Очевидно, что при постоянных к и х уравнение (23) совпадает с уравнением экспоненциального роста с постоянным темпом.
Процессы экспоненциального роста, соответствующие постоянным значениям относительных приростов, занимают важное место при моделировании развития биологических систем. Шмаль-гаузен [1935] указывал, что такие процессы характерны для роста при существенно постоянных внутренних и внешних условиях. В качестве примера системы, растущей при постоянных внутренних условиях, он приводил палочковидную бактерию Bacillus megatherium, питание которой происходит через поверхность, а площадь поверхности в связи с малой относительной толщиной бактерии пропорциональна длине, тогда как питать надо объем, в этих условиях также пропорциональный длине. Значит, питаемый через единицу поверхности объем не меняется в процессе роста палочковидной бактерии, в результате чего ее рост оказывается экспоненциальным (рис. 16).
Переменный характер внутренних и внешних условий развития приводит к тому, что при определенных характеристиках системы и среды, в которой она развивается, экспоненциальный рост с постоянным темпом прекращается. Тем не менее режимы экспоненциального роста оказываются в ряде случаев достаточно длительными, чтобы нельзя было не обратить на них внимание. Например, рост численности населения США в течение 60 лет с 1800 по 1860 г. оставался экспоненциальным (рис. 17).
Таким образом, экспоненциальные режимы роста соответствуют не только развитию при постоянных условиях, но и при изменении условий в определенных диапазонах, не превосходящих некоторых критических значений. Определение этих критических значений и представляет, по нашему мнению, одну из основных проблем моделирования развития, так как позволяет установить пределы использования результатов, характеризующих  сложив-
44

шуюся тенденцию роста. Броди [Brody, 1945] положил экспоненциальные на отдельных участках режимы роста в основу моделирования развития биологических систем, определяя в каждом случае пределы использования экспоненциальной зависимости но морфофункциональным признакам. Шмальгаузен [1935} отмечал, что процессы роста начинаются экспоненциальной фазой. Однако он полагал, что в дальнейшем экспоненциальный рост переходит в параболический. В связи с этим место экспоненциального роста в процессах развития биологических систем, их связь с процессами аллометрического типа и ранними стадиями развития представляет предмет дальнейшего рассмотрения.
Основная проблема, подлежащая исследованию,— являются ли_/ выделяемые по морфофункциональным признакам перестройки объективно связанными с классом используемой модели, т. е. можно ли, имея модель роста, определить момент развития, когда изменится характер роста системы.
Прежде всего с этой целью рассмотрим   свойства   решений уравнения воспроизводства (21).
Будем искать решение этого уравнения в виде (18), где
z=K+iu. Отсюда z — fce-**=0
или
X + iv — fce-^cos vx — i sin v%) = 0.
Разделяя действительную и мнимую части, получим
X — ке~х* соз их, v == ке~^ sinvx.
Обозначим
Хх = и;     voc = v;     кх = \х. (24)
Тогда система уравнений перепишется в виде и = \ie-u cos v, v = — p,e~usinv.
Режим роста будет соответствовать случаю, когда имеются корни системы уравнений (25) с действительными положительными частями и нулевыми коэффициентами при мнимой части. Если коэффициенты при мнимой части в данном случае будут отличны от нуля, то может реализовываться режим колебаний с экспоненциально возрастающей амплитудой. Такие режимы приводят к разрушению системы, поэтому будем считать их запретными для развития. В связи с изложенным представляет интерес
45

мнм 20 V
I   ,1
36
60 мин
о       '12       24
Рис. 16. Рост палочковидной бактерии Bacillus megatherium По абсциссе — возраст;  по ординате — длина,   шкала   логарифмическая   [По:. Шмальгаузен, 1935]
200 100 50
1900
1950 годы
Рис. 17. Рост численности населения США в 1800—1970 гг. По абсциссе — годы; по ординате — численность населения, шкала логарифмическая [По: Мендельсон, 1959; Народонаселение стран мира, 1978]
1,29 u
Рис, 18. Графическое решение уравнения и=3/2ле~м, определяющее границу области стабильного роста систем
По абсциссе — темп экспоненциального роста, умноженный на время запаздывания ; по ординате — значение функций левой и правой частей приведенного уравнения
исследование областей пространства параметров системы (25)* где и¥=0 и v=0. При этих условиях уравнения (25) приводят к соотношению
u=\ie-u. (26)
При и>0 здесь имеется единственный действительный положительный корень и при pi>0.
Для определения границы диапазона параметров уравнения (21), где возможна реализация единственного действительного положительного корня характеристического уравнения (26) при нулевом коэффициенте при мнимой части этого корня, что соответствует экспоненциальному росту, рассмотрим условия, при которых появляются корни с положительными действительными частями и отличными от нуля коэффициентами при мнимой части.
Из второго уравнения системы (25)
и=Ы\х+\п (— (sinv) /v). Это уравнение имеет смысл только в областях, где sin v<0, так как логарифм отрицательного числа не существует. Первый такой диапазон n<v<2n. Отсюда в области 0<v<2n заведомо не существует псевдоположительных корней характеристического уравнения, т. е. корней, в которых при положительной действительной части имеется ненулевой коэффициент при мнимой части. Положение этой границы и определяет величину диапазона, в котором отсутствуют псевдоположительные корни характеристического уравнения.
Нулевое значение действительной части корня характеристического уравнения получается из уравнений (25) при выпол-
46

нении условии
cosv=0,    [a=-v/sinv,
откуда
v=e/,*  и   ul=**/2k. (27)
Из (24) и (25) определим предельное значение параметра и, соответствующее единственному положительному корню характеристического уравнения и=3/2ле~"и, откуда численное или графическое решение (рис. 18) дает значение и=1,29, или с учетом обозначений (21)
ЯЖ=1,29. (28)
Полученная оценка границы существования режимов экспоненциального роста без колебаний с возрастающей амплитудой характеризует ограничение сверху на диапазон параметров уравнения (21), где возможен процесс экспоненциального роста с постоянным темпом или постоянство относительных приростов.
Выход на указанную границу при фиксированном значении темпов экспоненциального роста должен приводить к изменению темпов роста, либо величины запаздывания, обеспечивающих возвращение системы внутрь области, где могут реализовываться стабильные экспоненциальные режимы роста.
Рост амплитуды колебаний при выходе в критический диапазон, определяемый соотношением (28), фиксируется в экспериментальных данных в виде увеличения их разброса. Это общее свойство критических точек для процессов различной природы. Так, Паташинский и Покровский [1975] указывают, что по мере приближения к критической точке в веществе растут флюктуации и отмечают определенное сходство флюктуационных явлений в разнообразнейших критических точках. Урысон [1973] отмечает, что индивидуальные вариации размеров тела человека при постэмбриональном росте усиливаются в период увеличения скорости роста. В связи с этим «повышение изменчивости морфологиче^ ских признаков в возрастах, относящихся к периоду полового созревания, может быть объяснено определенными биологическими факторами, в частности влиянием различий в физиологическом возрасте в период прохождения индивидуумами через критические точки» (с. 26). Такую же точку зрения высказывал Д'Арси Томпсон [Thompson, 1945].
В критических точках наблюдается рост восприимчивости системы к внешним воздействиям [Паташинский, Покровский, 1975] и возникновение «дальнего порядка» — не исчезающей в макроскопических масштабах корреляции между элементами системы [Фишер, 1968]. Указанные особенности соответствуют критическому характеру границы (26). Иллюстрацией процесса, про

м
1fi\
% 60
40 20
"4   +'
WW
10
15 16 Годы
Рис. 19. Динамика роста мальчиков (а) и его коэффициента вариации (б) По абсциссе возраст; по ординате: а — длина тела, б— коэффициент вариации, %; по ординатам шкала логарифмическая [По: Куршакова. 1973]
исходящего при переходе системы через критические точки, когда между критическими точками рост является экспоненциальным, может служить динамика роста мальчиков и соответствующая кривая изменения коэффициентов вариации [по: Урысон, 1973]
где т — математическое ожидание (среднее арифметическое)» о — среднее квадратическое отклонение. На рис. 19 видно, что моментам изменения тенденций (темпов экспоненциального роста) соответствуют максимумы коэффициентов вариации, характеризующих величины разброса результатов измерений. Отметим принципиальное сохранение средних значений коэффициентов вариации за весь период постэмбрионального роста и повышение их значений лишь в критических точках, соответствующих изменению темпов экспоненциального роста. Такой характер разбросов ростовых характеристик достаточно типичен.
Граница области экспоненциального роста (28) отделяет пространство параметров уравнения (21), где процесс развития идет достаточно стабильно, что находит отражение в понижении в этой области коэффициентов вариации. В то же время имеется тенден^ ция к росту в области, предельно близкой к границе (28), так как это обеспечивает максимальные темпы роста, в результате чего система может достигать большего размера. Экспериментально зафиксированы механизмы отбора особей на максимальные темпы роста [Роус С, Роус Ф., 1964]: у головастиков лягушки Rana pipiens . выделен гормон, угнетающий рост медленно развивающихся особей.
В связи с изложенным можно предположить, что, зная характер изменения времени запаздывания в процессе развития биологической системы, можно считать его, в соответствии с соотно-
48

шениём (28), определяющим величину относительного прироста dy(x)/y(x)dx=Xa/£. ~~
Отсюда линейный рост времени запаздывания с возрастом систе-мы х=Ах приводит к параболическому росту dy(x)/y(x)dx** —а/Ах, тогда как линейное убывание времени запаздывания о возрастом x~AJx — к росту по кривой нормального распределения (кривой Гаусса) dy (х)/у (х) dx=aAtx или, после интегрирования,
В данном случае параболический рост представляет собой огибающую экспоненциальных режимов роста при убывающих* а кривая Гаусса — при растущих темпах.
Открытым в данном случае остается вопрос о том, насколько результаты моделирования зависят от точности определения параметров зависимости времени запаздывания от размеров системы и ее возраста при заданном общем виде такой зависимости. Этим может быть решен вопрос о том, какой моделью достаточно пользоваться, чтобы определить основные тенденции роста и критические точки, выделяющие участки развития при условиях, близких к постоянным.
Для этого необходимо провести исследование зависимости характера траектории роста системы от вида правой частц динамического уравнения модели. Такое исследование может быть выполнено в рамках понятия о грубости системы [Андронов и др.., 1967; Баутин, Леонтович, 1976].
3.2. Условия грубости уравнения роста
В предыдущем разделе показано, что для моделирования развития биологических объектов необходима информация о влиянии внутренних и внешних условий на время запаздывания. Трудности определения подобного рода данных для широких классов объектов приводят к необходимости выработки обоснованных требований к точности исходных данных и методической точное сти моделей, соответствующих цели исследования.
Очевидно, что на различных участках кривой роста ценность информации, доставляемой каждым экспериментом, при равномерном их распределении в диапазоне времен развития, различна. Так, на участках эволюционного роста (например, экспоненциального) новые экспериментальные точки практически бесполезны, тогда как критические моменты в развитии системы, связанные с изменением ее морфофункциональных структур, требуют существенно более детальных исследований. При проведении экспериментов по росту целесообразна именно "такая стратегия исследования — максимальное количество экспериментов проводить в''79ЖрЬс.таостях' точек, соответствующих изменению темпов роста,
49

Дня оценкя влияния функциональных зависимостей времени запаздывания от внутренних и внешних условий на траекторию роста воспользуемся методами, разработанными для исследования грубости динамических систем [Андронов, Понтрягин, 1937; Андронов и др., 1967; Баутин, Леонтович, 1976]. Система считается грубой, если при изменении в некоторой области вида правых частей дифференциальных уравнений характер траектории развития сохраняется. Потере грубости соответствует изменение характера развития системы. С этой точки зрения потеря грубости системой будет соответствовать критической точке, разделяющей два принципиально различных по характеру развития процесса.
В экспериментальной биологии проверке грубости соответствует исследование реакции биологической системы на повреждающее воздействие. Низкая чувствительность к изменению величины действующего стимула соответствует грубой системе. При моделировании процесса развития системы в этой области чувствительность результатов моделирования к изменению вида правой части дифференциальных уравнений процесса в определенном диапазоне (при сохранении основной тенденции развития) также должна быть низкой. При повышенной чувствительности биологических систем к внешним воздействиям говорят о наличии порога [Насонов, 1959], критической точки [Светлов, 1960]. Этот факт и должен найти отражение в модели процесса. Модель, адекватная процессу, должна быть грубой на участках развития системы с малой чувствительностью к внешним воздействиям и терять грубость там, где процесс развития системы проходит критическую фазу.
Произведем оценку грубости уравнения воспроизводства (21) по А. А. Андронову. Из уравнения (21) за счет смещения аргумента получим
dy(x+x)/dx=ky(x). 1 (29)
Разложение левой части уравнения (29) в ряд Тэйлора с использованием двух членов разложения дает соотношение
dy(x)/ dx+xd2y(x) /dx2=ky (х)
или
d2y(x)/dx2=l/x[-dy(x)/dx+ky(x)]. (30)
Уравнение (30) эквивалентно системе двух дифференциальных уравнений
(x)/dx=w (х), dw {х) dx=—lx • w (х) +к/х- у{х).      (31)
Пусть
1/Ж=ф(1/). ^>
Тогда из системы уравнений (31)
dy(x)/dx=w(x)=p,   dw(x)/dx=-y(y)w(x)+kq>(y)y(x)=Q.   (33)
50

Необходимые условия грубости динамической системы (33) по Андронову и др. (1967)
Д =
ру
qy
или применительно к рассматриваемой системе 10 1
ф0,
откуда
— w + ку + к
ФЫ 1 фо % (у) J
или
dyw фк
(34)
dx
(35)
Таким образом, наличие экстремума функции ф(г/), определяемой в соответствии с (32) величиной, обратной времени запаздывания, является одним из необходимых условий грубости системы.
Появление экстремума у функции <р(у) можно наблюдать на данных о динамике интенсивности обмена, которая является характеристикой времени, затрачиваемого системой на обменные и регуляторные процессы, что и определяет время запаздывания.
Полученное условие грубости (34) показывает, что на участках роста, где зависимость запаздывания от размера системы не имеет экстремума, модель является грубой. В этом смысле, считая, что зависимость размера от возраста может быть представлена некоторой функцией, можно получить аналогичный вывод для запаздываний, зависящих от возраста системы. Таким образом, принципиальным является вопрос о том, имеет ли функция, характеризующая изменение запаздывания со временем или размером системы, экстремум на анализируемом участке развития. Наличие экстремумов в таких функциях является экспериментально зафиксированным фактом.
Установлено [Зотин, 1974], что у зародышей лососей в ходе развития происходит медленное усиление интенсивности дыхания, которое достигает максимума к стадии гаструляции, а затем снижается (рис. 20). Как уже упоминалось, на критический характер стадии гаструляции указывал Светлов [1978].
В связи с тем, что на участках монотонного изменения времени запаздывания с возрастом или размером системы модель является грубой, если выполняется неравенство (35), требования к точно-
51

тЩ/чйс-т
OA
Ft*  t.    \-
10
20        30 Сушки
Рис. 20. Интенсивность дыхания изолированных от желтка зародышей лосося на ранних стадиях развития при 5°
По абсциссе — возраст; по ординате — интенсивность дыхания, ммЮ2/час мг сухой массы, б — равняя бластула, s — начало гаструляции, п — стадия половинного обрастания желтка [По: Зотин, 1974]
2 г кг
1
Ofi
0,4 01
+ +
I г » | | f i t i f г г i i г г i,-
кз
1
1
0,8 0,4
0,1
А
4
.1
10
15
20
Годы
i    ' 1 | mil ||""""
5      10   15 20 Лет
Рис. 21. Динамика изменения массы камбалы с возрастом
По абсциссе — возраст; по ординате — масса; шкала логарифмическая [По: Мина,
Клевезаль, 1976]
Рис. 22. Динамика изменения массы камбалы с возрастом По абсциссе — шкала логарифмическая
сти информации о конкретном виде функции ср(у) и точных значениях параметров этой функции будут низкими. В связи с этим может оказаться, что широкий класс моделей, которые характеризуются одной и той же тенденцией изменения функции <р(у) до критической точки, будет практически одинаково хорошо описывать результаты измерений характеристик роста.
Практическая реализация указанного пути исследования является в некотором смысле возвратом к исследованиям, которые проводились в 1920—1930 гг., но в дальнейшем не получили развития [Шмальгаузен, 1935; Brody, 1945; и др.]. Исключением являются работы Терсковой [1975, 1978], где параболический рост применен для определения критических стадий развития эмбрионов птиц.
Сложившаяся в настоящее время оценка этого направления представлена соображениями типа приведенных Миной и Клевезаль [1976]: «Поскольку уравнение экспоненциальной зависимо
52

сти, как правило, не может описать рост животных в течение жизни, мы этот случай рассматривать не будем» (с. 19). Далее анализируются модели, претендующие на описание роста животных лв течение всей жизни»,— логистическая кривая, кривая Гом-пертца, уравнение Берталанфи, которые даже в рассмотренных там же примерах оказываются не более содержательными, чём экспоненциальная и параболическая модели.
Так, на основе обработки данных о постэмбриональном росте камбалы Мина и Клевезаль делают вывод, что логистическая функция не подходит, кривая Гомпертца описывает процесс роста только с 2 до 8 лет, а уравнение Берталанфи — с 2 до 16 лет. На рис. 21 приведена экспоненциальная зависимость, построенная но данным Бивертона и Холта [Beverton, Holt, 1957], из которой видно, что она хорошо описывает процесс с 2 до 8 лет, т. е. на том же отрезке, где удовлетворительно работает кривая Гомпертца. На рис. 22 те же данные приведены в двойном логарифмическом масштабе, откуда видно, что с 2 до 17 лет процесс может быть представлен параболической моделью вместо уравнения Берталанфи.
Таким образом, в данном случае уравнение Гомпертца является в области грубости аналогом модели экспоненциального роста, а уравнение Берталанфи — параболического роста. Несомненным преимуществом экспоненциальной и параболической модели в данном случае является меньшее число параметров, которые необходимо определить по экспериментальным данным.
К разбору конкретных приложений свойств грубых моделей мы вернемся при анализе закономерностей развития биологических систем.
На основании полученных результатов, требующих для обеспечения грубости системы монотонного изменения времени запаздывания, выделим в кривых роста следующие основные этапы, соответствующие различным тенденциям изменения относительных приростов:
— относительные приросты растут с возрастом или размером системы [Бейли, 1970; Уатт, 1971];
— относительные приросты постоянны, что соответствует экспоненциальному росту [Шмальгаузен, 1935; Brody, 1945; Medawar, 1945];
— относительные приросты убывают с возрастом или размером системы [Шмальгаузен, 1935; Medawar, 1945; и др.].
В рамках каждого этапа изменение вида зависимости времен запаздывания от возраста и размера системы может рассматриваться как возмущение относительно вида правой части дифференциального уравнения [31], при этом запас грубости системы будет определять пределы изменения вида функций <p(i/), которые практически не изменят характера результатов модели развития. Экспериментальные проверки показывают, что практически достаточно знания качественного характера тенденции изменения величины времени запаздывания с возрастом и размерами системы,
53

хм *Х у
после чего достаточно использовать линейное приближение для функции <р(у).
Характер зависимости ф(г/) показывает, что имеет место взаимодействие двух факторов, противоположно действующих на время запаздывания: один из них уменьшает это время, а другой — увеличивает. В принципе время запаздывания, как и время цикла функционирования любой развивающейся системы, состоит из двух основных составляющих — затрат времени на использование системой материальных* потоков, которое тем меньше, чем больше размер системы, и затрат времени на управление процессом развития, которое растет с ростом размера системы. В связи с этим х—Хм+Ху, где х — время цикла функционирования системы, хм — время на использование материальных потоков, ху — время на управление.
Тогда хм~1/у и ху~у. Отсюда
x=ajy+a2y. В начале развития преобладают затраты времени на усвоение материальных потоков, размер системы небольшой и управленческие проблемы не являются определяющими. При достижении некоторого критического размера затраты времени на управление становятся больше затрат на утилизацию материального потока Ху*>хм, в результате чего с ростом размеров системы растет время запаздывания и лимитирующим становится время на управление.
Точка, где хм=хУ1 соответствует минимуму кривой (рис. 23), следовательно,  является  точкой  потери грубости. Отсюда для описания процесса развития достаточно рассмотреть три указанных выше случая.
В связи с грубостью уравнения воспроизводства исследование моделей экспоненциального и аллометрического роста, а также модели, соответствующей кривой Гаусса, будет описывать достаточно общий случай развития биологических систем. Основной проблемой остается определение критических точек, соответствующих переходам от одного режима развития к другому, т. е. пределов использования указанных видов зависимостей.
Определение пределов использования результатов любой модели представляет проблему первостепенной важности, так как позволяет определить диапазоны параметров системы и ее характеристик, соответствующие сохранению механизма, положенного в основу моделирования. В конечном итоге моделированию подлежат критические уровни, достижение которых приводит к морфофункциональным и организационно-структурным перестройкам: развивающейся системы.
Рис. 23. Изменение величины релаксаций с размером системы
По абсциссе — размер, х; по ординате — релаксация, х
54

Основные направления, в которых ведется исследование критических уровней, следующие:
—■ экстраполяция до предела (в частности, до абсурда), когда определяется время достижения бесконечного или нулевого размера [Бейли, 1970; Уатт, 1971] либо асимптотического значения, как при определении асимптоты логистического роста [де Вит, 1964]:
— исследование качественных свойств динамических систем {Молчанов, 1971а, 19716, 1973; Чернавский, Чернавская, 1972];
— теория размерности [Мигдал, 1975];
— теория катастроф [Том, 1970; Брёкер, Ландер, 1977].
Информационная значимость результатов, получаемых с использованием этих направлений, принципиально разная. Так, в случае экстраполяции до предела определяется верхняя оценка для используемой модели, при этом оценить степень близости предельного значения к реализуемому, как правило, невозможно. Однако получаемые в этом случае оценки часто оказываются достаточно хорошими практически, так как при приближении к предельному значению функции обычно меняются очень резко, в связи с чем ошибки по времени относительно истинного значения времени изменения тенденции невелики.
При исследовании качественных свойств динамических систем, как и в теории катастроф, используют нелинейные дифференциальные уравнения такого вида, что они содержат бифуркационные границы. В ряде случаев качественная картина, характеризующая положение бифуркационных границ, известна [Молчанов, 1971а], и подбирается такой вид правых частей дифференциальных уравнений, который соответствует этой фазовой картине. При этом растет количество параметров, которые необходимо определять. В то же время остается существенный произвол в выборе нелинейных уравнений, тем более, что первичные представления о фазовом портрете системы могут быть неполными.
Критические значения характеристик системы с использованием теории размерности основаны на эмпирическом правиле, что безразмерная переменная в критической точке — порядка единицы [Мигдал, 1975]. Это правило широко используется для получения оценок предельных значений. Однако его теоретическое обоснование отсутствует, в связи с чем остается неясной природа высокого качества получаемых результатов.
Целью приводимого ниже исследования является разработка подхода к моделированию, позволяющему обходиться минимальным количеством параметров, которые необходимо идентифицировать по экспериментальным данным для определения траектории развития системы на эволюционном участке и момента перехода к развитию в соответствии с другой моделью, что определяется изменением лимитирующих развитие условий.
Нам представляется, что три выделенные группы уравнений, являющиеся грубыми для соответствующих участков развития, могут представить интерес с этой точки зрения, если для них определить условия, соответствующие критическим периодам развития.

3.3. Критические уровни моделей экспоненциального  развития
Известно, что при постоянных внутренних и внешних условиях развитие происходит с экспоненциальным изменением характеристик системы. Выше показано, что при изменении характеристик развивающейся системы в некотором диапазоне внутренних и внешних условий, развитие также остается экспоненциальным, если выполняются ограничения на параметры процесса, соответствующего уравнению (21).
Для уравнения (21) условие, выполнение которого обеспечивает экспоненциальный характер развития при &>0, т. е. для процессов устойчивого типа, согласно известным результатам [Эльс-гольц, Норкин, 1971], .
kx<i/e (36)
приводит к соотношению |Х|#<1, где Я — величина текущих относительных приростов i/, соответствующая темпу экспоненциального роста.
Показано [Эшби, 1962; Турецкий, 1974], что для развития систем с запаздыванием характерно включение механизмов регуляции (адаптации к новым условиям развития) с периодом, равным времени запаздывания. Это соответствует известному способу оценки критических состояний, при котором в критической точке величина релаксации близка к возрасту системы [Паташинский, Покровский, 1975], т. е. х=х*, где х* — значение аргумента в критической точке.
Отсюда предельное соотношение
(37)
Величина Хх является безразмерной. Мигдал [1975] приводит широко используемое эмпирическое правило, согласно которому в окрестности критической точки величина безразмерной переменной — порядка единицы, что подтверждает справедливость последнего соотношения.
При выполнении соотношения (36) уравнение (21) имеет экспоненциальное решение
у/Уо=е-кх,
где г/о — состояние системы при х=х0. Тогда с учетом (37)
t/*/i/o=e-\ (38)
где у* — характеристика системы в критической точке.
Динамику релаксаций часто считают экспоненциальной функцией времени [Дей, 1974]. Для процессов с экспоненциальным изменением релаксаций со временем x=x0eh' *, откуда в соответствии с (38) х*/х0=е и из условия равенства величины релаксации аргументу развития процесса в критической точке
х*/х0=е. (39)
56

Отсюда соотношение, характеризующее последовательность критических значений аргументов, можно представить в виде
хк = еа*-\   А; = 1,2,... (40)
Если аргументом является возраст системы, значение Тк характеризует критический возраст (или критический рубеж). Таким образом, для экспоненциальных процессов развития систем критические значения аргументов (и функций) отстоят друг от друга так, что их соотношение постоянно и равно числу Непера, основанию натуральных логарифмов е=2,71828...
Число Непера как критическое хорошо известно и широко используется в теории и практике моделирования. Так, при моделировании роста с использованием уравнений Берталанфи и Гомпертца отношение предельного размера к размеру в точке перегиба равно е [Мина, Клевезаль, 1976]. Содержательность этих моделей, очевидно, связана с тем, что в них критический характер соотношения размеров развивающейся системы в е раз представлен самим видом используемой зависимости.
«Постоянные времени» вводят как величины, соответствующие интервалу времени, когда изучаемая характеристика убывает в е раз [Гинкин, 1962]. Так, при моделировании заболеваний расчет количества антител ведут по времени, за которое это количество убывает за счет старения в е раз [Марчук, 1979].
В теории информации известно, что максимуму помехоустойчивости соответствует величина вероятности, равная е-1 [Шеннон, 1963]. Этот результат согласуется с общим свойством функций типа энтропии, которые называют степенно-показательными [Саве-лов, 1960]. Энтропия как мера видового разнообразия использована для определения числа видов в коллекции с равнообильными видами [Песенко, 1978], в результате чего можно ожидать, что информационная статистика для такой модели зафиксирует соответствующее критическое соотношение.
Еще одна область, где величина е фигурирует как критическая, относится к построению иерархических структур из большого числа однородных элементов. Идеальной по времени распространения информации является иерархия с модулем, равным е [Флейшман, 1971]. В целочисленном варианте это соответствует троичной структуре, находящей широкое распространение в линейных иерархиях [Обэр-Крие, 1973].
3.4. Критические уровни аллометрического развития
Для моделирования процессов развития на стадиях (этапах), где относительные приросты убывают, используют степенные зависимости, которые получаются как огибающие последовательных стадий экспоненциального развития. На рис. 4, где представлены данные о весовом росте эмбриона курицы [Шмальгаузен, 1935], Шкала по ординате логарифмическая, в результате линейные участки роста представляют собой экспоненты. Здесь видны три
57

Рис 24 Изменение массы куриного эмбриона с возрастом По абсциссе-возраст; по ординате - масса. Шкала по обеим осям - логарифмическая [По: Шмальгаузен, 1935]
Рис. 25. Изменение массы леща в постэмбриональном развитии
По абсциссе-возраст; по ординате - масса. Масштаб по обеим осям - логарифми
ческий [По: Шмальгаузен, 1935]
последовательных экспоненциальных режима с убывающими темпами. Те же данные на рис. 24 представлены в двойных логарифмических координатах; в результате полученная линейная зависимость представляет собой параболу, т. е. степенную зависимость, которая линеаризует всю совокупность экспериментальных точек.
Шмальгаузен [1935] применил для моделирования процессов роста степенную зависимость, назвав ее параболическим ростом. Он показал, что параболический рост, так же как и экспоненциальный, в определенные моменты времени скачком меняет тенденцию, причем такие скачки оказались связанными с наиболее принципиальными перестройками в развитии животного. Например, для леща это моменты наступления половой зрелости и старости (рис. 25).
Розен [1969] отмечает, что существует большое число разнообразных функционалов, заданных на самых различных множествах объектов, как органических, так и неорганических, значения которых связаны уравнением вида
Y = a.X\ (41)
58

Выражаемую этим уравнением закономерность, систематическое изучение которой на примере живых организмов было впервые предпринято Дж. Гексли [J. S. Huxley, 1932], называют законом гетерогонии, или законом неравномерного (аллометрического) роста. Это единственный вид явной зависимости между функционалами, относящимися к биологическим объектам, который был подвергнут систематическому изучению.
Под термином «аллометрическое развитие» мы будем понимать класс функционалов, описывающих соотношение между двумя характеристиками процесса, соответствующий соотношению (41) вне зависимости от того, каков конкретный смысл аргумента.
Использование принятого Шмальгаузеном [1935] термина «параболический рост» в данном случае невозможно, так как будет рассматриваться общий случай, когда показатель степени может быть не только положительной величиной. Кроме этого, получаемые результаты относятся к общим свойствам процессов, модель которых имеет вид (41), в результате чего они могут представить интерес не только для анализа динамики развития во времени, но и для решения задач относительного роста частей биологической системы.
Примечательно, что в последнее время обращается внимание на принципиальную общность моделей аллометрического типа для процессов различной природы. Показательными в данном случае являются статьи Саважо [Savageau, 1979], одна из которых («Ал-лометрический морфогенез сложных систем: вывод основных уравнений из первых принципов») включает как рубрики организ-менную и популяционную биологию, социальные системы, экономику и организацию. Автор с удивлением отмечает, насколько большой ряд явлений описывается этим простым уравнением: соотношение характеристик у основных групп животных и высших растений, работы по морфологии, физиологии, фармакологии, биохимии, цитологии, эволюции, анализ данных об этиологии ряда заболеваний.
Длительное время аллометрия рассматривалась как эмпирический закон, однако появилось большое количество работ, где делаются попытки ее теоретического обоснования. Прежде всего это результаты, связанные с теорией подобия и размерности [Седов, 1977]. Уравнение (41) замечательно тем, что при его использовании отношение значений двух производных величин не зависит от масштаба основных единиц измерения. В связи с этим безразмерные критериальные зависимости в теории размерности используют как основное именно это соотношение.
Большие исследования в теории ветвящихся случайных процессов привели к выводу закона Ципфа [Zipf, 1949], который считают одним из основных эмпирических законов современного науковедения. Рассматривается набор N элементов, каждый из которых помечен меткой из некоторого множества. Число различных меток п (х, N), каждая из которых встречается ровно х раз в выборке из N элементов при достаточно большом N, п(х, N) =
59

«Л/л:7, где А — константа, определяемая объемом выборки, у — показатель закона Ципфа.
Этот закон носит имена Ципфа, Лотки, Мандельброта, Паре-то—в зависимости от области применения [Яблонский, 1978]. Закон Ципфа, как и нормальный закон, является устойчивым (композиция устойчивых распределений дает распределение того же вида).
Розен [1969] выводил уравнение аллометрического роста из принципа оптимальности, анализируя класс функционалов, для которых минималь типа уравнений Эйлера дает степенную функцию.
Сава'жо [Savageau, 1979] предложил выводить аллометрию из дифференциального уравнения вида
±х = ахх\1 — Оъх\ \
считая, что на данном аллометрическом этапе величина xt лимитирует развитие системы.
Таким образом, большое количество эмпирических фактов и теоретические исследования (как связанные с моделированием /процессов аллометрического развития, так и посвященные исследованиям механизмов их столь широкой распространенности) обращают внимание на необходимость детального анализа свойств и возможностей применения аллометрической модели. Предметом нашего исследования является установление пределов, в которых аллометрическая модель может содержательно характеризовать
процесс развития.
Можно полагать, что границы использования аллометрической зависимости связаны с пределами развития системы в рамках определенной морфофункциональной структуры, границы которых определяются эндогенными и экзогенными условиями развития, В результате в определенные моменты происходят качественные изменения в объекте, например метаморфозы, ввод в действие новых функциональных систем, фазовые переходы.
Отсюда возникает необходимость выяснения пределов применимости аллометрических зависимостей с целью выявления узловых периодов развития, характеризующихся определенными механизмами их формирования, и классификации этих периодов.
Выше говорилось, что аллометрическое развитие соответствует убыванию относительных приростов обратно пропорционально возрасту системы. Отсюда дифференциальное уравнение таких стадий развития имеет вид
dy/dx=B/x-y. (42)
Аргумент х в уравнении (42) при изучении временной динамики представляет собой время или возраст. Однако аргументом может быть и не возрастная, а размерная, энергетическая и другая характеристика. В этом случае остается принципиально неясным, какая из характеристик системы причинно обусловливает
00

изменение другой, т. е. какая из переменных в рассматриваемом процессе является функцией, а какая аргументом.
Интегрирование уравнения (42) дает соотношение (41), которое широко используется при изучении динамики размеров тела в процессе роста [Брей, Уайт, 1959; Автандилов, 1973], в органо-метрии и морфометрии [Автандилов, 1973].
Рассмотрение зависимостей типа (41) для процессов самой разной природы показывает, что при определенных значениях параметров происходит скачкообразная смена их значений в зависимости (41) при сохранении ее общего вида. Для ряда процессов, по которым собран обширный экспериментальный материал, установлено, что при достижении определенных узловых (критических) точек система теряет устойчивость и скачкообразно переходит в новое устойчивое состояние, которое характеризуется принципиально другой природой развития процесса. Так происходит переход от ламинарного течения крови в кровеносных сосудах к турбулентному. Значит, описание процесса (системы) на определенном интервале изменения его параметров должно содержать информацию о потере устойчивости при достижении некоторого критического уровня (задачей является определение этого уровня) и о характере развития после потери устойчивости. Дифференциальное уравнение (42) не содержит бифуркационных границ, т. е., описывая процесс на отдельных участках, оно не содержит информации о критических уровнях развития системы, за пределами которых параметры процесса меняются. Это значит, что описание процесса надо расширить так, чтобы при сохранении аллометрического развития на отдельных участках обеспечивалось определение узловых (критических) уровней развития.
Таким образом, основной особенностью аллометрического развития является то, что величины входящих в зависимость (42) переменных относятся к одному моменту времени и, следовательно, не могут учитывать эффектов, связанных с влиянием предыстории развития на его текущие характеристики.
Введем следующие допущения.
1. В уравнении (42) скорость изменения изучаемой характеристики у пропорциональна ее величине при значении аргумента, смещенном относительно рассматриваемого на величину релаксации, т. е. dy/dx=B/x-y{x—x), где х — фазовое смещение переменной у (величина релаксации).
Заметим, что в соответствии с ограничениями на пределы использования моделей экспоненциального роста (28) в предельном случае для х~х к=к/х, что учитывает релаксационный характер экспоненциальных процессов в уравнении аллометрического развития как огибающей экспоненциальных режимов.
Процесс аллометрического развития характеризует следующий, более высокий по сравнению с экспоненциальным, уровень агрегирования, в связи с чем ему должен соответствовать собственный механизм «памяти» о предыстории, что и учитывается данным допущением. ,
61

2. Величина релаксации х является функцией от аргумента
3. При выборе переменных, характеризующих процессы, обычно используют условия выпуклости или устойчивости (Фишер, 1968). В связи с этим будем рассматривать уравнение устойчивого типа
dy/dx=-B/x-y (х-х). (43)
Основной задачей исследования является установление пределов использования степенных зависимостей.
Будем искать решение уравнения (43) в виде
у = е-М*)я. (44)
В связи с тем, что параметры В/х и х(х) являются переменными, темп изменения размера системы также ищется как функция времени. Из (44)
у = — [X (х) X + К (х)] е-^х (45)
ж подстановка (44) и (45) в (43) приводит к уравнению
к(х)х + Х(х)=^е^х\
если принять
К(х—х)=Х(х). Тогда
Х(х)==^-^—1(х) + ~ eM»*J и (46)
Дифференциальное уравнение (46) определяет динамику темпов изменения размеров системы Х(х), если известна величина темпа в некоторой точке Xii lk(xi)==!kl. Обозначим
\х=\(х)х/В. (47)
Тогда после подстановки (47) в (46) получим
Выше указывалось, что для аллометрии характерно убывание темпов развития, т. е. к(х)<0. При К(х)=0 получим экспоненциальную модель, а при %(х)>0 развитие будет происходить с
растущими темпами.
В данном случае мы рассматриваем область темпов, в которой реализуется аллометрическое развитие, либо как его предельный
62

случай — экспоненциальное изменение размера системы, т. е.
£(х)<0. (49)
Из уравнения (48) условию (49) соответствует неравенство
вх(х)
— \х + е  *   *<0. (50)
Это соотношение выполняется, если
Вх(х)/х<1/е. (51)
Для граничного случая
в2(х)/х=1/е (52)
и
ц*=е. (53)
При изменении знака неравенства (51) %{х) меняет знак, т. е. режим падения темпов, соответствующий аллометрическому развитию, сменяется режимом роста темпов, что для аллометрии невозможно.
Таким образом, соотношение (12), обеспечивающее вместе с (13) i(#)=0, характеризует переход к экспоненциальному развитию и прекращение аллометрического развития.
При построении оценок решений дифференциального уравнения его заключают между двумя функциями, поведение которых известно [Чаплыгин, 1935; Брейн де, 1961]. Пусть Х(х) — решение дифференциального уравнения первого порядка
i=*F(x, Я),    (а<я<Ь), a Xi (х) — функция, удовлетворяющая условиям
\l(x)<F(x,Ki(x)) (a<x<b,        %i(a)<X(a)).
Тогда Xi(x)<X(x) (a<x<b). Отсюда условие (48) приводит к неравенству
— \i + e х М'< — (л + е е.
Тогда Xi>X, т. е. темп убывания величины % в уравнении (46) для оценки (53) является максимальным. Это обеспечивает при
y(x0)>yi(xo),   y(x)>yi{x).
Таким образом, из (53) величина \1*=Хх*/В=е определяет предельную оценку снизу для решений уравнения (43)

или
JL<J^L = e-*e. (55)
Уо ^ Уо v    '
В общем случае аллометрический режим развития обеспечивается при
(ХВх/х<1/е. (56)
Левая часть неравенства (56) соответствует уравнению (43) без запаздывания. В этом случае траектория процесса характеризуется степенной функцией (41). Из уравнения (46) видно, что при 3(х)=0 величина
М*) » JL < JL Г   я, (#) + — ем*)*<*>
х #2 ^ х  L х
в результате при Х(0)=Х0 уравнение без запаздывания дает оценку решений (43) сверху — уравнение (41).
Соотношение (55) характеризует оценку решений уравнения (43) снизу при фиксированной точке начала развития у(хо)=уо для верхней (41) и нижней (55) оценочных траекторий. П   = (**\-*< У*н =cJB
Уо        \ х0)    ^   Уо
откуда
**1*о <е* = 15,15426... (57)
Следовательно, соотношение между значениями аргумента, соответствующими началу и концу устойчивого аллометрического развития, является величиной постоянной и равной ее [Кузьмин, Жирмунский, 1980 а, б, г].
Исследуем устойчивость уравнения (43).
Уравнение (43) является линейным по г/, в связи с чем однородная часть его уравнения возмущенного движения совпадает € исходным уравнением.
Решение ^}Х\ уравнения (43) называется устойчивым, если для каждого е>0 существует такое 6(e)>0, что из "неравенства \(р(х)—^)(х) |<б(е) на начальном множестве следует
у(х) — у(х)\ ф (х)      (х)
при х<:Х0, где г|) (х) — любая непрерывная начальная функция [Эльсгольц, Норкин, 1971]. Не обладающие этим свойством решения называются неустойчивыми.
Исследование условий устойчивости тривиального решения уравнения (43) проведем на основе идеи Н. Н. Красовского [1959],
64

предложившего рассматривать функционалы, обладающие свойствами, аналогичными свойствам функций Ляпунова.
Тривиальное решение уравнения (43) устойчиво, если существует непрерывный определенно-положительный функционал
V[y(x),x]    (5?<$<0), snp М*)|<#, #>0, F (0, *) = <>,
производная которого вдоль интегральных кривых неположительна
-Тх-
где y9(x+s) — решение уравнения (43), определяемое начальной функцией
ф(ж)  (х0—Х<Х<Хо).
Рассмотрим функционал
о
V [у (s), х] = г/2 (х) + 2а I у* (х + s)ds =
—X
х
(58)
X—X
ГДеПр^;одГяТункЦионала (58) в силу уравнения (43) dV = - lS-у (х) у (х-£) + 2а[у*(х) -у*(х- *)].
dx х
или
iL = _ 2 \^-у (х) у (х-+    1У\«- *) ~ У2 (*)]} •
(-59)
(60)
Для того чтобы <2F/d#<0, достаточно
у2(х-х)-у2(х)>0. Это требование невозрастающего  характера  решения  уравнения (43).
Таким образом, условие (60) обеспечивает устойчивость тривиального решения в соответствии с приведенной выше теоремой.
Этот результат соответствует также теореме 11 Н. Н. Красов-ского [1968] об асимптотической устойчивости решений.
3  А. В. Жирмунский, В. И. Кузьмин 65

Устойчивое решение у w   называется асимптотически устой-<р(з)
чивым, если
Hm |— у(х)\ = 0
ф(*)   ФИ
для любой непрерывной начальной функции ^(2:), удовлетворяющей при достаточно малом б4>0 условию
\<р(х)-1р(х)\<01.
В соответствии с теоремой Н. Н. Красовского об асимптотической устойчивости решение уравнения (43) асимптотически устойчиво, если существует функционал V[y(s), х], удовлетворяющий условиям:
0 V
■ Yly(8),z]<Wily(z)] + Wt([l y2(s)ds] ') V[y(s),x]^W(y(x)),
где Wi(r) и Т^2(г) — монотонно возрастающие функции при г>09 причем Wi(0)=W2{0)=0, W(r) и ф(г) — непрерывные, положительные при г>0 функции.
Из приведенного выше функционала V[y(s),x] первые два условия теоремы выполняются:
0 7
l)Vly(*),z]<r* + cp*, т = у\ p=[ly*(s)ds] \    с>2а.
2) V[y(s),x]>r\
Из условия 3
dV/dx<—ф(г),    ф(г)>0   при   г¥=0.
Условие (60) определяет, когда dV/dx^O, что соответствует предельному случаю. В качестве ф(г) можно принять г2.
Таким образом, тривиальное решение уравнения (43) является устойчивым, если
y(x-2)-y{z)>0, (61)
т. е. величина у по траектории убывает.
При В>0 уравнение (43) представляет собой уравнение неустойчивого типа, для которого имеет место ограничение типа (52), разнящееся только величиной константы в правой части. Точное значение константы для произвольных переменных функций запаздывания в настоящее время неизвестно, хотя есть соображения
66

[Мышкис, 1972], что (В/х*) inf я <3,44.
Получаемая в этом случае оценка соотношения между критическими значениями аргумента не превосходит оценки, получаемой для уравнения устойчивого типа. В связи с этим можно предположить, что оценка (57) является предельной для степенных функций и, следовательно, должна характеризовать наиболее принципиальные ограничения.
Естественно, что внутри диапазона, определяемого соотношением (57), будут иметь место и последовательности критических соотношений меньшего масштаба. В то же время последовательные аллометрические режимы могут агрегироваться в более крупные критические уровни. Связи между критическими соотношениями для экспоненциальных и аллометрических режимов будут рассмотрены в разделе 3.6.
Таким образом, полученная величина соотношения между критическими значениями аргумента является одной из иерархии критических соотношений, которая, однако, должна встречаться достаточно часто в реальной действительности, по крайней мере там, где применимы степенные зависимости.
Соотношение (57) означает, что иерархически связанные процессы развития имеют среди критических параметров такие узловые параметры, последовательные значения которых различаются в ее раз.
Содержательность полученного соотношения можно оценить только на основе детального анализа большого экспериментального материала по процессам аллометрического развития. Конкретные материалы, связанные с проявлением критического соотношения (57), рассматриваются в разделах 4—7, а связям этого соотношения и критического соотношения для экспоненциальных процессов (в е раз) посвящен раздел 7.
Соотношение ее как критическое может быть получено при использовании теории размерности применительно к ряду частных задач, встречающихся в теории надежности, методах планирования экспериментов, кинетике роста биологических систем.
Так, в теории надежности используется двойной показательный закон с функцией распределения вида
F(x)
Этот закон применяется в задачах горячего резервирования и аналогичных им, в частности при схеме накапливающихся повреждений в весьма большом количестве элементов. Система состоит из большого количества параллельно соединенных элементов и функционирует до тех пор, пока остается хотя бы один исправный элемент. Если для критической точки х=1 [Мигдал, 1975], то получим величину вероятности [Кордонский, 1963]
3* 67

В теории планирования экспериментов введена преобразующая функция [Адлер и др. 1976]
d=l — е-е'
которую считают предпочтительной в связи с ее более высокой чувствительностью в средней зоне по сравнению с зонами, близкими к 0 и 1. Эта функция используется при поиске обобщенных показателей, характеризующих психофизическую деятельность [Сердюк, 1977].
Таким образом, проявляется интерес к использованию величины ее как самостоятельного основания логарифмического исчисления.
Изучение пределов использования зависимости типа (41), направленное на расширение возможностей анализа длительности стадий развития сложных систем, подчиняющихся аллометрической зависимости (41), может иметь ряд достаточно общих приложений, связанных со свойствами процессов, поведение которых описывается степенными функциями.
3.5. Критические уровни уравнения Гаусса
Выше показано, что случаю роста с возрастанием относительных приростов соответствует уравнение Гаусса
У = Уое^1\ (62)
В соответствии с приведенными в табл. 1 данными дифференциальное уравнение процесса в этом случае имеет вид
^L = Xe^y(x-ak),   а = -±?. (63)
Как отмечалось в разделе 2.1, при Я>0 область, где характеристическое уравнение для (63) имеет единственный положительный корень, определяется соотношением
w   —     3 2<j2 е     ^2    '
тогда
я<1,60о. (64)
Переход через границу (64) приводит к появлению псевдоположительных корней характеристического уравнения. Это указывает на принципиальное изменение характера протекающего процесса, в связи с чем при нарушении указанного соотношения необходимо использовать для описания процесса другую зависимость.
Из таблиц функции нормального распределения можно определить долю общего размера системы, принадлежащую диапазону,
68

соответствующему неравенству (64): /*т0=г0,9452, т. е. вне диапазона остается 5,48% общего диапазона развития системы. Отсюда соотношение 0,9452 : 0,0548=17,2 задает критический диапазон изменения характеристик процесса, который достаточно близок к критическому соотношению для аллометрического развития ее.
При Я<0 характеристическое уравнение в предельном случае приводит к соотношению и
Яа/2о2е so* < 1/е
и Я72о2-^0,28. Тогда
Я<0,75о, . (65)
и соответствующее табличное значение функции распределения
Fo=0,7734, (66)
т. е. вне диапазона (65) остается 23% длительности развития системы. Переход через границу (65) также приводит к качественному изменению характера процесса.
Выделение из общего количества измерений совокупности, содержащей ~1/4 общего количества, также достаточно характерно для имеющихся массивов экспериментальных данных. К этому вопросу мы вернемся в следующем разделе.
3.6. Условия синхронизации экспоненциальных и аллометрических критических соотношений
Анализ процессов развития показывает, что существуют переломные точки — рубежи в развитии. Исследование свойств экспоненциальных процессов, проведенное в 2.2, показывает, что для них существуют соотношения между последовательными критическими значениями аргумента, при которых меняются качественные свойства процесса. Последовательные критические соотношения согласно (40) можно представить в виде
хк = ех*~\ .4=1,2,(67Х
где xh — значение аргумента в критической точке с номером к. Если аргументом является возраст развивающейся системы, значение характеризует критический возраст (рубеж). Таким образом, последовательные критические точки (рубежи) должны быть представлены в экспериментальных данных геометрической прогрессией.
При выводе критического соотношения (67) в 3.3 мы принимали, что для процессов развития систем с запаздыванием имеет место цикличность, заключающаяся в том, что моменты адаптационных перестроек (управляющих воздействий) следуют друг за другом с периодом, равным времени запаздывания. Это значит, что существует последовательность циклов, которая будет проявляться

в результатах развития, однако менее значимо, чем в критических точках.
Можно ожидать, что сила проявления рубежей существенно увеличивается в тех случаях, когда они совпадают с рубежами, принадлежащими к другим системам длительностей циклов.
Имея в виду сосуществование в процессах развития экспоненциального типа критических точек, составляющих геометрическую прогрессию (67), и рубежей, представленных арифметической прогрессией, определим моменты совпадения рубежей, формируемых этими последовательностями. В такие моменты сила проявления критического рубежа должна быть максимальной.
Определим условия, при которых в один из интервалов времени, границами которого являются последовательные неравномерные критические рубежи, определяемые геометрической прогрессией (67), будет укладываться целое число рубежей, разделенных равномерными отрезками времени.
хк — возраст границы, соответствующей одному из неравномерно чередующихся рубежей, и при этом до следующего критического рубежа укладывается ровно «п» равномерно чередующихся рубежей длительностью Th. В соответствии с соотношением возрастов между последовательными неравномерными рубежами (67) интервал времени между последовательными критическими рубежами будет соответствовать соотношению
Л* _ a*-i =.■ (1 — е"1) хк = пТк. (68;
Пусть один цикл длительностью Тк находится за рубежом хк, т. е. при значении возраста рубежа хк~^—Тк будет приходиться на критический рубеж #ft-T, где к—у — номер критического рубежа.
Тогда xh~x-Th={n-l)Th. В этом случае
Лг1 - Л-v = (п - 1) Т\ (69)
Из (68) и (69) получим систему уравнений вида
я* (в _ 1) = епТ\      хк (ev-i — 1) = е? (п — 1) Тк. (70)
Отсюда
е — 1   _ е     п
е^"1 - 1 п-\
или
п    _У-1(е-1) п — 1~   е^"1 — 1
70

Тогда
е-1 (71)
п~ e-2 + e^-i)
lim п= -—2,39.
Значит, возможные варианты целочисленных значений п соответствуют только тг=2 или п=1. Определим значение, обеспечивающее целочисленное значение п.
По смыслу введенных соотношений f>2 и является целым числом. При ^=2
П= е-2 +е-* =1'58'
при 7=3
»=ТЙТ^ = 2-01- (72)
Дальнейшее увеличение значений 7 до <*> не приводит к появлению нового целочисленного значения п. Значит, единственное значение *у, удовлетворяющее поставленным условиям, 7=3, тогда п=2.
Синхронизация через т критических рубежей приводит к соотношению
fr 1        ™b
x = П
1 — е
В зависимости от т/г множитель, характеризующий отношение возраста рубежа к длительности равномерного промежутка времени Tk, xh/nTh имеет значения: 1,58 при т=1; 1,15 при ттг=2; 1,05 при ттг=3 и 1,0 при т=оо.
Отсюда следует, что имеет смысл рассматривать положение синхронизированных рубежей для пг=1 (разобранный выше случай синхронизации соседних рубежей) и /тг=2 (синхронизация через один рубеж), так как все последующие соотношения соответствуют положению рубежей, определяемых длительностью цикла Т\ ввиду того, что
Нтж* = 7гГ*.
71

Следовательно, полный набор критических рубежей задается следующими соотношениями.
для 771=1
основные синхронизированные рубежи до ft—3-го порядка,
\ 1 — е"1      / в — 1     '
промежуточные несинхронизированные рубежи, определяемые промежутком времени Г\
для 771=2
= 2*=»-А-7*.
1 — е 2 е2 — 1
Легко показать, что возможность синхронизации с более дальними рубежами в данном случае отсутствует.
В результате получается два крайних случая, первый из которых соответствует возможности синхронизации критических рубежей неравномерной шкалы (67) через 3 уровня от хк до хк~'6, а второй рубежам, определяемым равномерной шкалой интервалов Тк.
Случай синхронизации рубежей через один от хк до хк~2 может получаться относительно некоторых промежуточных фаз равномерных диапазонов Тк, так как в одном процессе развития основные критические точки можно объединить для случаев формирования рубежей процессов аллометрического развития и рубежей экспоненциального развития, но одновременно они не могут быть объединены для всех трех случаев в связи с соотношением (71). Это значит, что данный случай может соответствовать появлению некоторых вспомогательных рубежей, которые могут иметь своеобразные проявления в процессах развития.
Перейдем теперь к рассмотрению рубежей, определяемых соотношениями (70).
Из (68) с учетом (72) положение критической границы определяется соотношением
я* = 2е(в-1)"1ГЯс = 3,16Гк. (73)
Значит, если известна длительность равномерно расположенных рубежей Т\ то положение соответствующих критических рубежей при введенных условиях согласования неравномерно и равномерно расположенных границ рубежей может быть рассчитано. При этом оказывается, что между двумя последовательными неравномерно расположенными рубежами укладывается 2 равномерных диапазона, а между вторым и четвертым — один. Это избавляет от необходимости рассматривать отдельно характеристики неравно
72

Таблица 2. Положение критических рубежей для сишфошюнрованнога диапазона
Наименование
Обозначение
Зависимость для вычисления
Основной к Промежуточный к Основной к-1 Основной Ä—3
xh xh~l
**=3,16Г*
хь_2ГА=1,16Г* зЬ-ЗГ^ОДбГ*
мерных и равномерных последовательностей рубежей. Положение критических границ представлено в табл. 2.
Соотношение между возрастами основных критических рубежей для синхронизированного диапазона в рассматриваемом случае
я*-з = ез_ 20. (74)
В связи с тем, что положение последовательных основных критических рубежей определяется соотношением (67), а положенно основного критического рубежа А-го порядка через длительность формирующего его равномерно рубежа ^ — соотношением (73) „ то после подстановки (73) в (67) получим
Тк=Тк1Ч, (75)
т. е. длительности последовательных значений членов, формирующих равномерно расположенные рубежи, соотносятся между собой так же, как и возраста неравномерно расположенных критических рубежей — в соответствии с геометрической прогрессией
вида zY^-^rVT7*-^.
Таким образом, общий диапазон, в котором происходит синхронизация равномерных рубежей с неравномерно расположенными, не выходит за соотношение возрастов рубежей, соотносящихся в е3 раз. Это свидетельствует,, что неравномерно расположенные последовательные критические рубежи, в свою очередь,, агрегируются в более крупные критические группы. Такой вывод находит подтверждение в имеющихся данных по периодизации ряда природных процессов, которые мы рассмотрим в следующих разделах.
Агрегирование характеристик экспоненциальных процессов в диапазоны, объединяемые соотношением между критическими границами типа е3, показывает, что здесь имеют место процессы, более длительные, чем те, что происходят при небольших диапазонах изменения внутренних и внешних условий. В связи с этим может оказаться существенным учет характера изменения условий развития. Такие процессы, как показано выше, описываются алло
73

тлетрической моделью, которая является огибающей экспоненциальных режимов развития при убывающих темпах.
Критическое соотношение для аллометрических процессов характеризует развитие, в ходе которого не меняется тенденция убывания относительных приростов, тогда как полная синхронизация основных критических рубежей экспоненциального типа происходит при соотношении между рубежами ее.
Если большее соотношение соответствует полной длительности цикла, то меньшее выделяет одну аллометрическую фазу. В связи с этим весь диапазон синхронизации е3 распадается на аллометрическую фазу, в которой сохраняется тенденция убывания относительных приростов, и фазу перестройки — от момента завершения аллометрического участка развития до окончания периода полной синхронизации.
Оценим величину доли фазы перестройки в полном диапазоне синхронизации (е3—ее)/е3===24,6%. Эта величина достаточно близка к одной из границ, полученных для уравнения Гаусса (66).
Таким образом, большой цикл синхронизации можно представить состоящим из аллометрической фазы, занимающей 3/4 его длительности, и фазы перестройки, длительность которой составляет 74 цикла.
Оценим положение фазы перестройки в последовательности неравномерных критических рубежей
хк - 0,246z* — а** = 0,754 — 3,167* - 0,246.0,16Г* = 2,427*.
Следовательно, диапазон между последовательными неравномерно чередующимися рубежами хк и xk~l включает как границы рубежей равномерной последовательности TÄ, так и границу фазы перестройки.
В связи с тем, что в синхронизированный диапазон равномерно расположенных рубежей входят характеристики процессов равномерной шкалы до А:—3-го порядка, можно ожидать проявление между критическими неравномерными рубежами хк и хк~л границ равномерных рубежей по крайней мере ft—1-го порядка.
В соответствии с соотношением (68) границы равномерно расположенных рубежей связаны зависимостью (75). Тогда между границами неравномерного диапазона хк и хк~1 равномерно чередующиеся рубежи Тк~л зададут последовательность соотношений:
хк^+Тк-^уШк-1=1,53Т\    ^-1+4Г-1=7,16Г-1=2,63Г,
xk-i+2Th-i=5A6Tk-l=i190T\   a^i+5Tk-i=8916Th~i=3ßOTk.
Таким образом, на основе рассмотрения последовательностей равномерно и неравномерно чередующихся рубежей выявляется участок агрегирования рубежей, включающий основные границы, соответствующие неравномерным последовательным критиче
74

Таблица 3.  Положение рубежей
кя
Рубеж
Возраст в длительностях цикла Г

Рубеж
Возраст в длительностях цикла Т
1
Основной к-1

1,16
6
Перестройки
2,42
2
Равномерный к-
-1
1,53
7
Равномерный к-1
2,63
3
Равномерный к-
-1
1,90
8
Равномерный к-1
3,00
4
Равномерный к

2,16
9
Основной к
3,16
5
Равномерный к-
-1
2,27



ским значениям xh и xh~\ промежуточным равномерным рубежам, длительность которых составляет Th и Г*"1, а также рубеж фазы перестройки. Полная совокупность этих рубежей приведена в табл. 3.
Для определения положения основного рубежа xh и остальных, указанных в табл. 3 рубежей, необходимо знать длительность любого равномерного цикла Тк. По нему могут определяться как положения предыдущих и последующих границ по уравнению (73), так и длительности равномерных циклов различных порядков по зависимости (75).
Могут быть два основных подхода к определению величины Т\
1. Взять за исходное положение некоторого хорошо изученного рубежа и по нему по зависимости (73) определить Тк.
2. Принять в качестве исходной величины длительность некоторого хорошо изученного равномерного цикла, определяющего критические рубежи исследуемого порядка.
Таким образом, мы аналитически определили элементарное звено в цепи развития, в пределах которого происходит синхронизация как основных равномерных рубежей порядка к и к—1и неравномерно чередующихся рубежей. Определили и положение в этом звене времени завершения этапа перестройки, относящегося еще к предыдущему звену. Схема построения этого звена дана в табл. 3. Подобные типовые звенья, повторяясь, составляют целую цепь, отвечающую всему процессу развития рассматриваемой системы. В этой цепи каждое последующее звено связано с предыдущим соотношением (67).
75

4. КОЛИЧЕСТВЕННАЯ ОЦЕНКА УРОВНЕЙ БИОЛОГИЧЕСКОЙ ОРГАНИЗАЦИИ
Как отмечалось в разделе 1.1, в последнее время большое внимание уделялось классификации иерархически соподчиненных биологических структур. В связи с этим рассмотрим соотношение размеров биологических систем различных уровней.
В широко распространенном в наше время учебнике «Биология клетки» [де Робертис, 1973] приводится классификация разделов биологии, основанная на размерных диапазонах изучаемых объектов и способах их изучения (табл. 4).
Диапазоны, выделенные Робертисом и др., нанесены на рис. 26, где они сопоставляются с расчетными аллометрическими диапазонами, полученными в результате последовательных делений исходного размера, за который взят рост (длина тела) «условного мужчины» 170 см [Человек..., 1977], на ее:Х2=Х1/ее.
При этом нижняя критическая граница первого диапазона составит
1,70/15,15 см=11 см.
Следующие границы находятся аналогичным образом.
Как видно из рис. 26, диапазоны Робертиса и др. в большинстве случаев перекрывают два и более расчетных диапазонов. Это относится к «органам» (3,5 диапазона), клеткам и бактериям (1,5 диапазона), клеточным компонентам и вирусам (2 диапазона).
Следует заметить, что де Робертис и др. [1973] отмечает условность проведенных ими границ между различными структурами, так как эти границы установлены на основе разрешающей спо-
Таблица 4.  Разделы биологии [по: де Робертис и др., 1973]
Размеры объектов, м
Раздел биологии
Структура
Методы изучения
>ю-4
Анатомия
Органы
Глаз и простые линзы
10-4-10-5 10-5-2-10-7
Гистология Цитология
Ткани
Клетки и бактерии
Световые микроскопы, рентгеновский микроскоп
2-Ю-7-Ю-9
<ю-9
Субмикроскопическая  (ультраструктурная) биология Молекулярная биология
Клеточные компоненты, вирусы
Расположение атомов и молекул
Поляризационный и электронный микроскопы
Рентгеноструктур-ный анализ
76

Рйс. 26. Сопоставление диапазонов де Робертиса, расчетных аллометрических диапазонов и размеров структур человека <i — размерная шкала (м), б — диапазон де Робертиса и
ДР. [1973]. е — границы и порядковые номера расчетных диапазонов,
г — предельные размеры различных структур в организме человека
м
1Y
1<Г<
•10'
И
диапазоны де Робертиса
органы
ткана
клешни, бактерии
клеточные компоненты\ вирусы
атомы и молекулы
Расчетные диапазоны
Структуры человека
организм
органы
Эндокринные
железы
Y покровы
- гаметы ' клетки -I   ядра
органоиды
мембраны органические 1   молекулы
■атомы
г
собности используемых для исследований приборов, а не размерных пределов собственно биологических структур. Кроме того, привлечение весьма многообразного материала по структурам животных, растений и микроорганизмов дает довольно пеструю картину, которая не позволяет выделить четкие границы между ними.
Поэтому, учитывая, что у разных организмов могут быть свои, различающиеся по абсолютным величинам системы иерархии, мы привлекли для сопоставления с расчетными диапазонами данные о размерах структур одного, достаточно хорошо изученного вида, а именно человека4. Данные взяты главным образом из справочника «Человек. Медико-биологические данные» [1977], учебников гистологии [Елисеев и др., 1972] и цитологии [Трошин и др., 1970] и справочных изданий, содержащих количественные характеристики различных структур человека и животных [Блинков, Глезер, 1964; Автандилов, 1973; David, 1977]. Эти данные сведены в табл. бив обобщенном виде представлены на рис. 26.
Как следует из данных таблицы 5, размеры структур организма в значительной степени объединяются в определенные группы. Первая — это «основные» внутренние органы; в таблице они расположены по убывающим максимальным линейным размерам. Среди внутренних органов наибольшие средние размеры имеет желудок (0,37 м), минимальные — желчный пузырь (0,1 м)— рис. 27.
Затем идет вторая группа — эндокринные и некоторые другие железы,—начинающаяся с вилочковой железы (8,4 см) и закан-
При наличии раздельных данных, относящихся к мужчинам и женщинам, использованы данные, относящиеся к взрослым мужчинам. В некоторых случаях, при отсутствии данных по человеку, использованы материалы по другим млекопитающим, что оговаривается.

Таблица 5.  Размеры структур организма человека
Структура
Размеры

пределы
среднее
значение
(м)


1,7
1
29,5-49,5 см
0,37 0,20
2
27-30 см 20-30 см 16-20 см
0,28 0,25 0,18
3 4
16,5 см 14-18 см
0,16 0,16 0,23
5 6
10-14 см
9,7-14 см 10-12 см 8-12 см
0,12 0,12 0,11 0,10
8 7 9
10

0,084
11
5-8 см 3-7 см 4-6 см 4-5 см 1,7-4,7 см
0,06
0,05
0,05
0,045
0,03
12 13 14 15 16
10,5-17 мм 4-15 мм
0,014 0,01
17 18
5-9 мм
0,007
19
7-14 мм
ю-ю-3
20
2-12 мм
7-Ю-3
-
6-13 мм 0,5-4 мм
9-Ю-3 2,3-Ю-3
21 22
0,5-3,5 мм 1-2 мм
2-Ю-3 1,5-Ю-3
23
. 24
0,5-2,5 мм
1,5-Ю-3
25
1,3-ю-3
26
Источник
Условный мужчина
Органы Желудок
Легкие Печень
Мозг
Поджелудочная же-
Селезенка
Сердце
Почки
Желчный пузырь Эндокринные железы Вилочковая железа (3-14 лет)
Щитовидная железа
Надпочечники
Околоушная железа
Семенники
Предстательная
железа
Гипофиз
Паращитовидная
железа
Эпифиз
Покровы и стенки
(толщина)
Стенки желудочков
сердца
Стенки желудка Кожа с подкожной клетчаткой Стенки предсердий Стенки желчного пузыря
Слизистая оболочка желудка Кожа
Человек, 1977, с. 28
Человек, 1977, с. 140 Автандилов, 1973, с. 110
Человек, 1977, с. 165 Человек, 1977, с. 153 Автандилов, 1973, с. НО
Человек, 1977, с. 215
Человек, 1977, с. 157
Автандилов, 1973, с. НО
Человек, 1977, с. 114 Человек, 1977, с. 125 Человек, 1977, с. 178 Человек, 1977, с. 155
Человек, 1973, с. 114
Человек, 1977, с. 210 Человек, 1977, с. 206 Человек, 1977, с. 134 Человек, 1977, с. 186 Человек, 1977, с. 188
Человек, 1977, с. 210 Человек, 1977, с. 204
Человек, 1977, с. 208
Человек, 1977, с. 126
Автандилов, 1973, с. 109
Человек, 1977, с. 141
Елисеев и др., 1972, с. 402
Человек, 1977, с. 126
Елисеев и др., 1972, с. 155
Человек, 1977, с. 141
Человек, 1977, с. 65
78

(продолжение)

Размеры


Структура

среднее

Источник 
пределы
значение    
(м)


в т. ч. эпидермис
46-123 мкм
8,5-Ю-5

Человек, 1977, с. 59
Серозная оболочка
0,7—1,1 мм
0,9-10-3
27
Елисеев и др., 1972, с. 533
брюшной стенки




Половые клетки



Человек, 1977, с. 19
Овоцит Яйцеклетка
130-150 мкм
14- Ю-5
28  
89-91 мкм
9-Ю"5
29
Человек, 1977, с. 192
Сперматоцит
58-67 мкм
6,3-ю-5 5,45-Ю-5 4-Ю-5
30 31 32
Человек, 1977, с. 20 Человек, 1977, с. 187
Сперматозоид
  
Елисеев и др., 1972,
Ядра яйцеклеток
  
с. 50
Тканевые клетки
5-40 мкм
-
-
Свенсон, Уэбстер, 1980, с. 142
и клетки крови



Блинков, Глезер, 1964, с. 382 


 
Корковые пирами-
12-36 мкм
24-Ю-6
33
Блинков, Глезер, 1964, с. 383
дальные клетки Ганглиозные клетки
8-40 мкм
24-Ю-6
34
Блинков, Глезер, 1964, с. 382
головного мозга



Елисеев и др., 1972, с. 166
Тучные клетки
-
22-Ю-6
35 
Печеночные клетки
20-25 мкм
22-Ю-6
36
Елисеев и др., 1972, с. 513
Фибробласты
-
20-Ю-6
37
Елисеев и др., 1972, с. 161
Клетки крови
5-11 мкм

38
Елисеев и др., 1972,  
 
с. 131—141
Ядра




Печеночных клеток
7-16 мкм
11-10-*
39
Елисеев и др., 1972, с. 513
Почечных, печеноч-
6-14 мкм
ю-ю-6
40
Davis, 1977, с. 28, 30
ных и других клеток



Елисеев и др., 1972, с. 50
Сперматозоидов
4-Ю-6
41 
Органоиды и другие




ультраструктуры



David, 1977, с. 40
Ядрышки печеноч-
1,2-2,1 мкм
1,7-10-6
42 
ных клеток




Митохондрии пече-
0,8-2 мкм
1,4-Ю-6
43
Елисеев и др., 1972, с. 514
ночных клеток




Центриоли (длина)
0,3-2 мкм
1,2-Ю-6
44
Робертис де и др., 1973, с. 413
Лизосомы
_
МО"6
45
Трошин и др., 1970
Микротельца пече-
0,4-0,8 мкм
6-Ю-7
46
Грубан, Рехцигл, 1972, с. 27
ночных клеток


 
Базальные тельца
1,5-Ю-7
47
Робертис, 1973, с. 413
(кинетосомы)


 
Микроворсинки
0,05-0,1 мкм
7,510-8
48
Трошин и др., 1970, с. 86
(диаметр)
1

 
79

(окончание)

Размеры


Структура

среднее
значение
(м)
Источник 
пределы 


Рибосомы
15-35 нм
2,5-10~8
49
Трошин и др., 1970, с. 86
Микротрубочки
20-25 нм
2,2-10-8
50
Браун, Рахишев, 1975
(диаметр) Окаймленные
15-20 нм
1,7-Ю-8
51
Неворотин, 1977
пузырьки
10-15 нм
1,2-Ю-8
52
Трошин и др., 1970,
Десмосомы (межкле-   
с. 87
точные контакты) ^




Мембраны в клетках




(толщина)
Мембрана эритроци-
8,5-9,25 нм
8,9-Ю-9
53
David, 1977, с. 136
тов
Мембрана клеток
7,67-9,67 нм
8,7-10-9
54
David, 1977, с. 134
эпителия ротовой




полости
Ядерная оболочка
7-8нм(Х2)
7;5-10-9
55
Елисеев и др. 1972, с. 50
Мембрана аппарата
6-7нм
6,5-Ю-9
56
Елисеев и др. 1972» с. 43
Гольджи

6,0-ю-9
57
Елисеев и др. 1972,
Мембрана митохонд-
  
с. 41
рий




Молекулы Альбумин Глобулин Актин
РНК, ДНК (шаг
5-15 нм 4-12 нм
3,0-3,4 нм
10-Ю-9 8-Ю-9 5,5-Ю-9 3,4-Ю-9
58 59 60 61
David, 1977, с. 145
David, 1977, с. 145
Финеан, 1970, с. 128
Финеан, 1970, с. 192, 187
спирали)
2-3 нм
2,5-Ю-9 2,0-Ю-9
62 63
Финеан, 1970, с. 126
Миозин (диаметр) ДНК (диаметр спи-   
Финеан, 1970, с. 18в
рали)
Холестерин Полипептидная цепь,
0,6-1,8 нм 0,98 нм
1,8-Ю-9 9,8-Ю-10
64 65
Финеан, 1970, с. 244 Финеан, 1970, с. 119
диаметр Целлюлоза
0,476X1,828 нм
4,7-Ю-10
66
Финеан, 1970, с. 29&
Атомы и расстояния




между ними

3,88-Ю-10
67
Пиментел, Спартли,
Атом Na, диаметр
  
1973, с. 56
Атом Н, диаметр

1,42-Ю-10
68
Пиментел, Спартли, 
  
1973, с. 56
Ионные связи
0,2-0,3 нм
2,5-Ю-10
69
Робертис де и др., 1973, с. 58
Межатомные рассто-
0,096-0,269 нм
1,8-Ю-10
70
Коптев, Пентин, 1977, с. 22
яния
Расстояния между
0,98-10-10
70
Маррел и др., 1968, с. 83
Н и О в молекуле воды




80

чивающаяся эпифизом (7 мм). Нижние границы средних размеров органов и эндокринных желез практически совпадают с нижними границами первого и второго расчетных диапазонов (11 см и 7,4 мм).
Следующие размеры, обнаруженные в литературе,— толщина различных структур стенок, покровов, органов, тканей. Здесь резко выделяются среди прочих стенки желудочков сердца (14— 7 мм) и желудка (13—6 мм). Далее идет кожа с подкожной клетчаткой (2,3 мм) и т. д. Определение линейных размеров собственно гистологических структур затруднительно, так как «в длину» некоторые из них распространены вдоль всего тела, а в толщину они сильно варьируют, в нижнем пределе достигая размера клетки, например в однослойном эпителии.
В следующий — четвертый расчетный диапазон (4,9 • 10~4—3,2 • •10"5 м) четко вписываются половые клетки, как зрелые, так и гаметоциты. Здесь же располагаются ядра яйцеклеток (4-10~5м).
Подавляющее большинство тканевых клеток, включая клетки крови, имеют размеры 40—5мкм [Блинков, Глезер, 1964; де Робертис и др., 1973; Свенсон, Уэбстер, 1980; и др. J и приближенно соответствуют пятому расчетному диапазону (32—2,1 мкм). Сюда же относятся и ядра клеток (рис. 28).
Органоиды цитоплазмы и ядра — митохондрии, лизосомы, ядрышки, более мелкие ультраструктуры (рибосомы, десмосомы), а также недавно описанные окаймленные пузырьки занимают шестой и седьмой диапазоны (2,1 -10~6—1,4-10"7—9,3-10~9 м) V
Некоторые ультраструктуры в зависимости от функциональной активности в значительной степени меняют свои размеры. Это относится к аппарату Гольджи, различным вакуолям. Пиноцитоз-ные пузырьки, например, варьируют от 5 нм до нескольких мкм [Трошин и др., 1970]. Кроме того, мелкие субмикроскопические структуры изучены еще недостаточно. Поэтому более детальное рассмотрение ультраструктур представляется в настоящее время затруднительным.
Далее следует группа различных клеточных мембран, толщина которых колеблется в-очень узких пределах (9-10~9—6-10"* м). Близкие размеры имеют глобулярные белки, диаметры молекул нуклеиновых кислот, фибриллярных белков, липиды. Молекулы ДНК и белков имеют значительно большие размеры. Мембраны и органические молекулы занимают преимущественно   восьмой
Размеры элементов полипептидных цепей кератинов характе ризуются величинами, показанными на рис. 29 [по: Финеан, 1970], из которого видно, что имеется разрыв в диапазоне 9-10""11— 5* • 10~10 м, что соответствует рассчитанной критической границе.
Еще мельче атомы, размеры и расстояния между которыми соответствуют девятому диапазону (6,1 • 10""10—4,0 • 10~и м).
81

2 4 8 16 Порядковый номер
10~\ м
20
+
Рис. 27. Вариационный ряд длин основных органов человека
По абсциссе — порядковые номера, по ординате — длина. Шкала логарифмическая. Органы расположены в порядке убывания максимального линейного размера: желудок, легкие, печень, мозг, поджелудочная железа, селезенка, сердце, почки, желчный пузырь, вилочко-вая, щитовидная железа, надпочечники, слюнная железа, семенники, предстательная железа, гипофиз, па-ращитовидная железа, эпифиз [По: Человек, 1977]
Рис. 28. Вариационный ряд диаметров ядер клеток разных органов человека
По абсциссе — порядковые номера; по ординате — диаметры ядер. Шкала логарифмическая. Границы сверху и снизу — расчетные [Диаметры рассчитаны по данным Хесина, 1967]
Рис. 29. Вариационный ряд периодов а- и ß-ке-ратинов при больших углах рассеяния По абсциссе — порядковый номер;  по ординате — период на рентгенограмме [По: Финеан, 1970]


+++*
т+



1
г» \



1     +
'       I    N111   I
4 8 16 30 Порядковый номер



20


10


6


4

2


1

,         ,     I I | 11 1 1 1 I I—. 

2       4           Ю       20


Порядкодый номер
Проведенное сопоставление линейных размеров различных групп структур организма человека с рассчитанными критическими диапазонами показывает весьма существенное их соответствие6  Это позволяет предположить, что критические диапазоны
82

Таблица 6. Расчетные диапазоны и характерные для них структуры организма человека
Пределы (м)
расчетного диапазона
верхний
нижний
Характерные структуры
1
2 3 4 5 6 7 8 9
1,7* 0,11 7,4-1 4,9-1 3,2-1 2,1-1 1,4-1 9,3-4 6,1-1
0,11 IM 4,9-1 3,2-\
1,4-; 9,3-: 6,1-4,0-
основные органы
ю-3
10-4
ю-5
Ю-6
ю-3 ю-4 ю-5
ю-6 ю-7 ю-9
.10-ю
эндокринные железы покровы, стенки
гаметы
клетки, ядра
органоиды и другие ультраструктуры
ю-7 ю-9
.10-ю
»
мембраны, молекулы
атомы
* Значение роста «условного мужчины» [Человек, 1977], принятое за исходное для расчета.
В тех же случаях, когда одного диапазона оказывается недостаточно, происходит изменение зависимости, которое на логарифмическом графике выражается изломом прямой (см. рис. 27) или (и) резким отскоком значений величин в последовательном ряду (см. рис. 29).
83